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VISCOSITY  OF  m  D  WORKING  FLUIDS 


M -Nab 


21 .1  GENERAL  INTRODUCTION 

Considerable  and  increasing  use  of  alkali  metals  is  currently  being  made  in 

devices  concerned  with  the  generation  of  electrical  power  by  unconventional  means 

(thermionic  diodes  and  magnetoplasmadynamic  generators  (MPD)"’  )  and  electrostatic 

propulsion*4-.  The  attractiveness  of  the  alkali  metals  follows  from  their  low 

ionization  potentials  and  the  consequent  ease  with  which  they  may  be  ionized. 

Few  attempts  have  been  made  to  calculate  or  measure  their  thermodynamic  and  transport 

properties  (other  than  electrical  conductivity)  and  the  effect  of  ionization  on 

5 

such  properties  . 

The  element  cesium  is  of  particular  interest  in  closed  cycle  MPB  electrical 
power  generation  since,  although  it  has  a  high  capital  cost,  it  has  the  lowest 
ionization  potential  of  the  group.  Base-lead  MPD-steam  stations  envisage  the 
use  of  a  3mall  atomic  percentage  of  alkali  metal  in  an  inert  gas.  such  as  helium 
or  argon,  since  it  is  with  seeding  f ractions  of  abcut  one  atomic  percent  that  the 
maximum  electrical  conductivity  is  obtained.  For  space  powe  •  generation,  radiator 
size  and  weight  considerations  may  necessitate  the  use  of  a  Rinkine  cycle  with 
a  pure  alkali  metal. 

In  the  present  calculations,  the  viscosity  of  pure  cesium  has  been  evaluated 

theoretically  and  the  part  played  by  ionisation  assessed.  (Under  given 

temperature  and  pressure  conditions  ionization  effects  are  most  noticeable  in 

cesium  because  of  the  low  ionization  potential.)  The  viscosity  is  of  particular 

interest  at  the  present  time  since  consideration  of  'second  order*  effects  in 

MPB  generators  includes  a  study  of  the  boundary  layers  in  such  generators. 

6-8 

Recent  studies  indicate  that  the  power  loss  arising  from  re  erse  currents  in 
the  boundary  layers  may  reduce  the  power  output  of  an  MPD  generator  by  a  factor 
of  four. 

Following  this  calculation,  which  is  appropriate  to  a  kine  cycle  MPD 
generator  using  pure  cesium  as  the  working  fluid,  the  viscosily  cf  cesium-helium 


21 .1 


mixtures  has  been  investigated,  the  effects  of  seeding  fraction  being  of  particular 
interest  in  this  case  which  applies  to  terrestrial  base=-load  power  generation. 

The  general  behaviour  of  the  viscosity  of  binary  mixtures  is  then  examined, 
special  regard  being  paid  to  the  occurrence  of  maxima  in  the  viscosity  aiid  their 
dependence  on  temperature  for  a  large  number  of  gases. 

21.2*  VISCOSITY  OF  PARTIALLY -IONIZED  GASEOUS  CESIUM 

21.2.1  INTRODUCTION 

To  calculate  the  viscosity  of  a  gas,  it  is  necessary  to  determine  the 
velocity  distribution  function  f  (r,  v,  t)  for  each  species  of  particle  present. 
This  function  is  defined  such  that 

/j(r,  v,  t)  dx  dy  dj  dvx  dvy  d.va  dt  .....  (21.1 ) 

is  the  number  of  particles  of  species  j  which  nave  position  co-ordinates  between 

a-  and  x  +  da,  y  and  y  +  dy,  $  and  $  +  d$,  and  which  have  components  of  velocity 

between  v*  and  v*  +  dv« ,  Vy  and  Vy  +  dvy ,  and  Vj  +  dv5 .  If  the  dependence 

of  f  upon  the  variation  of  flow  velocity  with  position  in  the  gas  is  known, 

the  viscosity  can  be  calculated.  Unfortunately,  the  Boltzmann  -  Fckker  -  Planck 

equation,  which  according  to  the  statistical  mechanics  of  ionized  gases  gives  f  to 

9  10 

a  good  approximation  for  charged  particles  *  ,  is  mathematically  intractable. 

Data  concerning  the  velocity  distribution  function  can  be  obtained  from  the  simpler 

Boltzmann  equation,  although  the  use  of  this  equation  involves  the  assumption 

that  the  time  during  which  two  particles  are  in  contact  (this  is,  a  collision) 

is  infinitesimal  compared  with  the  time  between  collisions;  this  assumption  is 

9  11 

doubtful  when  charged  particles  are  present  *  .  The  use  of  an  approximate 

12  13 

solution  of  the  Boltzmann  equation  developed  by  Chapman  and  Enskog  is  discussed 
in  section  21 .2.4-. 

The  only  other  feasible  approach  seems  to  be  that  of  elementary  kinetic 
theory,  in  which  viscosity  is  treated  as  a  'free  path  phenomenon'.  The 
coefficient  of  viscosity  can  be  calculated  very  roughly  by  assuming  a  simple  form 
for  the  velocity  distribution  function,  and  evaluating  the  momentum  transported 
across  unit  area  in  the  ga3  by  the  thermal  motions  of  the  particles.  This  method 
is  developed  in  section  21.2.3  to  determine  the  nature  of  the  effect  of 
ionization  on  viscosity;  it  is  too  approximate  to  yield  a  great  deal  of  useful 
information. 

*  Based  on  ED  63-10  'The  viscosity  of  gaseous  cesium  at  temperatures  up  to 
3000°K*  C.A.  Robinson  and  I.R.  McNab,  and,  'Viscosity  of  partially  ionized  gaseous 
cesium'  CJU  Robinson  and  I.R.  McNab.  J.App.Phy3,  June  1964. 
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The  viscosity  of  an  ionised  gas 

Viscosity,  the  development  of  shear  stresses  in  a 
the  transport  of  momentum  from  point  to  point  by  the  ran 
the  molecules.  In  simple  kinetic  theory {.  iscosity  is 
phenomenon’ ,  since  the  free  path  is  the  distance  a  parti 
collisions  which  change  its  momentum.  In  more  rigorous 
with  the  free  path  becomes  more  obscure. 


moving  fluid*  is  due  to 
dom  thermal  motions  of 
known  as  a  'free  path 
cle  travels  between 
theory  the  connection 


A  monatomic  gas  in  which  a  small  percentage  of  the  atoms  are  thermally 
ionized  will  contain  not  only  charged  particles  and  neutral  atoms  in  the  ground 
state,  but  also  a  considerable  number  of  excited  particles.  with  the  result  that 
inelastic  ar  superelastic  collisions  will  c .  as lonai^y  taxe  pj-flce*  xne  e*i  e^t 
of  these  excited  particles  is  very  complicated  to  analyse  and  has  been  neglected 
in  this  treatment.  It  also  appears  from  considerations  of  elementary  kinetic 
theory  that  the  momentum  transported  by  ele  irons  Ln  a  1  per  cent  ionised  gas 
will  be  most  10“^  of  that  carried  by  neutrals.  owing  +o  the  small  mass  and 
large  charge/mass  ratio  of  the  electron.  A  eordingiy,  the  momentum  transport 
due  to  electrons  has  been  neglected,  as  has  the  small  effect  of  collisions  with 


electrons  on  the  motions  of  the  heavy  parti  s .=  The  plasma  is  thus  reduced,  as 
far  as  viscosity  phenomena  are  concerned,  to  a  finery  mixture  of  positive  ions  and 
neutral  atoms,  with  electrons  present  omy  to  restore  the  macroscopic  charge 
neutrality.  The  resultant  simplification  ;f  the  ths-ry  is  cnsiderable. 


21.2.3  Approximate  calculation  of  viaCySity  from  elementary  kinetic  theory 


An  expression  for  the  viscosity  of  a  sligntiy  ionized  gas  will  be  obtained 
from  'free  path*  considerations.  Collisions  between  atoms  and  atoms,  and  between 
atoms  and  ions  will  be  assumed  to  behave  as  encounters  between  rigid  elastic 
spheres  in  whirh  the  distances  of  closest  approach  of  the  centres  are  respectively 
p  and  p ' .  The  mean  free  path  of  an  atom  is  therefore 
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/2  n  (nap^  + 


n  c 

r 


.....  (21.2) 


The  force  between  two  ions  will  be  assumed  to  be  tnat  of  Coulomb  repulsion 
only,  since  it  is  unlikely  that  two  ions  in  a  plasma  at  the  temperatures  under 
consideration  ( 1 00C°  to  2b00°K)  will  have  sufficient  kinetic  energy  to  approach 
closely  enough  for  the  short-range  quantum  mechanical  forces  to  be  important. 

The  evaluation  of  the  mean  free  path  of  an  ion  is  complicated  by  the  fact  that 
any  one  ion  may  at  any  instant  be  ciose  enough  tc  each  of  several  others  to 
interact  with  them,  since  the  mean  distance  between  ions  may  be  considerably  less 
than  the  Debye  shielding  distance: 


h  = 


_  l  kT 


1/2 


4 tr«  e 

k  9 


(21 .3) 


The  mean  distance  an  ion  travels  before  its  velocity  is  deflected  through 

approximately  45°  may  be  evaluated.  This  will  be  used  ae  the  mean  free  path 

of  an  ion  among  ions.  The  mean  square  change  per  second  in  the  velocity  of  an 

ion  in  $  direction  perpendicular  to  its  initial  velocity  is  denoted  by  <(AVjJ2)  . 

The  time  after  which  an  ion  has  been  deflected  appreciably  from  its  original 

direction  by  other  ions,  called  by  Spitzer*^  the  deflection  time,  is 

v2 

=  y  ,  x2v  *****  (21‘^) 

<  (AvJ  > 

During  this  time  the  ion  will  have  travelled  a  distance  which  will  be  used  as 
the  mean  free  path  of  an  ion  among  ions: 


<(Av* X) 

Employing  methods  similar  to  those  used  by  Rosenbluth  et  al^ 

<uvj2>  =  “y/z  in4 


(21.5) 
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where  A  = 
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The  expression  (21.5)  now  becomes 

mi2vV(v2+?) 

”  4 

otr  n^  e  In  A 

which  substituting  root  mean  square  values  becomes  approximately 
9f2  k2T2 

„  4 
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When  the  effect  of  collisions  with  atoms  is  taken  into  consideration,  the  mean 
free  path  of  an  ion  in  the  plasma, 
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is  ob  lined. 

The  coefficient  of  viscosity  due  to  each  species  of  particle  is  the  momentum 
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transported  by  particles  of  that  species  in  unit  time  across  unit  area  in  the  gas 
per  unit  velocity  gradient  perpendicular  to  that  area*  If  the  velocity 
distribution  is  assumed  everywhere  locally  Maxwellian ,  the  approximate  expression 
for  the  coefficient  of  viscosity  is 


n 


M 


2m  kT 

!  V~ 


.....  (21  .11) 


Summing  over  the  two  types  of  particle 

(naAa  Vma  +  n,k.  Vm, 
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Values  for  the  mean  free  paths  may  now  be  substituted.  Assuming  ma  ■  m, 
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When  the  fractional  ionization  a  is  small,  by  far  the  greater  part  of  the 
viscosity  is  due  to  the  first  term  in  the  bracket.  The  Couiomb  interactions 
enter  the  expression  only  through  the  term  in  in  A.  and  their  only  effect  is  to 
reduce  further  the  already  small  second  term  in  the  bracket.  It  may  be  inferred 
that  the  viscosity  of  a  slightly  ionized  gas  is  almost  entirely  due  to  momentum 
transported  by  neutral  atoms,  and  that  the  effect  of  ionization  is  slightly  to 
reduce  the  viscosity,  through  the  effects  of  Coulomb  repulsion.  The  physical 
interpretation  is  that  the  long-range  Coulomb  forces  reduce  the  momentum 
transported  by  ions  by  shortening  the  mean  free  path  of  the  ions. 

21.2.4  Theory  based  upon  the  Boltzmann  equation 
21.2.4.1  General  theory 

As  stated  in  section  21.2.2,  elementary  kinetic  theory'  calculations  show 
that  virtually  all  the  viscosity  of  a  slightly  ionized  monatomic  gas  will  be 
due  to  momentum  transport  by  atoms  and  ions;  this  results  in  considerable 
simplification  of  the  calculation. 

To  compute  the  viscosity  of  a  gas  it  is  necessary  to  obtain  information 

about  the  velocity  distribution  function  of  each  species  of  particle  present. 

This  is  usually  done  'y  solving  the  Boltzmann  equation  by  the  approximate  method 

12  13 

due  to  Chapman  and  Enskog  ,  which  has  be*  highly  developed  by  Chapman  and 

1 1 

Cowling  ,  The  distribution  function  is  expanded  as  a  series  and  a  set  of 
integral  equations  derived  ''or  coefficients  describing  the  first  order 
perturbation  from  the  known  equilibrium  solution.  These  equations  are  solved 
by  an  expansion  in  orthogonal  polynomials.  Expressions  are  obtained  for  the 
transport  coefficients,  including  viscosity  of  a  pure  gas  or  gas  mixture,  in 


terms  of  the  collision  integrals,  which  depend  on  the  interparticle  forces  through 
the  angle  of  deflection  x  produced  when  two  particles  collide  with  relative 
velocity  g  and  impact  parameter  b.  The  results  are  expressed  most  conveniently 
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the  colliding  particles,  and  the  Debye  length  is  the  characteristic  Coulomb 
interaction  distance,  given  by  equation  (21 *3) • 

The  collision  integrals  have  been  normalized  in  several  ways,  the  notation 
followed  here  being  that  of  Kirschf elder,  Curtiss  and  Bird  ,  thus: 
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exp(-f^-jg*2s+3  ( *  “Cos  py)b*  db*  dg*  ....(21.14) 


where  the  deflection  angle  %  in  a  centre -of -mass  frame  of  reference  is  given 
by  classical  mechanics: 


X  =  -  2b*  g* 

in  which  u  =  Vr. 
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min 
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The  expression  for  the  first  approximation  to  the  viscosity  of  a  gas 


proximations  require 


(l  1  ( 2  2 V 

mixture  involves  only  Cr  *  y  andfr  ,  highs 

( r  a )  * 

evaluation  of  fM  *  '  for  larger  p  and  s.  The  viscosity  of  a  pure  gas  (to  the 
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first  approximation)  is 

t)  =  266.93  1(f7  (ifT)*  /fr2n^>2^  (T*)  .....  (21.16) 

where  r?  is  in  poise,  a  in  £  and  T  in  °K.  The  viscosity  of  a  binary  mixture 
of  components  of  the  same  molecular  weight  M  (such  a3  the  mixture  of  atoms  and 
ions  here)  is  given  in  terms  of  the  viscosities  nA  and  q  of  the  separate 

*  n 

components  aid  a  third  quantity,  jj  g,  expressing  the  atom-ion  contribution: 


P12 


"1  „  2n  {  J  ,2)* 


=  266.53  10  (Itr/  /  f12  fi 


(Ti2*) 
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c.  ! 
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in  which  ,  =  (c^+c^/Z,  The  viscosity  of  a  binary  mixture  of  atoms  and  i^ns 
with  fractional  ionization  a  then  becomes 


b  1  +  Z 
?mix  X  +  Y 


* •  <  . .  (21 .18) 


where 


+  Ml 3Lt  +  3l 


Y  -  ^  A*  ^ 
5  ®i  n 


.  \2  2a(l^*)n  .  2 
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These  formulae  may  be  used  to  evaluate  the  viscosity  when  the  appropriate  collision 
integrals  have  been  calculated  from  the  interparticle  potential  functions. 

21.2.4.2  Interparticle  potential  functions  and  collision  integrals 

In  the  absence  of  experimental  data  for  the  interaction  potentials  in 
cesium  the  following  potential  functions  have  been  assumed. 

Atom-atom 


The  Lennard-Jones  (6-12)  potential 

r/o  \  12 


■(=?)]• 


(21.21) 


which  takes  into  account  the  inverse  seventh  power  attractive  (van  der  Waal’s) 

forces  and  short  range  repulsive  forces,  has  been  used  with  considerable  success 

in  the  past  (see,  for  example,  Ref.  16)  and  will  be  assumed  here.  The  variation 

(2  2)* 

of  6  (r)  with  r  is  shown  in  Fig.  21 .1 ;  the  colision  integral  fl  1  for  this 
a 


potential  is  tabulated  in  Ref.  1b  and  shown  in  Fig.  21.4. 

Atom-ion 

For  this  case  the  potential  requires  '\n  inverse  fourth  power  attractive 
component  to  account  for  the  attraction  between  the  ionic  charge  and  the  induced 
dipole  moment  in  the  atom,  while  the  short  range  repulsive  forces,  which  vary 
rapidly  with  distance,  may  conveniently  be  represented  by  an  impenetrable  core. 
The  potential  takes  the  form 


<t  .  (r)  =  00  for  (r  <  a  ,  ) 
rai  ai 


and  <£ai(r)  = 


ai  r 


c  .  \  4 

“TV  for  ( 


.....  (21.22) 


r  >  «t  .)  . 

ai 
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A  graph  cf  d>  ,(r)  against  r  is  shown  in  Fig.  21.2, 
ax 

Writing  equation  (21.15)  in  the  form 

X  *  *  -26  .  (21  .23) 

and  substituting  the  above  potential,  yields  an  elliptic  integral  of  the  first 

kind  for  0.  Thi3  can  be  reduced  to  standard  form  as  shown  in  Appendix  21 .A , 

and  evaluated  by  references  to  tables  or  numerical  integration.  ^iile  no 

definite  information  is  available  for  cesium, experimental  polarizabilities  for 

other  monatomic  gas  molecules  suggest  that  the  interaction  potential  e^/k  will 

be  of  the  order  of  500°K,  Since  the  range  of  gas  temperature  under  consideration 

is  1000°  to  3500 °K  all  cases  of  practical  interest  here  will  be  covered  by  a  range 

of  values  of  TV  from  1  to  20.  The  required  integrals  ft  .  *  '  and  ft  '  *  y 
Si  &1  fil 

(Iq.  21 ,14)  have  been  evaluated  by  mechanical  quadrature  on  a  Pegasus  digital 
computer  for  several  values  in  this  range,  and  are  shown  in  Table  21.1  and 
Fig.  21 .5*  This  method  of  evaluation  appears  superior  in  this  case  to  the 
approximate  method  suggested  by  Chapman  and  Cowling  (see  Ref.  11 )  which  is  mainly 
useful  for  potential  functions  of  this  form  when  the  attractive  component  is  very 
weak. 

Ion-ion 

The  most  precise  form  of  potential  to  represent  the  ionic  repulsion  is  the 
shielded  Coulomb  potential  predicted  by  the  Dabye-Huckei  theory: 

-1  -r/h 


<P  ^  oc  r  e 


-18 


However  it  has  been  shown  by  Liboff  that,  even  in  the  case  of  a  highly  ionized 
gas,  only  a  very  small  error  in  viscosity  computations  is  incurred  by  using  an 
ordinary  Coulomb  potential  with  a  cutoff  at  the  Debye  length  instead.  The  use 
of  this  modification  greatly  simplifies  calculation.  The  short  range  forces 
may  again  be  represented  by  an  impenetrable  core,  although  it  is  very  unlikely 
that,  at  the  temperatures  under  consideration,  the  colliding  pair  of  ions  will 
have  sufficient  kinetic  energy  for  this  to  be  brought  into  action.  The 
potential  function  then  assumes  the  form  (shown  in  Fig.  21 *3 ) i 


^(r)  =  oo  for  (r  <  a/) , 

h(r^  =  §7  (^)  for  < r  <  h)> 


and  <£^(r)  =  0  for  (r  >  h). 

Neglecting  the  few  pairs  of  ions  which  approach  as  closely  as  the  impenetrable 
core,  the  deflection  ‘.tngle  becomes  that  given  by  classical  mechanics  for  a 
Coulomb  potential: 
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X  =  2  arc  sin  (4  ^'*b*‘~  +  1  j"*7 

This  will  be  used,  for  all  ion-ion  encounters  in  which  the  impact  parameter  is  less 
than  the  Debye  length;  when  the  impact  parameter  is  greater  than  h  the  deflection 
is  zero.  Using  this  value  of  x  in  equation  (21 =14)  and  replacing  the  upper 
limit  of  integration  over  b*  of  oo  by  h*  yiexds 


0(2»2)*  _  1_  j  exv 
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Hirschf elder,  Curtiss  and  Bird  recommend  that,  since  the  term  in  square  brackets 

2 

varies  comparatively  slowly  with  g*  it  may  be  taken  outside  the  integral  sign,  g* 
being  replaced  everywhere  by  mean  value  over  all  collisions,  2Th*»  If  this  be 
done  the  expression  integrates  immediately  to  give 


n42-2)' 

ii 


1 


VP. 


,  ,  ln(l  +  16T  *2h*2)  - — 

L  x 
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(21  .24) 


The  slight  inaccuracy  introduced  through  the  last  approximation  is  not  important. 

(2  2)* 

An  accurate  evaluation  of  Q, .  *  should  consider  the  variation  of  h*  with 

XX 

temperature  and  ion  density,  approximate  values  of  h*  of  practical  interest  here 

p  ii. 

being  in  the  range  10  t  h*  t  10.  However,  owing  to  the  form  of  equations 

(21.18),  (21.19)  and  (21.20),  UmiYdoes  not  depend  critically  on  the  value  of 

n. provided  this  is  large  and  the  fractional  ionization  is  *<  0.1,  and  thus  h* 

may  be  assumed  constant  without  appreciable  loss  of  accuracy;  this  has  been 

-  (2  2)* 

confirmed  by  numerical  calculations.  Fig.  21.6  shows  the  variation  of  * 

with  as  given  by  equation  (21.24)  for  h*  =  10“,  10^  and  lO"4". 


21.2.4.3  Results  and  discussion 

(op)*  ill)*  (22)* 

From  the  values  of  Q  ’ *  given  in  Ref.  1b,  of  fi  ,  *  and  0  ,  * 

aa  al  (3 

computed  by  numerical  integration  and  shown  in  Table  21  .1  ,  and  0^  ’ 

calculated  from  equation  (21.24),  the  viscosity  coefficient  of  the  slightly 

ionized  cesium  may  be  calculated  for  any  set  of  values  of  the  force  constants, 
fractional  ionization,  Debye  length  and  temperature. 


The  values  which  have  been  used  here  are 


0a  “  rai 


=  0  =  3.5 X,  f  Jk  =  1 00°K,  <  ./k  =  500°K  and  „ ,/k  =  c2/k(T,  =  4.77.1 04  °K 
1  a  ax  1 


The  fractional  ionization  (a),  which  depends  cn  temperature  and  pressure,  is 
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obtained  from  Saha’s  equation  : 
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log,  (  H  PV  "  +  I  log  T  -  6.483,  . (21  .25) 

\l-q  / 

using  the  ionization  potential  of  3*893  V  for  cesium.  A  constant  reduced  Debye 
length  (h*5  of  10^  has  been  used. 

Fig.  21.7  shows  the  resulting  viscosity  of  ionized  cesium  (rj  ,  )  for  pressures 
o  i  i  c  mix 

of  10  ,10  ,  10*^  and  10  3  alai  and  temperature  in  the  range  1000°  to  3500°K  as 

found  from  equations  (21 .18)  (21.19)  (21.20),  and  the  viscosity  of  un-ionized 
cesium  (rj  )  over  the  same  temperature  range,  from  equation  (21.16).  It  is  clear 

a 

that,  over  the  greater  part  of  the  temperature  range,  the  viscosity  increases  with 
temperature;  the  viscosity  of  un-ionized  cesium  may  be  expressed  as 

r\  *  9*5  10~^  T^^  g/crn  sec 

It  is  also  clear  that  the  effects  of  ionization  are  small  until  the  fractional 
ionization  is  greater  than  1 at  fractional  ionizations  of  about  1C0S  the  reduction 
in  viscosity  due  to  ionization  is  over  20$.  Thi3  drop  in  viscosity  with 
increasing  ionization  may  be  explained  by  saying  that  the  ions  have,  owing  to 
their  mutual  Coulomb  interactions,  a  much  shorter  mean  free  path  than  neutral 
atoms,  and  hence  contribute  little  to  the  momentum  transport  process.  It  is 
doubtful  whether  the  Enskog-Chapman  solution  of  the  Boltzmann  equation  can  be 
used  to  evaluate  transport  coefficients  for  ionized  gas  having  a  fractional 
ionisation  greater  than  about  10??:  solution  of  the  Boltzmann-Fckker-Planch 
equation  becomes  the  appropriate  procedure^. 

The  neglect  of  electrons,  diatomic  cesium  and  excited  atomic  states  in  the 
partially  ionized  gas  introduces  negligible  error  in  the  computed  viscosity,  as 
do  the  approximations  made  in  deriving  the  ion-ion  collision  integral  (equation 
21 *24)  and  the  assumption  of  constant  h*  in  the  numerical  evaluation.  The  most 
likely  sources  of  error  are  the  critical  interaction  energies  (f),  effective  atomic 
diameters  (a)  and  the  exact  form  of  the  interparticle  potential  functions  (#). 

It  is  clear  that  more  accurate  calculations  of  the  viscosity  of  ionized  cesium 
require  experimental  data  on  e ,  a  and  <p  which  is  at  present  lacking. 

In  principle  the  above  treatment  may  be  applied  to  any  partially  ionized 
monatomic  gas  provided  few  excited  and  molecular  particles  are  present  and  the 
appropriate  potential  functions  and  numerical  quantities  are  employed. 

21 .3  THE  VISCOSITY  OF  HELIUM -CESIUM  MIXTURES 

Detailed  evaluation  of  the  performance  of  an  MPD  generator  using  helium- 
cesium  mixture  as  the  working  fluid  requires,  amongst  other  factors,  basic  physical 
data  relating  to  the  gas  mixture.  Few  experimental  results  are  available  for  such 
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mixtures,  or  even  for  the  pure  gases  over  the  temperature  range  required  (about, 
1000v  to  3000°K )?  thus  theoretical  estimates  must  be  obtained.  In  this  section 
the  helium-cesium  viscosity  is  evaluated  as  a  function  of  temperature  and  cesium 
concentration.  The  viscosity  is  important  for  determining  the  flow  parameters 
in  an  MPD  generator  and  thus,  for  example,  the  power  losses  sustained  through 
reverse  currents  in  the  boundary  layers- 

21.3.1  Evaluation  of  viscosity 


The  previous  section  has  considered  the  effects  o*'  ionization  on  the  viscosity 
of  pure  cesium  and  shown  that,  provided  the  fractional  ionization  is  less  than  l£, 
the  presence  of  ions  and  electrons  does  not  significantly  alter  the  viscosity. 

Since  the  factional  ionization  of  a  helium-cesium  mixture  in  thermal  equilibrium 
at  the  temperatures  and  pressures  of  most  MFD  generators  is  less  than  1%,  the 
effects  of  ionization  will  be  neglected  in  the  present  evaluation  of  viscosity. 

(Note  that  the  neglect  of  ionization  may  not  bo  valid  if  non-thermal  ionization 
occurs  in  the  gas  mixture;  in  this  case  the  method  of  treating  the  problem  will 
depend  on  the  ionization  level.)  Thus  the  problem  is  reduced  to  evaluating  the 
viscosity  of  a  binary  mixture  of  helium  and  cesium.  The  first  approximation  to 


the  viscosity  is  then  given  by  equation  21.18,  where 

x  _  (1“*2)2  +  2(1-x2)x2  +  *jL 

^1  **12 


(21.26) 
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In  the  above  expressions,  the  subscripts  1  and  2  refer  to  helium  and  cesium 
respectively,  and  x ^  is  the  mole  fraction  of  cesi  a,  defined  by 


*2  = 

requires  values  of  the  collision  integrals  Q 


•  *  *  * »  (21 .33 , 
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which  are  obtained  from  equations  (21.14)  and  (21.15). 


Evaluation  of  77  . 

O  (1,1)*  . 

ni2  ni2 

Information  on  helium-helium  interactions  shows  that  the  Lennard-Jones  6-12  potential 

reproduces  experimental  results  very  accurately  for  the  present  data  above  10°C^ 

(p560);  this  potential  has  the  form  given  in  equation  21.21.  No  experimental  data 

are  available  for  the  interaction  potential  for  cesium-cesium  encounters,  while 
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for  cesium-helium  collisions  the  interaction  potential  has  the  asymptotic  form 
of  the  van  der  Walls  force: 


<p(r)a  .....  (21.34) 

r 

In  this  evaluation  the  Lennard-Jones  6-12  potential  has  been  assumed  for  these 
interactions  also.  The  collision  integrals  for  this  potential  have  been 
evaluated  and  are  given  in  Hirschf elder  et  al1^  (Table  1M)  as  functions  of  the 
reduced  temperature  T*  =  kT/V .  Thus  evaluation  of  T^*,  T^*  ^2*  *'rom  bbe 

critical  interaction  energies  and  e2  enab^es  tbe  collision  integrals  to  be 

determined  and,  knowing  the  atomic  diameters  and  masses,  the  viscosity  of  the 
binary  mixture  may  then  be  obtained. 


Information  on  c .  and  0.  for  helium-helium  interactions  is  rather  varied  at 
present.  According  to  Chapman  and  Cowling  (p.395) '• 

c J .  -  6. 03^  n.  =  2.70  1  from  viscosity  measurements 

1  k  1 

and  c4 /,  =  6.03°K,  a.  =  2.63  from  virial  coefficient  measurements 
Ik  1  16 

while  from  Hirschf elder  et  al  (p.1110): 

=  10.22°K,  <7^  =  2.576$  from  viscosity  measurements 


«nd  «,/k  =  10.22°K,o,  =  2.556  *  (qwmt™  oalqulation)  ^from  ylrlal  coefficient 

J  measurements 

a  6.03°K,o^  =  2.63  $  (classical  calculation) 

Initial  calculations  carried  out  here  with  =  iO°K  and  c ^  =  2,6$  gave 

theoretical  values  of  77  consistently  less  than  those  observed  experimentally. 

A  better  fit  with  experimental  values  is  obtained  using  =  10°K  and  cr-2,5  $. 
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Although  no  information  is  available  for  cesium,  examination  of  data  for  other 

monatomic  elements  shows  that  c-A  =  200' K  and  =  4&  will  not  be  unreasonable 

2  k  <i 

values,  (The  value  of  <r?  used  here  will  be  more  a  curate  than  the  3»5A  used 

22  ^  22 
previously*.  The  conclusions  reached  by  Robinson  and  McNsb  D  namely  that  the 

effects  of  ionization  are  small,  are  not  affected  by  this  change.)  Following 

the  usual  procedure  of  defining  the  interaction  energy  and  effective  diameter  for 

helium-cesium  encounters  by  geometric  and  arithmetic  means  of  the  values  for  the 


pure  gases,  that  is 


*12  =  (*1*2)T 


-  (21.35) 


0^2  =  (cr^  4-C?2  )/ 2 


(21.36) 


yields  (  ^  k  “  K  and  r ^  =  3.25  A 

For  this  interaction  it  is  possible  to  determine  the  van  der  Waals  constant 
C  in  equation  (21.35)  using  the  above  values  of  s,,p  an<^  ^u3  ^rom  equation 

(21.21)  and  (21.35) 

C12  *  **12  °126  *  29,2  x  10'60  "*■ 

This  value  of  C  may  be  compared  with  these  found  experimentally.  Rosin  and  Rabi 

and  Estermann,  Foner  and  Steri?  using  molecular  techniques  found  p  =  29.4-  and 

370  x  10“^°  cm^  respectively,  while  from  measurements  of  the  polarizabilities  of 
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alkali  metals  by  Scheffers  and  Stark  ,  and  polarizabilities  of  rare  gas  atoms 

obtained  from  refractivity  measurements  a  value  of  44  is  obtained.  Thus  the 

value  of  used  here  is  in  excellent  agreement  with  Rosin  and  Rabi's  results 

and  in  fair  agreement  with  polarizability  measurements  but  is  very  different  from 

oi 

the  value  obtained  by  Estermann  et  al.  Further  measurements  on  this  interaction 
are  required  to  clarify  the  situation. 


For  the  present  calculation  the  values  of  f .  and  a  „  which  gives  C,p  - 
-60  6  £.  c  c 

29.2  x  10  cm  have  been  used,  together  with  the  values  of  <4,  and  Op 

previously  quoted.  With  these  values  and  equations  (21.18)  (21.26-33)  the 

viscosity  of  tho  mixture  has  been  obtained  as  a  function  of  temperature  and  mole 

fraction  of  cesium  Xp. 

21.3.2  Results  and  discussion 

The r esults  of  the  calculations  are  presented  graphically  in  Fig.  21.6,  7 
and  8.  Fig.  21.8  illustrates  the  variation  of  viscosity  with  temperature  for 
pure  helium  (xp  =  0),  pure  cesium  (xp  =  1.0)  and  binary  mixtures  with  mole  fractions 
of  cesium  (x>)  of  0.01,  0.1  and  0.5-  The  viscosity  of  pure  cesium  may  be  expressed 

aa  =  6.8  x  10~k  gs/cm  se 
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V 


0.6^  t 

-  5.1  x  10  T  -  gn/cm  see 

over  the  range  200°  to  3000°K,  T  being  in  “K.  A  comparison  between  experimental 
^123 

results  *  and  the  present  theoretical  values  (Fig.  21 .9) ,  shows  that  the 
agreement  for  helium  is  within  about  6$  over  the  range  of  measured  values.  Better 
choice  of  the  form  of  the  interparticle  potential  function  may  remove  this 
discrepancy. 

From  Fig.  21.6  the  addition  of  cesium  to  x ^  -  0.01  increases  the  viscosity 
by  about  while  for  =  0.10  the  viscosity  is  about  4QJfc  greater  than  that  of 
pure  helium  above  2000°K.  As  the  mole  fraction  of  cesium  is  increased  still 
further  the  viscosity  of  the  mixture  becomes  "eater  than  that  of  *pure  cesium,  as 
exemplified  by  the  curve  x ^  -  0.5  of  Fig.  21  .  according  to  Chapman  and  Cowling 

(p.23l)  the  explanation  of  this  effect  is  as  follows.  The  addition  of  a  quantity 
of  a  heavy  gas  to  itself  causes  two  opposing  effects  (on  the  viscosity)  -  a 
reduction  in  the  mean  free  path  and  an  increase  in  the  number  of  carriers  of 
momentum  -  which  just  balance.  If,  however,  a  lighter  gas  is  added  to  the  heavy 
gas,  the  small  additional  transport  of  momentum  by  the  lighter  rarticles  more  than 
outweighs  the  decrease  in  the  mean  free  path  of  the  heavy  particles  caused  by  the 
addition*  (The  alteration  in  mean  free  path  of  the  heavy  particles  by  much 
er  particles  is  very  small.)  Thus  the  viscosity  of  the  mature  may  be 
higher  than  that  of  either  of  the  pure  gases.  This  effect  may  be  seen  more 
clearly  in  Fig.  21.10.  According  to  Chapman  and  Cowling'  (p.233)  the  maximum  in 
the  binary  viscosity  should  become  less  pronounced  as  the  temperature  increases, 
the  viscosity  of  a  hydrogen-hydrochloric  acid  mixture  being  given  as  an  example. 
However,  this  does  not  occur  for  a  cesium-helium  mixture  (Fig.  21 .10);  the  maximum 
in  fact  becomes  more  pronounced. 

21 .4  VISCOSITY  MAXIMA  IN  BINARY  GAS  MIXTURES 

21 .4.1  Introduction 

The  calculations  reported  in  section  21.3,  using  the  first  approximation  to 
the  viscosity  of  a  gas  mixture,  showed  that  for  helium-cesium  mixtures  the 
viscosity  of  the  mixture  was  higher  than  either  of  the  individual  components; 
that  is,  the  addition  of  a  light  gas  to  a  heavy  gas  increased  the  viscosity  of 
the  latter  (in  general,  but  not  always,  the  viscosity  of  a  heavy  gas  is  greater 
than  that  of  a  light  gas).  The  explanation  '>f  this  effect,  put  forward  by 
Chapman  and  Cowling11 ,  was  outlined  in  section  21.3.2. 
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According  to  Chapman  and  Cording'  the  maximum  in  the  binary  viscosity 
should  become  less  pronounced  as  temperature  increases,  the  viscosity  of  a 
hydrogen-hydrochloric  acid  mixture  being  given  as  an  example.  However  the 
calculations  reported  in  21.3  showed  that  this  was  not  true  for  helium-cesium 
mixtures  -  the  maximum  becomes  more  pronounced  as  temperature  increases.  In 
view  of  the  doubt  associated  with  the  data  used  for  cesium  in  these  calculations, 
further-  work  has  been  carried  out  on  this  effect  for  other  gases  for  which  more 
reliable  data  is  available. 


21 .4.2  Theory 

The  effect  of  the  introduction  of  a  light  gas  into  a  heavy  -”'s  is  most 
conveniently  examined  by  the  use  cf  the  1 impurity'  approximation  for  the  viscosity 
of  a  binary  mixture  given  in  Ref.  1  I .  This  expression  is  derived  for  the 
condition  that  only  a  small  amount  of  impurity  gas  (the  light  gas  in  this  case) 
is  added  to  the  main  gas. 


Writing  Chapman  and  Cowling' a  impurity  approximation  in  the  nomenclature 

1  g 

employed  by  Hirschfelder,  Curtiss  and  Bird’  : 

r3i*  (m1+m2)2  .  x*  .  2ti  2 

5^12  2mTMq  ^'2  + 


1  -x_ 


''mix  =  n2  *  *. 


,  +  iA' 

•  *  c  '  — — 


*'■¥.2  ^  -  Sri 


(21 .37) 


4  R 


2 


The  lighter  gas  is  species  1,  the  heavier  gas  species  2. 

It  is  apparent  that  the  viscosity  of  gas  2  wi^i  be  increased  by  the  addition 
of  the  lighter  gas  only  if 

A  >  0  .....(21.38) 

where  n 


A  = 


_  1 


(W 


5  “12  2  U^2 


'1 2 


+  2(l A 


12*^ 


2n, 


2  f?l2 


....  .(21  .39) 


21 .4.3  Calculations  and  results 

A  Pegasus  autocode  computer  programme  has  been  used  to  investigate  equations 
(21.38)  and  (21.39)  for  increasing  temperatures  for  various  gas  mixtures.  The 
viscosities  an^  ^ave  been  evaluated  assuming  that  the  Lennard-Jones  (6-12) 
potential  is  operative,  thus  an  unequal  internal  interpolation  routine,  together 
with  the  data  of  Table  21.2  here  and  Table  1R  cf  Ref,  1b,  was  u3ed  to  evaluate  the 
appropriate  collision  integrals.  (Tho  data  for  the  Lennard-Jcnes  force  constants 
in  Table  21.2  is,  with  the  xception  of  cesium,  taken  from  Ref,  24). 

The  results  (values  of  A  at  different  temperatures)  fall  into  three  main 
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ategories: 

mixtures  which  show  no  viscosity  maximum; 

mixtures  which  have  a  viscosity  maximum  at  low  temperatures 
which  disappears  as  the  temperature  increases, 
rd  a  third  group  of  mixtures,  unpredicted  in  Ref  11,  having  a  viscosity  maximum 
hich,  in  general,  becomes  more  pronoun  ed  as  the  temperature  increases.  These 
ixtures  are  3howr.  in  Table  21.3*  In  all  these  cases  the  first  component  is  the 
ighter  gas  (species  l)  and  >  1  •  The  mixtures  in  the  first  column  are 

uoted  in  Ref,  24  as  having  no  maximum  and  those  in  the  other  two  columns,  with 
he  exception  of  A-Cs  and  He-Cs  as  having  a  maximum. 

The  quantity  A  has  been  evaluated  over  a  wide  range  of  temperature  (200°  to 
000°K)  and  typical  results  for  several  gas  mixtures  in  each  group  are  shown  in 
ig.  21.11.  No  account  has  been  taken  of  dissociation,  ionization,  vibrational 
r  rotational  energy  modes  or  chemical  interactions  in  these  calculations  although 
hese  effects  will  undoubtedly  become  important  at  high  temperatures  for  some  of 
ae  gases.  It  is  probable  from  this  point  of  view  that  the  results  for  the  inert 
as  mixture^  will  be  the  most  re.' .  able. 


21.4.4  Discussion  and  conclusions 

It  has  been  found  that  for  a  number  of  h-'  ary  gas  mixtures  the  viscosity 

lows  a  maximum  which  becomes  more  pronouncer.  ^s  temperature  increases.  This  is 

1  1 

intrary  to  the  example  quoted  by  Chapman  and  Cowling  (H^-HCl).  The  criterion 
iven  in  equation  (21 .38 3  and  (21.39)  for  a  viscosity  maximum  is  written  by 
lapman  and  Cowling  in  a  rather  different  form: 


4A  J?2 


(21 .40) 


E  =  |  nmo  (D12)1  and  A  =  |  A12* 


. (21.41) 


“3*o'  1  .  512  .  ' 

nee  A  is  almost  temperature  independant  the  dominant  temperature  dependant  term 

r  constant  pressures: 


g  =  Z  £E£  (p  ) 
}  iff  '  12'l 


3 /  2  1/  2 

!>  C  It  ^  f 


a,  for  hard  sphere  collisions,  (D. 0).  ex  T  ,  thus  S  a  T  .  Typically 

1  <-  1  Q  gC  Q  ^  £ 

rrati-ns  of  r with  temperature  are  approximately  as  T  *'  ,  thu3  S/ri^  oc  l/T  ’  . 

rce  if  the  inequality  of  equation  (21.40)  is  true  at  a  given  temperature,  as  the 
mperature  increases  it  will  eventually  cease  to  be  true,  that  is,  the  viscosity 
ximum  w*ll  vanish,  as  stated  in  Ref.  11. 
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However,  the  diffusion  coefficient  varies  »a  i 
collisions , 


3/- 


jnly  f o’ 


sphere 


For  unlike  molecules  repelling  each  utner  witn  a  torce  varying 


inversely  as  the  v ,  „  ^  power  of  the 


12 

Di  2  K  T 


di  3  tan 


.  a  is  t  m  1 


1+S 


where 


s  =  -  +  — ~* 
2  w.-.-l 


12 


iA 


nowever  v 


12 


For  hard  spheres,  v  ^  =  oc,  S  =  '2  and  _  «  T"‘  2,  as  given  above, 
is  usually  in  the  range  3  to  11  (Ref  11,  p  243.  25 /)■  Pgr  ^2  *  ®  ~  ^  ar,a 

D  oc  T1  thus  for  77^  oc  T0*^  oc  T0*^  -  Thus  is  the  reverse  case  of  that 

quoted  ,.oove  since,  as  the  temperature  increases  the  viscosity  maximum  becomes 
more  pronounced. 

2  £7 

A  physical  explanation  has  recently  been  put  forward  by  Cowling  .  To 
take  further  the  comments  in  ref.  11  (p  269)  ?  there  is  a  maximum  in  the  viscosity 
only  if  the  lighter  gas  has  a  smaller  equivalent  molecular  radius  than  the 
heavier  an(i  the  maximum  disappears  with  increasing  temperature  only  if 

the  difference  between  the  two  molecular  radii  becomes  smaller  as  the  temperature 
increases.  For  gases  with  no  maximum  o,/ t  at  0  C  (N  b.  Tnese  c  are  not  tnose 
of  Table  21.2)  is  greater  than  about  V3,  for  gases  with  a  maximum  always  cycg 
is  less  than  about  0.6;  for  gases  with  between  these  values  the  maximum 

may  disappear  with  increasing  temperature.  It  ha=  therefore  been  shown  that 

some  gases  exhibit  a  viscosity  maximum  whicJr  becomes  more  pronounced  as  temperature 
increases,  contrary  to  the  more  usual  cases  where  tne  maximum  disappears.  At 
the  present  time  no  experimental  data  is  known  wt ich  supports  this  conclusion, 
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APPENDIX  21 -A 


CALCULATION  OF  DEFLECTION  ANCLE  FOR  ATOM -ION  COLLISION 


When  the  form  of  the  potential  4> 


bstituted  from  (21,22),  equations 


(21,15)  and  (21.23)  give 

8  =  b*g* 


ai 


(r  *  ) 
min 


-1 


_ du^ _ 

7^2  ,  *2  m2  .2  *i-l72 

(g*  -  b*  g*  u*  r  u*  ,1 


(21 .A.1 ) 


If  the  force  law  were  that  of  a  simple  attraction  with  no  impenetrable  core, 
the  reduced  distance  of  closest  approach  of  the  centres  would  be  the  reciprocal  of 
the  smallest  root  u0*  of  the  equation 

0  .....  (21. A. 2) 

2  2 

which  has  real  roots  if  and  only  if  b*  *  When  the  roots  are  real,  the 

reduced  distance  of  closest  approach  will  be  u$*  or  1,  whichever  is  greater, 


0  ?  ?  2 

-  b*  g*‘  u*  +  u* 


since  under  no  circumstances  can  the  pair  approach  more  close _lv  than  a  , 


cases  in  all  may  be  distinguished: 


ai 


-1 


Three 


(l)  if  b*2  >  1  +  "”~r  ,  equation  (21  .A. 2)  has  real  roots  and  uD*  >1.  Thus 

g 


r  ,  *  ~  u0* 
min 


-1 


(2)  if  Jr  <  t*2  <  1  .  -n  , 

8  g* 

therefore  r  ,  *  =  1 . 
mm 


-1 


equation  (21 .A.?)  has  real  roots  and  u0*  <  1: 


2  2  + 

(3)  if  b*  <  —  ,  equation  (21, A. 2)  has  no  real  roots,  so  =  1 


In  cases  (l )  and  (2)  the  integral  may  be  simplified  by  the  substitution  u*  = 

20 


u0*  sin  $  employed  by  Thomson  .  In  case  (1 )  the  result  is 
6  =  b*u*  I  - — T - T77 


while  case  (2)  yields 

8  =  b*Uo' 


where  §  -  arc  sin  (u|  ^). 


(1-aC  Bin2,)" 

g*2 


M. 


(21  .A. 3) 


*4  ,  V2 

(1  -  ^  sir/  *,) 


(21 .A .4) 


.  (21 .A. 5) 


In  case  (3)  the  internal  is  reduced  to  standard  form  by  the  substitution 

2  u*2 

tan  B  i— 

2  g * 


0  1  1 

U  !  •  i  y 


♦  *  *  *  * 


(21 .A. 6) 


which  gives  $  =  -  b*^ 

where  r  =  2  arc  tan  (^7  )  . (21  *A,7) 

Thus  in  all  cases  6,  and  hence  y  (eqn.  (21.  1$)  can  be  expressed  in  terms  of 

an  elliptic  integral  of  the  first  kind.  Equations  (21 .A. 3),  (21  .A.4)  and  (21 .A. 6) 

(1  4 )+ 

were  used  to  obtain  values  of  y  from  which  the  collision  integrals  Q  ;  *  7  and 
(2  2)*  ai 

°.i  ’  given  by  equation  (21.14)  could  be  evaluated  by  Gaussian  quadrature. 
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TABLE  21 ,1 


COMPUTED  COLLISION  INTEGRALS  (Q 


ai 


AM 


ai 


FOR  ATOM-ION  INTERACTION 


TABLE  2*  .2 


DATA  USED 

IN  VISCOSITY  CALCULATIONS 

GAS 

<j»* 

U 

H2 

36.7 

2.958 

2.016 

Ha 

6.03 

2.63 

4.004 

C\ 

137 

3.882 

16.047 

kh3 

320 

2.60 

17.036 

w* 

33.74 

2.756 

20,177 

CO 

88.0 

3.706 

28,019 

*2 

79.8 

3.749 

28.022 

C2\ 

205 

4.232 

28.063 

C2H6 

230 

4.418 

30.080 

°2 

88,0 

3.541 

32.008 

HC1 

360 

3.305 

36.480 

A 

124 

3.418 

39.960 

C°2 

213 

3.897 

44,024 

C3H8 

254 

5.061 

44.111 

S02 

252 

4.290 

64.671 

Kr 

166.67 

3.682 

83,830 

It 

225.3 

4.070 

131.341 

Cs 

250 

4.1 

132.95 

* 

Units  of  °K 

*• 

Units  of  Angstroms 

+ 


Estimated 


TABLE  21 .3 


RESULTS 


Wo,  maximum  Maximum  which  disappears  Maximum  always 


ch4-co2 

H2-HC1 

H2'C2H4 

C2H4*A 

H2-Ch4 

H2-C2H6 

C2V°2 

H2-NH3 

H2-S02 

co-o2 

H2^°2 

V°2 

H2-Xe 

He -A 

H  -0 

A  -Cs 

rie-Kr 

2  2 

H2-Ne 

He-Xe 

He-Cs 

FIG. 2  POTENTIAL  FOR 
ATOM-  ION  INTERACTION 


FIG.  3  POTENTIAL  FOR  ION-ION  INTERACTION. 

FIG  21-1.2  &3 


COLLISION  INTEGRAL 


REDUCED  TEMPERATURE  T* 


COLLISION  INTEGRAL  iL^*  AS  A  FUNCTION  OF  REDUCEO 
TEMPERATURE  T£  FOR  LENNARO  -  JONES  POTENTIAL 


FIG  214 
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REDUCED  TEMPERATURE  T  * 


COLLISION  INTEGRALS  JL1^*  AND  Jl  (Il;*AS  A  FUNCTION 
OF  REDUCED  TEMPERATURE  T^*  FOR  ATOM-ION  INTERACTION 


REDUCED  TEMPERATURE.  T 


ION- ION  COLLISION  INTEGRAL  AS  A  FUNCTION  OF  THE  REDUCED 


TEMPERATURE 


viacosrrv  (»o  q/cm 


COSITY  OF  IONIZED  CESIUM  AS  A  FUNCTION  OF  TEMPERATURE 


VISCOSITY  AS  A 


VISCOSITY,  am/c^scc 


VISCOSITY  AS  A  FUNCTION  OF  MOLE  FRACTION  OF  CESIUM 

FIG  2110 
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PARAMETER  A  AS  A  FUNCTION  OF  TEMPERATURE 


FIG  21  II 
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Defined  by  equation  (21  ,4  l) 

coefficient  defined  in  equation  (21 .20)  or  (21 .32} 

impact  parameter 

van  der  Waals  constant 

diffusion  coefficient 

defined  by  equation  (21  .41 ) 

electron  charge 

velocity  distribution  function 

relative  velocity 

Debye  length 

ionization  potential 

Boltzmann  constant- 

molecular  weight 

particle  mass 

number  density 

pressure 

separation  of  two  interacting  particles 

temperature 

time 

deflection  time 

-1 

r 

smallest  root  of  equation  (21  Jt.2) 
particle  velocity 

defined  in  equations  (21.19)  or(21 .26)  to  (21.28) 

mole  fraction 

position  co-ordinates 

fractional  ionization 

arc  sin  (uo*  ) 

?  arc  tan  (  /vg*) 
defined  by  equation  (21 .39) 

critical  value  of  energy  of  interaction  of  two  particles 
viscosity 

viscosity  of  two-component  mixture 
is  defined  in  equation  (21 .23) 

ratio  of  critical  impact  parameters  (equation  21.,’) 


mean  free  path 


V 


force  power  constant 

p  min.  separation  of  centres  of  two  atoms 

p *  lain,  separation  of  centres  of  atom  and  ion 

o  critical  value  of  separation  of  two  particles 

<P  potential  energy  of  two  interacting  particles 

X  angle  of  deflection  of  each  of  two  colliding  particles  in  a 

centre-of-aass  frame  of  reference 

tp  »  arc  sin  u*/uc* 

q(p»3)  co3j_islon  integral  defined  in  equation  (21 .14) 

Subscripts 
a  atoms 

e  electrons 

i  ions 

direction  co-ordinates 
min  minimum  value 

1,2  components  of  a  binary  mixture,  particularly  helium  and  cesium,  respectively 

-  vector 

Superscripts 

*  reduced  (dimensionless)  quantity  (see  equation  (21.16)) 


mean  value 


22.1  INTRODUCTION 


COLLISION  CRO 


I.R.  MeNab 


Many  experiments  presently  being  performed  on  methods  of  generating 
power  directly  with  thermionic  diodes  or  magnetoplasmadynamic  devices  use 
cesium,  pure  or  diluted,  as  the  working  fluid.  To  evaluate  the  electrical 
conductivity  of  the  gas  or  plasma  the  contribution  made  by  collisions  of 
electrons  with  neutral  cesium  atoms  must  be  known.  An  increase  in  the 
number  of  such  collisions  per  unit  time  reduces  the  mobility  of  the  electrons 
and  hence  the  electrical  conductivity.  The  relationship  is 


(22.1) 


where 


M  mu 


,(22.2) 


c  is  the  scalar  conductivity,  n  ,  s  and  m_  are  the  number  concentration, 
charge  and  mass  of  the  electrons,  y  is  the  collision  frequency  for  momentum 
transfer  between  electrons  and  cesium  atoms  and  u  is  the  electronic  mobility. 

Equation  (22.1)  is  valid  only  where  electrons  and  cesium  atoms  are  present; 
if  other  particles  exist  in  appreciable  numbers  the  collision  frequency  is 
modified. 

For  many  gases  y  may  be  replaced  by  v  ,  the  elastic  scattering  collision 

m  ”0 

frequency,  since  the  values  agree  to  within  10  per  cent;  however,  recent 
theoretical  calculations  (discussed  below)  indicate  that  this  is  not  valid  for 
cesium. 

The  elastic  scattering  collision  frequency  (y  )  between  electrons  and 
neutral  atoms  is  given  by 


=  N  Q  v 
o  e 


(22.3) 


where  N  is  the  number  concentration  of  neutral  atoms,  v  Is  the 


n  o  * 


0=0  $=0 

where  dc/dQ  is  the  differential  cross  section  for  scattering  the  electron  into 
the  solid  angle  dO  =  sin  0  d0  d$. 

The  momentum  transfer  collision  frequency  is  given  by  a  similar  expression; 

V|(  =  N  Qjj  v  ••-•(22.5) 


where  Q..  is  the  momentum  transfer  cross  section; 


,  2w  .  2w 

/  1 

0  =  0  <£=0 


d? 


^  (1  -  cos  sin  0  d0  d# 


and  the  mean  fractional  energy  loss  per  collision  is 


S  = 


2  me 


*  Q, 


(22.6) 


(22.7) 


M  being  the  atomic  mass. 

The  cross  sections  and  collision  frequencies  averaged  over  a  Maxwellian 
distribution  may  be  found  from  the  following  equations; 

average  elastic  scattering  cross  section 

\  -  I  f0  (v)  %  (v)  dv  . .... (22,8) 

Jo 


average  momentum  transfer  cross  section 
-  I  fQ  (v)  Qy  (v)  dv 


(22.9) 


average  momentum  transfer  collision  frequency  per  atom 


v 

JT  =  f0  (v)  Qy  (v)  V  dv 

Jn 


(22.10) 


and  the  average  total  collision  frequency  is  defined  in  a  similar  fashion. 


22.2 


4  7?V  1 


.(22.11 


is  the  Maxwellian  distribution  of  velocities  for  the  electrons;  it  as  the 
distribution  most  often  assumed  in  plasmas,  although  not  always  valid. 
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LXFE  RI  MENTAL  I K VEST I CATI OK S 


Until  recently  experimental  data  on  the  cross  sections  of  cesium  (and  other 

t 

elements)  at  electron  energies  below  one  electron  volt  have  been  sparse.  Erodes 
measurements^  cover  a  large  energy  range,  from  electron  energies  greater  than 
lOOeV  to  about  O.S  eV;  the  low  energy  values  obtained  are  shown  in  Fig.  22.1. 

The  results  obtained  by  Boeckner  and  Kohler  have  recently  been  corrected  by 
Phelps  in  the  light  of  later  experiments  by  Mohler',  Phelps  found  that  the 
collision  frequency  per  atom  was  constant  over  the  range  'Mj.2  to  0.4  eV 
(Fig,  22.2)  ;  a  collision  cross  section  derived  from  these  results  using  the 
relation 


v  =  jNQ  v 


.(22.12) 


is  shown  in  Fig.  22.1.  Recent  experiments  at  low  electron  energies  in  cesium 
appear  to  have  confused,  rather  than  clarified  the  situation. 

Using  a  conductivity  cell  with  concentric  and  parallel-plate  electrode 
geometry  ,  Mullaney  has  derived  collision  cross  sections  which  are  lower  than 
those  of  Erode  or  Mohler  (Fig.  22.1).  The  cross  sections  are  averaged  over  some 
electron  distribution;  as  with  most  experiments  of  this  type  the  electron  energy 
distribution  is  not  known  accurately. 

Ro ehling"  has  used  a  similar  apparatus  with  a  parallel-plate  electrode 
geometry  to  measure  cross  sections  up  to  about  2000  K.  The  results  agree  closely 
with  Mullaney1  s,  being  lower  than  those  of*  Erode  and  Mohler;  a  Maxwellian  velocity 
distribution  for  the  electrons  is  claimed. 


Mongolia  and  Korchevoi b  ,  using  an  arc  discharge  between  tungsten  cathodes 

^  o  2 

in  cesium  vapour  measured  an  electron-cesium  atom  cross  section  of  --50x1u  cm* 

This  value  is  similar  to  those  found 


at  an  electron  temperature  of  about  5QQQ°K. 


by  Mullaney  and  Roehling,  but  at  a  much  higher  temperature  (Fig.  22.1  . 

Chen  and  Raether?  conducted  a  series  of  experiments  on  helium-cesium 


mixtures  and  pure  cesium  plasmas  using  microwave  technique 


In  the  temperature 


range  450°  to  550°K  they  find  that  the  average  momentum  transfer  cross  section  is 


22.3 


best  represented  by 


<4 


1.61x10 

T 


-20 


9.63x10 


-12 


+  2.03x10  ^  cm' 


e 


These  collision  cross  sections  are  shown  in  Pig.  22.1  where  it  can  be  seen  that 
they  are  much  higher  than  other  values  previously  found, 

g 

From  electrical  conductivity  measurements  in  various  gas  mixtures  Harris 
finds  an  average  momentum  transfer  collision  cross  section  for  cesium  of  about 
300  x  10  ^  cm^  and  a  collision  frequency  of  8  x  10  ^  cmV3ec  at  a 

temperature  of  1750°K. 

Using  the  method  of  phase  shifts  and  including  the  effects  of  polarisation, 

9 

Robinson  has  evaluated  theoretically  the  average  elastic  scattering  cross  section, 
momentum  transfer  cross  section  and  the  momentum  transfer  collision  frequency  for 
cesium;  the  results  are  shown  in  Pigs. 22.1  and  22*2.  In  view  of  the  wide 
scatter  of  experimental  data  Robinson’s  calculations  will  be  useful  in  estimating 
the  electrical  properties  of  cesium  plasmas  in  the  energy  range  belcw  1  eV.  As 
Robinson  points  out,  the  calculations  are  of  a  preliminary  nature  and  more 
accurate  theory  may  be  required  to  obtain  complete  agreement  with  experiment. 

Measurements  in  a  static  cesium  plasma  by  Ralph1 ^  gave  values  for  the 

—1 6  2  —1 5  2 

cesium  cross  section  of  2390  x  10  cm  and  1370  x  10  cna  at  a  temperature 
of  1050°K.  These  are  much  higher  than  other  results  reported  by  Ralph11  at 
1500°K,  which  yielded  200  x  10  ^  cm^. 

22.3  CONCLUSIONS 

Although  several  measurements  of  collision  cross  sections  and  collision 

frequencies  have  recently  been  made  for  cesium,  experimental  data  is  still  widely 

divergent.  If  all  the  experimental  data  is  correct  the  cross  section  of  cesium 

atoms  at  low  energies  exhibits  a  hitherto  unexpected  energy-dependence.  The 

9 

theoretical  calculations  of  Robinson^  li^.  within  the  range  of  the  experimental 

values.  For  most  of  the  calculations  in  other  Chapters  of  this  report  a 

-1 6  2 

constant  ce3ium  atom  cross  section  of  300  x  10”  cm  has  been  used;  reference 
to  the  experimental  values  shows  that  this  is  a  reasonable  approximation, 
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SEED  FRACTION  OPTIMIZATION 


by 


I.R.  Me-Nah 


23.1  INTRODUCTION 


In  analyses  of  MPD  generators  in  which  the  electrical  conductivity  of  the 

parent  gas  is  increased  by  the  addition  of  a  seed  element  the  scalar  conductivity 

is  usually  quoted  as  ? 

n 


where , 


<7 

o 


if  electron- 


e 


n  interactions  can  be  neglected 


•  •  *  t  • 


(23.1) 


=  V 


e**r 


+  v 


e-a 


'  ve  ( “P^P  *  Va  na)  . (25,2) 

If  the  seed  fraction  ( y )  is  small  (y  =  n E/np)  maximizing  conductivity  with 
respect  to  y  shows  that  the  optimum  seed  fraction  is 

X  .  =  q  r/o  .....  (23.3) 

'  opt  Me-P  'e-a 

is  usually  about  0.01  for  alkali  metal  seeded  inert  or  combustion  gases. 

In  view  of  the  rather  sparse,  and  often  suspect,  data  on  collision  cross  sections 
at  the  low  electron  energies  relevant  to  MPD  generation,  y  ^  is  often  assumed 

Op  v 

cons  tant. 


Considerable  attention  is  presently  being  focussed  on  the  production  of 
extrathermal  ionization  in  closed  cycle  MFD  generators  by  various  means,  but 
particularly  on  magnetically-induced  ionization  by  which  interaction  of  the 
electrons  with  the  electric  field  m  the  generator  causes  the  electron  temperature 
to  be  considerably  increased  over  the  gas  temperature.  When  this  occurs  ionic 
scattering  of  electrons  becomes  an  important  factor  in  the  determination  ot  the 
electronic  collision  frequency  and  hence  electrical  conductivity.  Since  the 
electron-ion  collision  frequency  is  (quite  strongly)  dependent  on  the  electron 
temperature  the  optimum  seed  fraction  is  no  longer  the  unique  quantity  defined 
by  equation  (23.3)  but  is  a  function  of  the  working  conditions  of  the  generator 


23.1 


even  when  q  and  q  are  constant.  The  selection  of  the  optimum  seed  fraction 
Sr  ©a 

for  this  case  and  the  consequences  of  the  temperature  dependence  are  examined  here. 
23.2  SIMPLE  THEORY 

It  will  be  assumed  that  the  Saha  equation  can  be  used  to  obtain  the  fractional 
ionization  of  the  gas  at  the  electron  temperature  so  that 

x2 

— P  pw.  =  AT  )  . (23.h) 

1-X  As  6 

5/0  2wm  3/ 2  g  g.  eV. 

where  A?*)  =  (kT^)  "  ( - )  .  exp  (-  -r^r  )  . (23. 5) 

h  sa  e 

In  the  simple  theory  developed  in  this  section  it  will  be  assumed  that  the 

fractional  ionization  is  small  (that  is  X<<1  or  n  <<  n  ).  Under  these  conditions, 

0  & 

Saha's  equation  becomes 

x2  ft,  »  /(T  )  . (23.6) 

Making  the  additional  assumption  that,  while  T  is  greater  than  T,  n  T  «  n  T, 

6  €*  0  & 

then 

Ptj  =  nK  .  (23.7) 

s 

Prp  =  PT  +  Py  =  (na  +  np)  kT  ;  . (23.8) 

s 

n 

and,  by  definition  x  -  — .....  (23.9) 

na  P 

so  that  p  =  xVm  . (23.10) 

AS 

and  Saha's  equs  tion  becomes 

X^  =  - —  .....  (23.11) 

XPT 

The  electron  concentration  may  be  f^und  from 

%  ■  -zr  . ^-12> 

Combining  equations  (23. 5), (23*8), (23. 11 )  and  (23.12)  yields 

lA  V  T  Vi.  72.0  k  \Vv/g  g  V/2  /  eV  \ 

■*  2  7/2  {-f)  Hr)  ^Nsr)--*23-13) 

where  (Ea+np)  can  8e  inserted  as  data  appropriate  to  chosen  gas  temperatures  and 
total  pressures.  The  seeding  atom  concentration  and  parent  atom  concentration 


23.2 


are  obtained  from 


a 


.....  (23.14) 
.....  (23.15) 


3  *t  na  +  V 

np  =  (l-tf)  (nft  +  np) 

For  a  seeded  gas  in  which  ionic  scattering  of  electrons  is  important  the  total 
collision  frequency  is 

v  m  =  v  (q  _rL  +  q  n  +  q.  ,  n  )  .  (23.16) 

eT  e  e-P  P  e-a  a  e-i  i 


where,  from  Spitzer^ 


I  |  4  -*•  \JIU  —  X  /  0 

A  /B'Te 

*--«  =  “2  ^  (  I  1/2' 


e- i 


. (23.17) 
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e 


where 


A  = 


256i'(eok)< 


B»  = 


,  ,3/2 

12 v  (fQk) 


and,  for  a  Maxwellian  distribution, 

1/2 


'  e  \ wm  / 


(23.18) 


Combining  equations  (23.1 ),  (23.13)*  (23.14),  (23*  7)  and  (23.16)  an  expression 
for  conductivity  in  t^ras  of  seeding  fraction  is  obtained 

1/2 


where 


JLx 


°°  F(l-x)  +  Gy  +  Hyl/2  [I  -  j-  In  x] 


(Mil 

'  'ga  1 


1/2 


'2ffm  k 
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3/i 


D  =  eV^/2k 

e2  (VM 

*  =  T  I 
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v  o 


D  \ 
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(23.19) 


F  =  v  q  ^  (n  +rO 
e  e-r  a  r 


G  =  v  q  (n  +n_) 

e  e-a  a  P 


1/2 


D 


v  A  / n  +n  \ 

e  ,  *•  P  j  -  / 

A— y  c  exp 

e  \  e 


I  = 


B»T  T1^ 


.  |  e  |  _D_ 

ln  I  7  xl/4,1/2  J  +  2T 

L(na+np)  0  J  e 


23.3 


Differentiating  conductivity  in  equation  (25.19)  with  respect  to  seeding  fraction 
and  setting  the  result  equal  to  zero  yields  a  quadratic  for  optimum  seeding 
fraction  for  which  the  solution  is 


^opt 


F  +  H2  +  H2 
C’~?  8(0-F)2  8(G-F)2 


I6F1G-F) 


H*" 


l  i 


(23.20) 


For  the  rase  when  ionic  collisions  can  be  neglected  (H  =*  0,  occurring  at  low 
fractional  ionization)  equation  (23.20)  yields 


‘opt 

»  q- 


F 

G~? 


AeP 


. (23.21) 


^eP*  *ea 

so  that  this  reduces  to  equation  (23.3).  (The  difference 


O.ualiy  qea  ..  ,ep 

between  equations  (23*3)  and  (23.21)  is  due  to  the  different  definitions  of 
seeding  fraction  used  in  the  two  cases.) 


When  ionic  collisions  are  important,  equation  (23.20)  shows  that  the 
optimum  seeding  fraction  is  increased  over  the  value  given  in  equation  (23.21): 
the  greater  the  effect  of  ionic  collisions,  the  greater  the  optimum  seeding 
fraction,  up  to  a  maximum  of  unity. 

However,  equation  (23.20)  is  not  very  useful  in  examining  the  effects  of 
extrathermal  ionization  on  seeding  fraction  since  such  effects  generally  cause 
high  fractional  ionizations,  so  that  the  approximations  used  in  deriving  equation 
(23.20)  become  invalid  and  a  more  accurate  approach  is  required. 


23.3  ACCURATE  TREATMENT 

If  the  fractional  ionization  is  high  it  is  necessary  to  use  the  exact  Saha 
equation  (equation  23*4)  in  which  the  total  seed  pressure  is  expressed  in  terms 
of  known  constants  of  the  system  (total  gas  pressure  and  temperature). 


Now 


Pm 

=  n  kT  +  n  kT  +  n  kT 

“x 

ls 

e  e  i  a 

=  nekTe  *  V1' 

.  (23.22) 

since 

n 

e 

”  ni 

and 

n_ 

-  n  +  n 

S 

e  a 

Using 

X 

=  n0/ng  equation  (23.22)  becomes 

PT 

=  tog  (XTe  ♦  T) 

23.4 


cr,  using  a  more  accurate  definition  of  seeding  fraction  than  that  of  section  23*1 , 
that  is  n„ 


X  = 


°T 


where  hj,  =  rip  +  ng  =  +  n&  +  n0, 


(23.23) 


At  constant  pressure  n^,  is  a  variable  so  that  an  expression  is  required  relating 
lPj,  and  p,p.  Using  the  previous  conditions 


=  n  kT  +  n  ,kT  +  n  kT  + 


e  e 


np: 


kT 


=  n  kT  +  n_kT . 
e  e  T 

Now 

ne  ' 

~  Xyrtj, 

so  that  PT  =  (Xx^  +  T) 

Eliminating  from  equations  (23.23)  and  (23.25)  yields 
*(XTe  +  T)pT 
%  =  x^e  +  T 

30  that  the  Saha  equation  (23.4-)  becomes 


(23.24) 

(23.25) 


(23.26) 


X~y(XTe  +  T)pT 

(l-X2)(XyT  +T) 

w 


/(t9) 


(23.27) 


Thus  for  given  T  ,  T  and  p_  this  cubic  equation  must  be  solved  for  X,  When 
X  is  found  ng  may  be  obtained  from  equations  (23.24)  and  (23.25): 


fypj 

ne  "  k(X*Te+T  ) 

.  (23.28) 

and  the  remaining  required  quantities  follow. 

Thus  rVp  =  n/Zx 

.  (23.29) 

np  =  r\p(l-x) 

.  (23.30) 

na  =  V(1"X) 

.....  (23.31) 

The  scalar  conductivity  of  the  gas  is  again  obtained  from  equation  (23. 1 )  and  the 
total  electronic  collision  frequency  from  equation  (23*16).  Writing 
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V.P  = 

F4T 

1  e 

1/2 

v  q  = 

GT 

1/2 

e  e-a 

1  e 

and 

Av  = 

H4T 

1/2 

e 

1  e 

gives 

n  (eV»  ) 

_® _ 8_ 


w-  H  n  /B'T  3//2  \  ..  (23.32) 

T.  ^  i"  ( — w 
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1/2, 


where  A  and  B'  have  been  defined  previously#  Substitution  in  this  equation  of 
the  values  of  ne»  r^,  and  other  parameters  previously  calculated  will  give  the 
conductivity. 

In  Section  23*2  it  was  shown  that  in  the  simple  evaluation  of  electrical 
conductivity  an  explicit  relation  for  the  optimum  seeding  fraction  can  be 
obtained.  In  the  general  case  a  simple  explicit  expression  for  X  ^  cannot  be 
obtained  although  equations  governing  the  value  of  X  ^  may  be  set  up  as  follows, 

Rewriting  equation  (23*32)  in  the  form 

X( e2/ m  ) 


o  = 
o 


w  y  /tw  3/2 

t/V-*)/*  +  T.V%(1-X)  *  -V2  ln(-r/2 

i  \  n 

e  \  e 


and  differentiating  yields 


e2/m 


d*o  = 


P4T 
1  e 


1/2 


Xd*/X2+(1  -x  )dX/*  +G1  Te1^2dX+H1  X2dn 


ivrty 


(23.33) 


Using  a  =  T0/T  6X1,1  b  =  )/ p,p  equation  (23.27)  becomes 

aX3X  (l+b)+X2(y+b)  -  abXX  -  b  =  0 
which,  when  differentiated,  yields 

j^aJC2  ( 1  +b* )  +2X(  1  +b )  -abX  dX  =  £ 

Differentiation  of  equation  (23.28)  gives 

(Xd*  +  XdX)TpT 


a(l+b)X3+X2-abX 


ax 


dn  = 
e 


k(XXTe  +  T)* 


(23.30 


.  (23*35) 


.  (23.36) 


Klinina tion  of  dX,  dn  ,  X  and  n  from  equations  (23.28)  and  (23.33)  to  (23.36) 

0  0 
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will  give  an  equation  for  d o/dX  m  terms  ot  A  ana  var iuua  onjLants, 
condition  do ^/dX  =  o  will  then  give  an  equation  having  as  one  ot  i" 

roots.  Since  it  is  nossible  to  solve  expli  ity  the  general  oubxo 


roots , 


bis  equation, 


it  is  also  possible  to  perform  the  elimination  explicitly  and  obtain  an 

equation  d>(X  )  -  0,  but  the  extent  of  manipulation  required  usually  renders 

H  opt 

this  impracticable,, 

23.4  COMPUTATION  AND  RESULTS 

The  Saha  equation  (23.34)  was  solved  for  fractional  ionization  by  Newton’s 
method  using  (b/l+b)1^  as  a  first  approximation  to  X  This  converged 
satisfactorily  to  the  correct  root  for  the  range  of  values  of  X  and  T$  considered. 
The  particle  number  concentrations  can  then  be  obtained  from  equ  tions  (23*28)  to 
(23.31)  and  the  electrical  conductivity  from  equation  (23*32).  In  view  of  the 
probable  length  of  calculation,  numerical  solution  of  equations  (23*33)  to  (23*36) 
to  obtain  the  optimum  conductivity  was  not  attempted* 

The  resulting  values  of  conductivity  for  helium -cesium  mixtures  are  shown 
in  Fig.  23*1  as  a  function  of  seeding  fraction  and  electron  temperature.  The 

conductivity  is  shown  as  the  ratio  of  the  conductivity  at  an  arbitary  seeding 

*•20  2 

fraction  to  that  of  the  pure  cesium*  Cross  sections  of  b01  x  10  m  ana 
320  x  10~c°  m2  were  used  for  electron-helium  and  electron-cesium  collisions 
respectively  (McNab3  and  Chapter  22  of  this  report).  The  calculations  were 
performed  for  constant  gas  temperature  (1000‘K)  and  pressure  t1  ata)  with  varying 
electron  temperatures  (1000°  to  5000°K.)c  Fig.  23  =  2  shows  the  magnitude  of  the 
optimum  and  pure  alkali  metal  conductivities  as  a  function  of  electron  temperature. 
Using  Fig.  23.1  and  23.2  the  actual  magnitude  of  the  scalar  conductivity  may 
easily  be  found  for  any  seeding  fraction. 

The  specific  power  of  a  segmented-electrcde  generator,  neglecting  ion 

2 

slip,  is  simply  proportional  to  cqB  ,  thus,  since  is  xndependant  of  magnetic 
field  (at  local  static  temperature  and  pressure  conditions),  the  specific  power 
curves  for  a  segmented-electrode  generator  as  a  function  of  seeding  fraction  are 
exactly  the  same  shape  but  are  displaced  by  an  amount  depending  on  B. 

2 

For  the  transverse  conductivity  (j^^  )  this  is  not  ^rue  since  Pe  is 
a  function  of  the  magnetic  field  strength*  Fig.  23*3  and  23*4  shew  oj  1+0 Q  as 
a  function  of  seeding  fraction  for  several  values  of  B  and  electron  temperatures 
of  2000°  and  50CD°K,  the  static  gas  temperature  and  pressure  being  the  same  as 

before . 

The  power  density  of  a  solid-elec troae  generator  is  proportion  1  to 

a  B2/1 +B  2*  this  is  shown  in  Fig.  23*5  and  23,6  for  the  same  conditions  as  in 

O  re  ' 


Fig.  23.3  and  23.4. 
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The  specif ic  power  of  a  Hall  generator  is  proportional  to  c^B  /l+£ 
and  this  is  shown  in  Fig.  23.7  and  23.8,  again  for  the  same  conditions. 

In  addition  to  the  dependence  of  specific  power  (p)  on  conductivity  and 
magnetic  field  there  is  also  a  gas  velocity  (u)  dependence: 

U2 

M2  rfa  1  .  (23.37) 

where  and  R^  are  the  ratio  of  specific  heats  and  gas  constant  for  the  mixture 
respectively.  Assuming  that  both  helium  and  cesium  are  ideal  monatomic  gases 
with  f  =  V3,  the  specific  power  for  constant  expansion  is  proportional  to  the 
gas  mixture  composition  only  through  the  R^  term.  The  mixture  gas  constant  is 
proportional  to  the  fraction  by  weight  (a)  of  each  constituent,  so  that,  in 
general 

=  }  a  ,R . 

Lj  J  0 


and,  for  helium-cesium  mixtures 


\  “  “CS^S  +  “He^e 

or,  if  a  is  the  fraction  by  weight  of  cesium, 

hi  *  0  RCS  +  (1^  hie 
Using  the  quantities  of  Section  23.4  it  can  be  shown  that  the  relation  between  a 
and  X  is 

\/ 

.  (23.39) 


(23.38) 


a  = 


X  ms 


*Z»s-ap)  ♦  *p 

Combining  equations  (23.38)  and  (23.39),  R^  may  be  evaluated  as  a  function  of 
seeding  f raction;  the  result  is  shown  in  Fig.  23.9.  Fig.  23.10  shows  the  product 
upon  which  the  power  density  of  a  segmented  electrode  generator  depends, 
while  Fig.  23.11  to  23*16  show  the  transverse  conductivity  and  specific  power 
dependence  for  solid-electrodo  and  Hall  generators,  including  the  effects  of  , 

23.5  DISCUSSION  AND  CONCLUSIONS 


Fig,  23.1  shows  the  ratio  of  mixture  conductivity  to  pure  seed  metal 
conductivity  (a/e^)  as  a  function  of  seeding  fraction  (X)  for  several  electron 
temperatures  in  cesium-helium  mixtures-;,  the  optimum  seed  fraction  increases  with 
electron  temperature  fairly  rapidly  above  airoyt  4000°K.  At  an  electron 
temperature  of  5000°K  the  optimum  seeding  fraction  is  very  high  (-'■0.8),  so  that, 
for  the  conditions  examined  here,  the  maximum  conductivity  is  obtained  with  puie 
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cesium.  As  the  electron  temperature  increases  the  variation  of  afo  with  X 
becomes  less  rapid  and  off -optimum  operation  has  leas  effect,  and  factors  other 


than  electrical  conductivity  (e.g.  engineering  problems 
could  determine  the  seed  fraction. 


or  the  cost  of  cesium) 


Fig,  2j ,2  shows  the  optimum  for  helium-cesium  and  pure  cesium  conductivities 
as  a  function  of  electron  temperature  the  steep  initial  rise  is  apparent. 

r\ 

Fig.  23-3  and  23*4  show  the  transverse  conductivity  (o^/l  */3  as  functions 

of  seeding  fraction  for  various  magnetic  field  strengths  (B) .  As  B  increases  the 

transverse  conductivitv  decreases  and.  for  T  =  2000°K  the  optimum  X  tends  to 

“  e 

unity  (pure  alkali  metal). 

P 

Figi  23„5  and  23»0  show  the  product  of  the  transverse  conductivity  and  B 

for  the  same  values  of  magnetic  field  as  previously.  Again,  as  B  increases, 

the  optimum  seeding  fraction  tends  to  unity  for  T  -  2GQ0°K-  and  for  T  =  5000°K 

$  e 

the  optimum  seeding  fraction  is  always  unity, 

n  o  o 

Fig.  23.7  and  23*8  show  the  term  a  B  as  functions  of  B  and  Xm 

u  e  e 

As  B  increases,  this  term  (which  is  proportional  to  the  specific  power  of  a  Hall 
generator)  increases  and  the  optimum  tends  to  unity.  The  optimum  seeding  fraction 
is  very  much  lower  than  in  the  previous  cases. 

As  stated  in  Section  23.4,  in  addition  to  the  dependence  of  specific  power 
for  an  MFD  generator  on  conductivity  and  magnetic  field,  there  is  also  a 
dependence  on  seeding  fraction  through  the  velocity  term.  This  can  be  reduced 
to  a  dependence  on  the  mixture  gas  constant,  R^,  which  is  given  as  a  function  of 
seeding  fraction  in  Fig,  23*9. 

For  a  segmented-elec-trode  generator,  neglecting  ion  slip  and  using 
equation  (23.37), 


f 


(23.40) 


where  f  =  P/iT^li^TKO-X)  (23.41) 

The  product  c Q is  shown  in  Fig.  23.10  as  a  function  of  seeding  fraction  and 

electron  temperature.  The  optimum  seeding  fraction  is  now  relatively  insensitive 

to  the  electron  temperature  (compare  with  Fig,  23-0  and  is  close  to  one  atomic 

percent.  Since  is  independent  of  B  (at  local  static-  temperature  and  pressure 

conditions)  the  normalized  specific  power  curves  \P  versus  x)  for  different 

magnetic  field  strengths  have  the  same  share  but  are  disc laced  by  amounts 
2 

depending  on  B  . 

The  manner  in  which  the  transverse  conductivity  (shown  in  Fig.  23.3  and  23*4) 
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is  modified  by  the  variation  of  gag  cons tan 

r~* 

where  is  giver*  as  a  function  of 

is  obtained  nith  smallest  magnetic  fields, 
optimum  seeding  fraction  tends  to  unity,  it 


it  is  shown  in  Fig.  23*11  and  23.12, 
X  and  B,  Again  the  maximum  value 
but.  although  as  B  increases  the 
:•  is  less  than  that  of  Fig.  23*3  and 


23.4. 


For  a  solid-electrode  generator  the  specific  power,  neglecting  ion  slip, 

p  =  ffoV2'/l+^e2  . (23.42) 

This  is  shown  in  Fig.  23*13  and  23.14  as  a  function  of  seeding  fraction  and 
magnetic  field  for  =  2000°  and  50QQ°K.  In  both  cases  P  and  the  optimum 
seeding  fraction  increase  as  B  increases,  the  optimum  being  in  the  approximate 
range  0.1  <  X  <  0.01 . 


For  a  Hall  generator,  again  neglecting  ion  slip,  the  specific  power  is 

t/.V 


p  = 


1*. 


.  (23.43) 


This  is  shown  in  Fig.  23*15  and  23.16  for  the  same  conditions  as  in  the  previous 
figures.  The  optimum  seeding  fraction  is  much  smaller  in  these  two  cases  than 
for  the  segmented-  or  solid-electrode  generators,  being  <0.0la/°. 

For  the  calculations  reported  here  =  ^/3  and  T  =  1000°K;  taking  M  =  0.5 
and  K  s  0.5  gives  jr^TKO-K)  »  102  ,  so  that  the  actual  specific  power  may  readily 
be  obtained  from  ‘P. 


In  any  practice. 1  MPD  generator  local  gas  parameters  will  vary  with  position 
in  the  channel;  the  gas  parameters  used  here  (T  =  1000°K  and  p  =  1  ata)  may 
represent  typical  values  for  a  large-scale  nuclear  MPD  generator.  The  techniques 
used  here  may  easily  be  adapted  to  give  the  best  seeding  fraction  for  maximum 
electrical  conductivity  or  specific  power  for  any  other  gas  mixture.  No 
mechanism  whereby  electron  temperature  elevation  is  achieved  has  been  postulated; 
different  methods  of  extrathermal  ionization  yield  different  dependence  of 
electron  temperature  on  gas  parameters. 

The  conclusions  found  here  relating  to  the  dependence  of  scalar  conductivity 
on  seeding  fraction  are  similar  to  those  quoted  by  Zimin  and  Popov^  except  that 
they  relate  to  elevated  electron  temperatures  and  low  gas  temperatures,  thus 
being  directly  applicable  to  closed-cycle  nuclear-MPD  generators. 
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24,1  INTRODUCTION 

Neglecting  Hall  effects  and  ion  slip,  the  specific  power  (power  output 
per  unit  volume,  P)  of  a  segmented  electrode  KID  generator  is 

P  =  oV^^Kd-K)  ....  (24.1) 

All  MFD  generators  built  or  under  construction  obtain  the  directed  (flow) 
motion  from  the  initial  thermal  energy  of  the  gas  by  expansion  through  a 
convergent  or  convergent-divergent  nozzle,  In  this  pi:.„ess  the  static 
temperature  and  pressure  fall  and  the  velocity  is  increased. 

When  thermal  equilibrium  exists  in  the  gas  the  scalar  electrical 
conductivity  depends  on  temperature  as 

F  (p,T)  exp  (-D/T)  ....  (24.2) 

where  F(p,T)  is  a  function  of  static  gas  „emperature  (T)  and  pressure  (p)  and 
D  is  constant  for  a  given  gas  (mixture).  The  temperature  dependence  of 
conductivity  in  the  range  of  interest  for  MFD  generation  is  determined  primarily 
by  the  exponential  term. 

Thus,  for  constant  magnetic  field  strength  and  loading  factor,  o  may  be 
maximized,  giving  maximum  specific  power.  This  calculation  has  been  performed 
by  several  workers  (for  example,  Sutton  )  and  the  results  indicate  that  the 
optimum  Mach  number  for  maximum  eV^  lies  below  unity  \0,6<k.  M  aC.3)  for  both 
closed  cycle  (seeded  inert  gas)  and  open  cycle  (seeded  combustion  products) 
devices. 

The  basic  limitations  on  the  power  obtainable  from  a  closed  cycle  nuclear 

heated  MPD  generator  is  the  working  temperature  of  the  reactor-  this  limits  the 

maximum  total  enthalpy  of  the  working  fluid  ant  hence  the  maximum  value  of 
2 

eV  .  The  working  temperatures  of  present  generation  inert  gas  cooled  reactors 
are  not  high  enough  to  permit  economic  KPD  power  to  be  generated  from  a  device 
working  at  thermal  equilibrium  conditions.  Future  developments  of  nuclear 


24.1 


reactors  to  yield  higher  working  temperatures.,  while  resulting  in  higher  MFD 
generator  outputs,  will  introduce  more  severe  materials  problems ,  capable 
only  of  long  term  and/or  expensive  solution.  This  paradox  may  be  resolved 
if  the  electrical  conductivity  of  the  working  fluid  can  be  increased  without 
appreciable  increase  of  the  gas  temperature;  high  specific  powers  ccr.ld  then 
be  obtained  without  severe  materials  problems. 


Several  methods  of  obtaining  thin  extra  thermal  ionisation  in  closed 
cycle  MFD  generators  have  been  suggested,  for  example:  high  voltage 

2  2  *  "o  ^ 

electron  beams  ,  high  current  arcs  ,  fission  product  P  rays*  ,  r.f,  , 

4  s 

magnetically  inauced  ionisation  and  photoionisation".  While  some  of  these 
methods  will  undoubtedly  prove  valuable,  none  has  yet  been  conclusively 
demonstrated  in  an  MFD  generator  of  substantial  power  output. 


>4.2  HON-wUILIBRIUH  FLOW 


A  farther  method  of  producing  extra  thermal  ionisation,  which  has 
received  little  attention  since  its  vi  iginal  suggestion  (independently)  by 
Lindley°  and  Sschenroeder  and  Daiber' ,  is  to  utilise  an  extremely  rapid 


expansion  of  the  working  fluid  through  a  large  area  ratio  nozzle  to  produce 
a  high  rate  of  fall  of  gas  temperature.  When  the  temperature  falls  rapidly 


in  a  gas  many  internal  adjustments  must  be  made  by  n aleouiar  collisions. 

For  example,  the  energy  of  molecular  vibrations  must  be  reduced,  chemical 
reactions  must  take  place  at  a  different  rate,  a  new  balance  must  be  found 
between  atoms  and  molecules,  and  ions  and  electrons  must  recombine  to  form 
neutral  particles.  All  these  adjustments  require  a  large  number  of  collisions 


between  particles  before  a  new  equilibrium  is  established.  If  the  temperature 
changes  more  rapidly  than  equilibrium  is  established,  non-equilibrium  flow 
results. 


*  :  +  r\ 
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Many  investigations,  both  theoretical  and  experimental  O'. 

non-equilibrium  flows  have  been  published  recently;  most  are  concerned 
primarily  with  molecular  gases  or  gas  mixtures  where  vibrational  ani 
dissociative  nun-equilibrium  are  the  dominant  factors.  This  inter?  t  arises 
becc „ne  the  usual  domains  in  which  such  flows  are  encountered  relate  to  real 


or  simulated  (shock  tube  and  wind  tunnel)  high  altitude  hypersot ic  and  missile 
re-entry  conditions  in  air,  and  similar  conditions  in  combustion  .gases  (for 
example,  ignition  in  hypersonic  ramjets).  In  general,  the  theoretical 
solutions  of  such  problems  are  complicated.  Because  the  changes  of  energy 
during  the  relaxation  processes  usually  form  a  significant  fraction  of  the 
total  enthalpy  of  the  gas,  the  quasi-one-dimcnsional  flow  equations  must  be 
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solved  simultaneously  with  the  rate  equations  which  determine  the  population 
densities  cf  the  (many)  species  present.  In  many  cases  the  required  rate 
coefficients  are  not  fully  known. 

The  flow  of  a  monatomic  recombining  gas  is  rather  more  simple  since 
valid  simplifying  assumptions  can  be  made  concerning  the  number  of  species 
present,  consequently  fewer  rate  equations  are  required  to  characterize  the 
flow.  For  example,  the  recombination  mechanism  may  be  assumed  to  be 
controlled  by  only  one  excited  atomic  state.  One  significant  difference 
between  a  recombining  atomic  gas  and  a  relaxing  molecular  gas  is  that  the 
electron  temperature  may  be  considerably  higher  than  the  temperature  of 
the  heavy  particles  in  the  atomic  gas.  This  leads  to  reduced  recombination, 
amongst  other  effects. 

In  order  to  examine  the  consequences  of  non  -equa librium  nozzle  flow  cn 
the  design  and  power  output  of  KPD  generators  the  following  investigation  lias 
been  carried  out:  while  later  investigations  may  show  that  this  non-equilibrium 
effect  is  insignificant  in  large  scale  devices,  it  will  undoubtedly  be 
important  in  analyzing  the  performance  of  small  devices,  such  as  that  at  IRD. 

24.3  TH30RY 

Th  e  mot'  common  concept  of  a  closed-cycle  nuclear  fuelled  MPD  generator 
envisages  the  working  fluid  to  be  an  inert  gas  seeded  with  small  amounts  of  an 
alkali  metal;  flow  through  the  expansion  nozzle  of  this  type  of  generator  will 
be  considered  here. 

The  assumption  of  small  fractional  seeding,  which  is  generally  valid, 
enables  considerable  simplification  to  be  made  in  the  analysis  since  the 
energy  involved  in  the  ionization  and  recombination  processes  in  the  seed 
element  form  only  a  small  fraction  of  the  total  enthalpy  of  the  working  fluid 
and  thus  may  be  neglected. 

Consequently,  neglecting  boundary  layers  and  shock  phenomena,  the  qunsi- 
one-dimensional  flow  equations  for  the  parent  gas  only  can  be  assumed  to 
approximate  adequately  the  flow  along  the  nozzle  axis.  An  additional 
simplification  resulting  from  the  assumption  of  small  fractional  seeding  is 
tnat  the  thermodynamic  properties  of  the  mixture  approximate  closely  those 
of  the  parent  gas,  ana  are  taken  as  such  in  this  analysis. 

The  flow  will  be  assumed  ' isentropic ' ,  so  that  it  can  be  described  by 
the  following  equations 


*  *  •  * 
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and,  for  a  process  with  an  isentropic  efficiency  q  =  /dT 


dT  f-1 
T  " 71  t 


&  =  0 
P 


....  (2k,1*) 
....  (24,5) 

....  (24.6) 


*  •  •  • 


(24.7) 


where  dT  is  the  actual  temperature  change,  dT'  is  the  isentropic  temperature 
change  and  r/  accounts  for  friction  effects  in  the  nozzle  expansion. 

For  a  given  nozzle  profilet 

A  =  b  f(x)  _  (24.8) 


Combination  of  equations  (24.3-8)  gives 


d  H 
dx 


TjVi  (1  +  M  )  d f(x) 

[  rK2  -  r?(y- 1)  -  rjl  f(x)  dx 


dT  _  r)Hd  T  (r-1)  d  f(x) _ 

dx  [  fM2  -  n  (r-D  M2  -  t]}  fix) 


....  (24.9) 


....  ( 24 .10) 


ox 


Solution  of  these  equations  in  conjunction  with  the  explicit  relations 

V  =  M  (rHT)^  ....  (24.11) 

and 


m 


P  =  bf(x)V 


....  (24,12) 


enables  the  flow  through  the  nozzle  to  be  described  if  the  nozzle  profile  and 
starting  conditions  are  known. 

To  complete  the  description  of  non-equilibrium  flow  two  further  equations 
are  required,  describing  the  rate  of  production  and  loss  of  electrons  and  the 
electron  energy  balance  respectively. 

The  most  general  form  of  the  electron  energy  equation  will  include  the 
energy  interchange  between  electrons  and  heavier  species  (atoms  and  ions)  by 
elastic  and  inelastic  collisions,  the  energy  gain  during  recombination,  the 
energy  lost  during  ionization,  radiation  and  thermal  conduction  in  the  electron 
gas.  */hile  equations  have  been  formulated  which  contain  all  these  terms  ,  they 
are  subject  to  considerable  doubt  regarding  their  accuracy  and  require  large 
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computers  for  their  r'lation,  Approximate  solutions  hare  been  found  inaccurate 

21 

in  certain  regions  (predicting  negative  temperatures^  ‘ } .  in  the  present 
analysis  the  electron  temperature  is  assumed  equal  to  the  gas  temperature 
throughout  the  expansion.  This  assumption  should  cause  the  electrical 
conductivity  at  the  nozzle  exit  (that  is,  inlet  to  the  KPD  generator)  to  be 
lower  than  the  true  value  since  higher  electron  temperatures  would  reduce  the 
rate  of  recombination  and  give  a  high  electron  density  and  conductivity. 

(This  will  be  partially  offset  by  the  increased  electron  collision  frequency.) 


The  equation  governing  the  rate  of  production  and  loss  of  electrons  is 
obtained  from  the  recombination  and  ionization  rates.  The  recombination 
mechanisms  considered  here  are  three  body  electron-electron-ion,  three  body 
electron-atom-ion  and  radiative: 

.  h  - 

A  +  e  ^  A  +  2e 

a 

1 


^2  + 

A  +  B  —Zk  +  e  +  B 

*2 

A  A+  +  e~ 


the  rate  coefficients  being  a1t  and  respectively. 

The  corresponding  ionization  coefficients  (/ 3^/9 ^  and  /9X)  are  obtained 
from  the  recombination  coefficients  using  the  condition  that  the  rates  must  be 
equal  at  equilibrium,  thus 


1 


etc. 


The  net  rate  of  production  of  electrons  per  unit  volume 

°  “  P  1nAne  +  ^2nA**B  +  ^3nA  “  *1ne  “  a2ne  "b  “  a3^e 


then  becomes 


s  *  %  ♦  ‘f*  *  “3: 


1 

iV) 

n. 

1  -  i 

L  i 

\nk  i 

K  11 

6  e 

....  (24.13) 


The  net  rate  of  production  of  electrons  per  unit  volume  in  a  volume  V  containing 

N  (=  r.  V  )  electrons  is 

e  e 

„  1  dNe  1  dtneV1) 


V*  dt 


r1  dt 
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For  a  moving  system  and  when  the  fractional  ionization  and  seeding  fraction  are 
small |  this  becomes 

d(vp)  x„v 

S  -  p 


dt 


dX 

dx 


Using  this  expression  for  S equation  (24.13)  becomes 


a  n  ^  /n  ^ 

Si;  -  (an  +an+o)  ( 

dx  p  VX  k  1ne  2"B  3  \  n. 


n. 


1-  '  “A  •  s  nt2 


-  1 


....  (29.14) 


-  (24.15) 


Simultaneous  solution  of  equations  £4, $,(24.10) and  (24.1 5\  in  conjunction  with 
equations  £4.1l\ (24.12)  and  the  relation 


n  mX  pV/mp 


....  (24.16) 


enables  the  flow  through  the  nozzle  and  variation  of  electron  density  to  be 
described  when  the  nozzle  area  is  specified  as  a  function  of  distance  and  the 
starting  conditions  are  known  (See  Appendix  24. A). 


The  scalar  electrical  conductivity  may  then  be  obtained  from 

>2 


n  e£ 
e 

m  v 
e  e 


....  (24,1?) 


where 

v  =  v  £,  q  .  n. 
e  e  “3  4ej  j 

and  the  cross  sections  employed  are 

Vh.  *  6'1-10'20  «2 
%-c.  *  io'20  «2 


1,84.10”^  De  In 

- 37T-  € 


|j.24. 


10?tV2 


....  (24.18) 


from  references  22,  25  and  24  respectively. 

This  value  of  conductivity  may  then  be  compared  with  the  equilibrium  value 
of  conduct Jvity  obtained  when  n^  is  determined  from  the  Saha  equation; 


2s.)  - 

n,  / 


A  's 


CT  exp  (-  D/T) 


....  (24.19) 


where 


«  .  3/2 

2vm  k  g  g. 

-I-  _|_i  ,  D  -  ,VlA 


are  constant  for  a  given  gas  (mixture). 
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24,4  CALCULATIONS 


The  theory  outlined  above  has  been  used  to  construct  a  computer 
programme  to  determine  the  electrical  conductivity  through  a  nozzle.  The 
calculations  have  been  performed  for  the  operating  conditions  of  the  MFD 
generator  under  construction  at  ltd).  The  working  fluid  is  helium  seeded 
with  one  atomic  percent  of  cesium;  the  appropriate  data  are  given  in 
Table  24.1. 

The  IRD  MF’D  generator  is  designed  to  operate  supersonically  and 
consequently  has  a  convergent-divergent  nozzle  of  which  two  sides  are 
parallel  and  two  are  axially  symmetric.  In  equation  (24,8),  b  is  twice  the 
width  of  the  nozzle  and  f(x)  is  the  height  from  the  centre  line.  The  nozzle 
profile  (rig24-.l)  is  composed  of  three  regions,  being  cosinusoidal  at  inlet 
and  exit  and  circular  near  the  throat: 


<a. 

i 

f(x) 

=  b1 

COS 

(b2x)  + 

i  a 

....  (24.20) 

f  (x) 

b4 

“  i 

b..  -  4(x‘  -b-.x  + 
-5  o 

yj* 

....  (24.21) 

>a2  * 

f(x) 

=  b8 

COS 

+  b*i 

....  (24.22) 

Difficulties  have  been  encountered  in  the  solution  of  the  equations  because  of 
this  profile:  where  the  curves  join,  discontinuities  in  the  second  derivative 
exist  and  these  cause  unsmooth  variation  in  the  theoretical  flow,  for  example  in 
the  derivative  of  the  Mach  number,  temperature,  etc.  In  practice,  such 
variations  will  probably  be  smoothed  out  by  boundary  layers. 

Additional  difficulties  are  encountered  in  the  solution  of  the  differential 
equations  for  dM/dx  and  dT/dx  since  these  become  indeterminate  near  M=1  (for  r?=1). 
These  have  been  overcome  by  the  use  of  equation  (24.A6)  near  the  throat,  Mach 
number  is  substituted  into  this,  giving  the  nozzle  area,  use  of  equations  (24.8' 
and  (24,21)  tren  enables  nozzle  distance  to  be  obtained.  Thus  all  the 
thermodynamic  and  flow  parameters  are  known  at  each  point  along  the  throat. 

Initial  versions  of  the  prograovne  extrapolated  electron  density  across  this 
throat  region,  however  it  was  later  found  more  accurate  to  extrapolate  fractional 
ionization  (so  that  dX/dx  remains  constant).  The  region  over  which  this 
extrapolation  is  performed  is  a  small  fraction  of  the  nozzle  length  and 
appreciable  error  is  unlikely  to  be  incurred. 

The  isentropic  efficiency  in  the  present  calculations  is  assumed  to  be 
unity  to  the  throat  and  0.9  afterwards. 
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24.5  RESULTS  AND  DISCUSSION 


Figures  24.2  to  24.8  show  typical  results  obtained  for  two  sets  of  data 

(T  =  1500°  and  2000°K,  p  =  1  atm)  in  a  helium-cesium  mixture*  (Other  data 
o  ro 

for  these  results  being  given  in  Table  24.1.) 

Figure  24.2  shows  the  variation  of  Mach  number  through  the  supersonic 
nozzle;  as  equation  (24.A6)  indicates,  the  Mach  number  depends  only  on  the 
profile  of  the  nozzle  provided  the  ideal  conditions  are  obeyed.  As  with 
most  of  the  results  given  here,  rapid  variation  of  the  Mach  number  in  the 
region  of  the  nozzle  throat  (x  -  1.905  cm)  is  evident. 

Figure  24.3  shows  the  variation  of  static  gas  temperature  through  the 
nozzle  for  the  two  sets  of  data  presented  here.  For  both  these  cases  the 
static  gas  pressure  at  the  nozzle  exit  is  0.05  atm,  consequently,  if 
equilibrium  conditions  (as  given  by  Saha's  equation*-  at  the  gas  temperature) 
were  obeyed,  the  electron  concentration  at  the  nozzle  exit  would  be  very  low. 
That  is,  the  conditions  at  the  generator  inlet  would  be  far  from  the  optimum 
(found  by  maximizing  <rV^)  mentioned  in  Section  24.1.  The  non-equilibrium 
electron  concentrations  (found  from  equation (24. 15)  and (24, 16))  is  shown  in 
Figure  24.4,  and  while  the  concentrations  drop  on  expansion,  in  both  cases 
the  nozzle  exit  value  is  not  more  than  an  order  of  magnitude  less  than  the 
inlet  value. 

Figures  24.5  and  24.6  compare  the  equilibrium  (as  given  by  Saha)  and 
non-equilibrium  (as  found  using  the  values  shown  in  Fig.  24.4)  electrical 
conductivities  as  a  function  of  the  distance  through  the  nozzle  for  T  =  1500° 

Q  O 

and  2000  K  respectively.  (Note  the  ordinate  scale  in  both  cases.)  The 
non-equilibrium  conductivity  for  Tq  =  1500°K  exhibits  a  slight  but  definite 
rise  through  the  nozzle.  This  occurs  because,  even  though  the  electron 
density  decreases,  the  opposing  effect  caused  by  the  decrease  in  temperature 
has  a  greater  effect.  In  Fig.  24.6  the  two  effects  approximately  cancel 
and  the  non-equilibrium  electrical  conductivity  is  sensibly  constant.  In 
both  cases  the  equilibrium  conductivity  drops  rapidly  (near  and  after  the 
throat  region)  to  very  low  values. 

The  effects  of  the  non-equilibrium  conductivity  on  the  power  output  of 
an  MPD  generator  are  shown  in  Figs.  24.7  and  24.8  where  the  quantities 
are  shown  as  a  function  of  Mach  number  in  the  nozzle.  For  a  segmented 
electrode  generator,  neglecting  Hall  and  ion  slip  effects 

P 

B^Kd-K) 
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hence  if  s  =  d  webers/m  and  K  =  0.5,  the  quantity  -rt~  is  numerically  equal 
to  the  specific  power  (in  MWe/ffi  ). 

For  comparison  the  specific  powers  obtainable  using  equilibrium 
conditions  are  shown  in  both  Figures.  These  equilibrium  power  densities 
exhibit  the  previously  mentioned  maximum  at  a  Mach  number  of  about  0,5. 

The  non-equilibrium  power  densities,  however,  increase  monotonically  with 
Mach  number  o  values  greater  than  10^  and  10'  MWe/m*?  for  T  =  1500°  and 
2000  K  respectively.  The  ratio  of  non-equilibrium  to  equilibrium  specific 

IQ  H 

power  at  the  nozzle  exit  in  the  two  cases  are  about  10  and  10  .  A  more 
practical  comparison  of  equilibrium  and  non-equilibrium  specific  power  at 
nozzle  exit,  with  the  equilibrium  specific  power  at  the  optimum  Mach  number; 
these  ratios  are  about  45  and  20  for  Tq  =  1500U  and  2000  K  respectively,  that 
is,  they  are  worthwhile  increases.  Comparison  of  Figs,  24,7  and  24,8  shows 
that  the  non-equilibrium  power  density  at  nozzle  exit  for  =  1500°K  is 
approximately  the  same  as  the  equilibrium  power  density  at  the  optimum  Mach 
number  for  Tq  =  2000°K.  That  is,  operation  of  a  nozzle  to  give  semi-frozen 
flew  (as  far  as  electrical  conductivity  is  concerned)  permits  the  same  power 
density  to  be  achieved  at  a  stagnation  temperature  about  500°K  lower  than  for 
the  equilibrium  case  (for  these  particular  conditions). 

The  recombination  coefficients  used  in  the  evaluation  of  equation  (24.15) 


three  body  electron-electron-ion 


a 


1 


1.1.10 
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three  body  electron-atom-ion  (Thomson  coefficient) 


a 


2 


and  the  radiative  recombination  coefficient 


u3 


—lb  — X/4 

=  3.10  T 


Of  these,  for  the  conditions  examined  in  the  present  calculations  has  the 

l 

greatest  effect,  thus  the  energy  loss  by  radiation  implied  by  the  inclusion 
of  the  radiative  process  is  not  significant. 


In  the  immediate  future  further  calculations  will  be  performed  using  the 
present  computer  programme,  with  slight  modifications  (for  example  using  the 
plasma  mixture  gas  constant).  It  is  hoped  that  a  more  refined  calculation 
including  tne  electron  energy  balance  will  eventually  be  undertaken,  although, 
as  mentioned  in  Section  24.3,  the  calculations  wxll  be  more  complicated  tnan 
present  ones. 
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APPENDIX  24. A 


STARTING  CONDITIONS 


The  procedure  used  here  to  obtain  the  starting  values  for  the  nozzle  is 
as  follows: 


From  the  equations  of  state 

p  =  pRT  , 

....  (24. A1) 

continuity 

m  =  P  AV  , 

-  (24. A2) 

enthalpy 

T  ■  T  +  V2  , 

0  2Cp 

....  (24. A3) 

and  Mach  number 

M2*  V2 

J'HT 

....  (24.A4) 

the  criterion  for  M=1  to  be  achieved  at  the  throat  of  a  convergent-divergent 
supersonic  nozzle  working  isentropically  is 
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(24. A5) 


Thus  for  a  given  nozzle  throat  area  and  given  gas,  selection  of  the  inlet 
stagnation  temperature  and  pressure  detmines  the  mass  flow, 
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From  the  relation  given  by  Shapiro  for  the  isentropic  flow  of  a  perfect 

gas 


A 

A* 
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X  +  1 


4  +1 
2  (  y-1) 


....  (24.A6) 


The  inlet  Mach  number  may  be  obtained  from  the  nozzle  inlet  area. 

Knowing  the  inlet  Mach  number  the  static  temperature  and  pressure  may 
easily  be  determined  from  the  above  equations,  and  the  remainder  of  the  entrance 
conditions  follow. 


TABLE  24.1 


DATA  USED  IN  CALCULATIONS  (MKS  UNITS) 

»1  =  1,635592232.1 O'2 
a2  =  1. 98069310. 1G~2 
b1  =  9.46798320. 10"^ 
b2  =  1. 3962884. 10? 
b?  *  9.5820168. 10"' 
bk  =  6.0833. 10"3 
bc  •  1. -.8025156. 10*"4 
bg  =  3.81.IO 

b?  «  3.85392939. 10"4 
bg  =  3.99766848. 10”3 
b9  =  1.28817056 
b10  =  2.54. 10“2 

=  2. 35233152. 10” 3 

b  =  0.0762 

C  =  2,415. 102‘ 

D  =  4,518. 104 

CD  «  5.2285.103 

R  =  2,075. 103 

t  =  1.658 

*  =  10"2 

Tj  =  0,9  for  £>0,0195 


NOZZLE  PROFILE  FOR  IRD  GENERATOR 
(ASSEMBLY  l) 
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VARIATION  OF  ELECTRICAL  CONDUCTIVITY  WITH  DISTANCE  THROUGH 
NOZZLE  FOR  EQUILIBRIUM  AND  NON- EQUILIBRIUM  FLOWS 
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VARIATION  Of  ELECTRICAL  CONDUCTIVITY  WITH  DISTANCE 
THROUGH  NOZZLE  FOR  EQUILIBRIUM  AND  NON- EQUILIBRIUM 
FLOWS  (T.«  2000*K) 


VARIATION  OF  <rV*  WITH  MACH  NUMBER  FOR  EQULBRLfM 
ANO  NON-  EQUILIBRIUM  FLOWS  (i;  s  I500°K) 
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VARIATION  OF  (TV1  WITH  MACH  NUMBER  FOR  EQUUBRJUM 
AND  NON-EQLALBRUM  FLOW  (\-  200CTK) 


NOMENCLATURE 
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A* 


1.2 


1-11 


c, 
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5e,i,a 

h 

K 

k 

M 
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P 

P 

q 

R 

s 

T 

T 

V 
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nozzle  cross  sectional  area 
sonic  tnroat  area 

see  equations  (24,20^  (24, 2i)  and  (24.2 2) 

ma/gietic  field  strength 

twice  nozzle  breadth 

see  equations  (24,20^  (24, 2 1)  and  (24.22) 

constant  in  Saha  equation  (24.19) 

specific  heat  at  constant  pressure 

eVi/k 

electronic  charge 
F(p,T)»  see  equation  24,2 
f(x/  -  half  nozzle  height 
atomic  statistical  weight 
Planck’s  constant 
loading  factor 
Boltzmann's  constant 
Mach  number 
particle  mass 
particle  concentration 
specific  power 
static  pressure 
collision  cross  section 
gas  constant 

net  rate  or  production  of  electrons 
static  gas  temperature 
stagnation  temperature 
flow  velocity 
ionization  potential 
random  particle  velocity 
fractional  ionization 


a  recombination  coefficient 

ionization  coefficient 
f  specific  heats  ratio 

ct  electrical  conductivity 

p  gas  density 
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n  isentropic  efficiency 

Af  seeding  fraction 

v  collision  frequency 

Subscripts 
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seed  gas  (cesium) 
parent  gas  (helium) 
denotes  Saha  equilibrium 
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‘5.1  INTRODUCTION 


An  investigation  of  some 


e  01 


riH  processes 


cur;  ing  in  a  aagne toplasmadynamic 
(MPD)  channel  to  assess  their  effect  on  flow  ;  ammeters  has  been  attempted, 
quantitative  estimates  being  obtained  for  the  IED  e* ceriment.  The  first  test 

-4 

section  of  the  IED  generator  (described  fully  elsewhere ‘ )  differs  from  other 
KPD  devices  in  several  respects:  the  working  fluid  is  very  pure  (f<  30  ppm  total 
gaseous  impurity)  inert  gas  (helium)  seeded  with  p  to  three  atomic  percent  of 
an  alkali  metal  (cesium);  the  accelerating  nc-ssie  is  de~igned  for  a  high  (20:1 ) 

(Mach  ; 
art 


issure  rat: 


l w  ?  §2. v ing  v ery  ni gn 
exit;  and  the  generator  channel  is  of  con 


f.'“  f1  : 


rcer  -  .  5 ;  at  the  nozzle 

with  r -enented  electrodes. 


Althongn  the  problems  under  investigation  here  refer  to  flow  in  the  channel, 

consideration  must  be  given  to  nozzle  processes  since  these  determine  the  inlet 

conditions  to  the  channel.  Thermodynamic  equi  ibrium  is  assumed  at  inlet  to 

the  nozz tej  that  is ,  the  components  of  the  gas  r  xture  *,_  alien  gas  atoms,  seed 

atoms,  ion3  and  electrons)  are  uniformly  dispersed,  have  the  same  temperature, 

and  the  fractional  ionisation  of  the  alkali  metal  is  determined  by  the  Saha 
2 

relation  .  It  is  also  assumed  that  there  is  no  fringing  of  the  magnetic  field 
into  the  nozzle. 

In  the  absence  of  a  magnetic  field,  electrons  and  alkali  metal  ions  and 
atoms  have  negligible  effect  on  the  bulk  flow  processes  since  the  seeding 
fraction  is  small  (generally  less  than  1%).  The  solution  of  the  flow 
equations  for  the  expansion  nozzle  yields  expressions  for,  inter  alia,  the 
temperature  and  pressure  at  tne  exit  I rom  the  nozzle.  From  these  two  parameters, 


for  a  given  seeding  material,  the  electron  density  may  be 
thermal  equilibrium  ionization  as  given  by  the  Saha  aqu  ft 
to  the*  extremely  rapid  expansion  of  the  plasma  in  ‘he  nos 
non-equilibrium  effects  may  occur.  The  first  of 
electron  temperature  at  the  exit  from  the  nosz.e. 


determine a ,  assuming 
n.  however,  owing 
-  u  at  least  two 
the  as  it  an  above -equilibrium 
Th Is  effect,  has  been 


C. 


i 


previously  examined^  and  found  to  be  negligible  for  the  oon-di tions  of  the  IRD 
generator*  Secondly*  owing  to  the  rapid  e xpansion  of  the  plasma*  recombination 
may  not  take  place  rapidly  enough  to  assure  the  thermal  equilibrium  electron 
concentration.  Under  these  conditions  the  flow  is  said  to  be  f semi -il osen*  and 
the  electron  c .noentration  at  the  exit  of  the  nosale  is  considerably  higher  than 
that  given  by  application  of  the  Saha  equation.  (The  electron  temperature, 
however,  is  the  same  as  the  gas  temperature).  The  thermodynamic  parameters  found 
from  solution  of  the  flow  equation  for  the  nozzle  see  Chapter  24,  together  with 
the  electron  concentration  for  the  semi-frozen  flow  uase,  are  used  as  the  inlet 
conditions  to  the  channel  for  the  present  evaluation. 

In  the  constant  cross-sectional  area  generator  channel  the  flow  equations 
are  modified  by  the  if  forces  arising  from  interaction  with  the  magnetic  field. 
These  modifications  consist  of  a  term  in  the  momentum  equation  for  the  magnetic 
decelerating  force,  and  a  term  in  the  energy  equation  for  the  electrical  power 
extracted.  Both  these  terms  contain  the  electrical  conductivity  and,  hence, 
the  electron  concentration.  Because  of  this,  in  contrast  to  the  nozzle, 

(which  has  no  applied  magnetic  field),  the  electron  concentration  may  have 
considerable  effect  on  the  flow  parameters. 

The  electron  concentration  may  be  governed  by  one  of  several  mechanisms; 
two  only  are  considered  here.  The  first  is  thermal  ionization  a *■  the  gas 

C 

temperature  and  pressure ;  in  the  second,  postulated  by  Kerreorosk”1’,  Hurwitz  et  al  , 

/  —  g 

Ken  .  Brother ‘  and  Wright  and  Swift-Hock  .  the  els  toons  have  a  higher  temperature 
than  the  remainder  of  the  gas  by  virtue  of  the  induced  electric  field-  Replacing 
the  gas  temperature  in  Saha’s  equation  by  this  electron  temperature  yields  an 

Q 

electron  density  greater  than  the  equilibrium  value,  BenDar'^1  and  Tamor^  have 
shown,  by  consideration  of  the  rate  equations  for  ionization,  excitation, 
recombination  and  de-excitation  ex’  cesium,  that  the  use  of  the  electron  tempera  tun 
in  Saha?3  equation  is  valid  for  MPD  generator  conditions. 

In  Section  25.2  (largely  based  on  MoNab  and  Cooper1  )  th^  flow  equations 

are  solved  assuming  ionization  equilibrium  a"  an  elevated  electron  temperature 

5 

which  is  obtained  by  the  methods  outlined  in  Hurwitz  et  al"h  Section  25.3 
details  more  recent  calculations  using  equations  giving  the  rate  of  change  of 
electron  temperature  and  electron  concentration, 

25.2  IONIZATION  EQUILIBRIUM  AT  THS  ELECTRON  TEFi -ERATURE 

25.2,1  Solution  of  the  flow  eguatiors 

In  the  solution  of  the  flow  equations  the  following  assumptions  have  been 

made: 


25.2 


(i) 


>n 


id 

ji 

und 

:h 


ore 

ing 


B 


ture 


the  flow  is  steady-state  and  one-dimensional  ii.e.  boundary  layers 
are  small  compared  with  the  duct  width;  reverse  currents  which  are 
thought  to  occur  in  the  boundary  layer  are  neglected) ; 

(ii)  the  properties  of  the  bounding  surfaces  are  not  considered  except 
for  the  estimation  of  friction  and  heat  loss  at  the  walls,  both  of 
which  are  assumed  uniform  along  1:  channel; 

(iii)  the  power  output  of  the  generator  is  determined  by  processes  in  the 
gas  only  (for  example,  possible  current  saturation  due  to  thermionic 
emission  characteristics  of  the  electrode  material  is  ignored); 

(iv)  end  loop  losses  are  neglected; 

(v)  the  energy  r eleased  or  absorbed  in  the  association  or  dissociation 
of  the  seed  material,  together  with  that  released  if  r eeombination 
occurs,  is  neglected; 

(vi)  impurities,  which  generally  reduce  the  electron  density  and  hence 
power  output,  are  absent. 

The  validity  of  these  assumptions  is  discussed  in  Section  25.2,5. 

25.2.1*1  Non-zero  magnetic  field 

For  the  system  shown  in  Fig.  25*1  the  following  equations  apply  (the 
symbols  being  explained  in  the  Nomenclature 

Energy 

Change  in  total  energy  =  mCjiT 


This  is  composed  of: 


(a)  heat  losses  to  walls 


4 h  (^aw~xw)  Adx 


(b)  net  electrical  power  extracted  =  ~E, J  Adx  -=  -cB2V^K(' -K)Ad 


hence 


mCpdTe 


=  -sbVkO-X)  ix  -  fa 


*  *  * » *  (25.1 ) 


Momentum 


Rate  of  change  of  momentum  =  mdV 
This  is  due  to: 

(a)  force  due  to  pressure  change  =  Adp 


25.3 


(b)  force  due  to  friction 

(c)  force  due  to  MFD  interaction 


2fpY2Adx 

D 

-  J  x  B  Adx  =  -oB* V(l -K)Adx 


hence 


Continujl 


=  _dp  -  Z&yffc  -  oB^V(l  -K)dx 


.....  (25.2) 


p  -  oR T 
Enthalpy 

V2 

T  =  T  +  -L. 
o  2C 


.....  (23.3) 


.....  (25.4-) 


(25.5) 


Prom  equation  25.1  to  25*5  an  expression  may  be  obtained  giving  the  stagnation 
tempera  ure  as  a  function  of  th«  gas  velocity  and  overall  parameters  of  the  system 
only  (see  Appendix  25.A,  equation  25.A.7}. 

The  velocity  may  be  obtained  as  a  function  of  distance  in  the  following 
Banner.  Using  equations  (25.1 )  and  (25.A.2)  to  eliminate  dTo,the  following 
expression  is  obtained: 

» ■  5^  . <„.« 

vOt^-Kg) 

On  integration,  assuming  Tq  constant,  this  yields: 

jr 

x  +  constant  =  In  [ (V2  -  Kg)  ^2  .  V  13]  .....  (25.7) 

K,T  E_ 

where  K.  =  - 2.  _ 

12  2K3  2Kp 

K  T 

and  K  ,  =  —  n--° 

13  2K3 

The  constant  in  this  equation  is  evaluated  from  the  initial  conditions. 

The  assumption  of  constant  Tq  is  valid  if  the  integration  is  performed  over 

small  velocity  intervals  since  Tq  varies  slowly  compared  with  the  other  flow 
parameters. 

Having  obtained  Tq  as  a  function  of  velocity  and  distance, all  the  remaining 
thermodynamic  quantities  may  be  obtained  from  equations  (25.3),  (25.4)  and  (25.5), 


1 

V  13i 


(25.7) 


25.4 


together  with  the  definition  of  the  Mach  number: 


St !  _  JL 
18  rRT 


25.2.1  2  Zero  magnetic  field 

When  BsO,  the  functions  /*  (V)  and  g!(V)  in  equa' 


f\W  =  -:Vj 

r-Kg) 


.....  (25.8) 


(ii'jiAiU’j  become 


g:  (v)  =  — 2-l — 

qlKj'V^-Kg) 


where  *  =  2fnyAd,  the  other  symbols  having  been  previously  defined. 

Equation  (25»A.4)  may  then  be  integrated  in  the  same  manner  as  before,  using  an 
integrating  factor.  In  this  case  the  integration  say  be  performed  directly, 


without  a  binomial  expansion,  to  give: 


^•(V-tKg') 


[in  V  +  (V^KgO^j 


+  constant 


.  (25.9) 


where 


[TawwTw)(r-1 

JTfmC 


From  equations  (25. S)  and  (25.7) ,  T  and  the  remaining  thermodynamic  quantities 

o 

may  be  obtained  as  a  f  unction  of  distance  for  the  case  when  3-0. 

25*2.2  Electrical  conductivity 


The  scalar  electrical  conductivity  of  a  plasma  is 


a  =  n  m  +  n.ee, 
e  i 


(25.10) 


according  to  simple  theory.  In  most  cases  the  contribution  fh'om  the  ions 
may  be  neglected  because  of  their  large  mass  and  consequent  low  mobility 
compared  with  that  of  the  electrons.  Thus  equation  (25.10)  becomes 


n  e 
e 

a  =  n  eu  =  - 

e  6  Vtot 


—  2  «8lq,  x  1 0 


.  (25. 11) 


where  n  is  the  electron  concentration  (per  &)  and  v__,_  is  the  total  electron 
e  f  TOT 

collision  frequency  for  momentum  transfer.  For  a  gas  mixture  containing  j 


atomic  species 


vmnm  -  2,v  . 
TOT  j  ej 


.....  (25.12) 


hence  for  a  helium-cesium  plasma  under  the  conditions  of  interest  for  KPD 
generators  (i.e.  no  Cs,-,  or  H  +): 

b  6 

-■’*£  c; 


'tot  ■  V-c.*  +'e-C»  *  Vh.  .  («.i3) 

(Note  that  there  is  no  f  term  in  this  expression,  because  electron-electron 

e-e 

soliiaions  have  negligible  effect  on  the  total  electron  momentum*) 

The  collision  frequency  is  given  in  terms  of  a  tomic  cross  section  by 

(8kT  *^V£  1 

B  9,.  n 4  s  jr—2,  I  Q  .  n.  =  6.211  *  105  X  2  Q  ,n  sec'1 
®«  •  1  J  ej  j  e  ej  J 

•  •  •  •  •  (23  *14) 

(in  mks  units)  for  ft  Maxwellian  velocity  distribution. 

Pres  a  previous  assessment  of  experimentally-obtained  low  energy  cross- 

sections  a  constant  value  of  6.1  x  10  2Sn2  has  been  used  for  Q  .  to  rive 

e-He  * 

Vk,  =  5'7?  *  10'16  t/"h,  . v«.i5) 

Little  data  is  available  for  the  electron-cesium  atom  collision  cross-section 
at  low  energies.  Erode12  gives  Q  ^  =  3*1  *  10* 5 8  m2  at  0.40  eV,  and  this 

figure  has  been  used,  assumed  constant,  to  give 

A 

* 

%-C,  ■  ’•»  *  1°'14  \2  . (25.16) 

The  equivalent  cross-section  for  electron-ion  collisions  is  given  by  Spitzer1^: 

4  —10 

Qei  =  25Sr(e  kT  )?  k  A  s  10 -  In  A  m2  . (25.17) 

0  6  T 

9 

T  3/2 

wh.«  A  =  1 .24-  x  107  J,/2  .  (25.18) 

e 

*“  *  Itvt121  *.3/2  ]  . (».i») 

«  © 

Equations  (25*11 ),  (25*13),  (25.15),  (25*16)  and  (25*19)  may  be  used  to  obtain 
the  conductivity  of  the  gas  when  the  electron  concentration  is  known;  n  is 
evaluated  for  two  different  cases  below. 

25.2.2,1  Equilibrium  at  the  gas  temperature 

The  electron  concentration  is: 

re  =  ^  =  7*24  *  1q22  tT  m"3  . (25*20) 

p  being  the  partial  pressure  of  the  cesium  vapour  in  newton,  m  .  The  fractional 
ionization  (X)  for  cesium  obtained  from  Saha’s  equation2  at  the  gas  temperature, 
is  given  by: 


25*6 


1.956  x  10" 


/'  O  C  *5 

...  -j  ,s£  1  } 


25.2.2,2  Equilibrium  at  the  electron  temperature 

v 

According  to  Hurwit-s  ©t  al^,  the  electron  field  induced  by  the  sags© tie 
field  in  an  MPD  generator  channel  will  cause  an  elevated  electron  temperature . 
For  a  segmented  electrode  generator  the  relative  temperatures  when  equilibrium 


is  obtained  are: 


o  9  9 


»  1  - 


...  (25.22) 


where 


1  .  7n  x  10  B/ 1 


.....  (25.23) 


In  equation  (25.22).  8  is  a  facto 


or  account!: 


K  i.-va.-wuU  :  fc- 


ineissti-o  collisions, 


deviations  from  the  Maxwellian  distribution  function,  and  nor -constant  collision 


frequency.  For  a  monatomic  gas  8  <%  1  and  an  elevated  electron  temperature  i 
fairly  readily  achieved.  (In  combustion  gases  0  may  be  1(f  or  10  ,  owing  to 
inelastic  collisions,  and  small  electron  temperature  elevation  is  likely.) 

The  electron  temperature  is  obtained  from  Saha's  equation  when  the  gas 


p  7~:  n  i  fm-n 


iron  temperature  is 


temperature  is  r enlaced  bv  the  electron  temneratui 


According  to  BenDsniel 


j  -  \ 

and  Tamor  ,  this  is  valid  for  MPD  generators.  Equation  (25.22)  is  actually  a 
simplified  version  of  Sutton’s  equation. 

25.2.3  Computations 


Equation  (25. A. 7)  is  obtained  by  integration,  assuming  e  to  be  constant 
over  the  velocity  interval  (V).  ,  a  mean  value,  o.  v-eing  used.  The  temperature 
and  pressure  obtained  as  a  result  of  this  intexrati  Oil  C*  Z'lt  li5  Sti  Xii  uf>6 


caxcuxatxcn  ox  a  new  vaxue  of  tne  conauotiv; 


end  of  the  velocity 


interval. _  The  process  is  then  repeated  using  a  new  value  of  conductivity 

(p  —  J/T+  ^n  1  f  oe.  A  Wl  A/>4-_-.„  t  . 


—  JfT" 


;  in  equation  (25.A.7) •  After  several  iterations. 


1+1  1  1+1  s+1 

o  the  values  of  ?,  u  p,  etc.  constitute  the  solution  at  the  end  of  the 
n+  * 

velocity  interval. 


The  electrical  conductivity  is  a  function  of  the  electron  concentration 
and  the  electron  temperature.  In  the  case  considered  in  section  25.2.2.1 
(Tp«T^ne  could  be  found  directly  using  equations  (25*20)  and  (25.21 ).  However, 
for  the  case  of  elevated  electron  temperature  (section  25.2.2.2)  an  iterative 
solution  for  n  and  T  we  necessary,  taking  place  within  that  described  above 
for  a.  This  consisted  of  solving  simultaneously  the  equation  obtained  by 
combining  equations  (25.20)  and  (25.21 )  (T  replaced  by  ?_).  which  is  of  * 


the  form: 


%  *  *  (*) 

and  an  equation  of  the  fora 

Te  ’  *  (E@)  *  0>) 

resulting  from  the  combination  of  equations  (25*13),  (25*22)  and  (25*23)* 

Since  (b)  is  of  the  fora  represented  in  Fig,  25*2 ,  difficulty  was  experienced 
in  obtaining  the  solution  when  it  lay  on  the  almost  vertical  part  of  the  curve. 

The  procedure  finally  used,  although  cumbersome,  did  not  give  rise  to  the 
difficulties  experienced  in  attempting  other  methods  (in  some  cases  the  solution 
oscillated  between  two,  or  even  three,  values)* 

Equations  (a)  and  (b)  were  combined  to  give  an  implicit  function  of  T 

6 


only: 

t?  (Te)  =  0  .....  (25.24) 

and  increasing  values  of  T  were  substituted  in  ?j( T  )  until  it  became  zero* 

J+l  I’tl  ^  ® 

Thus  T  and  0  n  were  obtained,  so  that  a  could  be  calculated.  A  flow 
0  © 

diagram  for  the  procedure  is  shown  in  Fig.  25.3* 


The  input  values  and  constants  used  in  this  calculation  are  given  in 
Tables  25.1  and  25*2. 

25.2.4  Results 

Investigations  were  made  for  three  conditions;  zero  magnetic  field; 
equilibrium  at  the  gas  temperature;  and  elevated  electron  temperature.  The 
first  case,  being  that  of  no  power  extraction,  consists  of  the  solution  of 
the  equations  of  supersonic  flow  in  a  constant  area  duct,  including  friction 
and  heat  transfer.  With  power  extraction  for  both  quilibrium  and  elevated 
electron  temperature,  the  effect  of  varying  the  magnetic  field,  B,  and  the 
load  factor,  K,  was  investigated. 

The  results  for  equilibrium  at  the  gas  temperature  show  very  little 
dependence  on  chese  parameters,  mainly  because  of  the  very  low  electrical 
conductivity,  so  that  the  friction  and  heat  transfer  terms  in  the  flow 
equations  dominate  those  describing  the  MFD  interactions.  Thus,  since  the 
effect  of  power  extraction  is  negligible,  the  solution  is  almost  identical 
with  that  for  zero  magnetic  field,  the  only  effect  being  that  the  length  for 
supersonic  flow  decreases  slightly  with  increasing  magnetic  field,  because  of 
the  contribution  of  the  magnetic  drag  term  in  the  momentum  equation. 

In  contrast,  the  results  for  elevated  electron  t emp era ture  show  considerable 
dependence  on  the  magnetic  field  (see  Table  25*5):  the  electrical  conductivity  is 


25*S 


high  C^l  0*  Hi-fi--  ./a)  so  that  electrical  power  ext  a  tion  and  magnetic  drag  control 
the  flow.  The  latter  is  particularly  important .  the  length  for  supersonic  flow 
being  significantly  altered  for  varying  magnate,  field.  ^  Pig«  25*9)  *  Pig*  25*10 

shows  the  variation  of  Mach  number  with  distance  for  various  alues  of  the  load 
factor,  K.  A  reduction  in  K  causes  a  significant  decrease  in  the  length  for 
supersonic  flow,  in  accordance  with  the  expression  for  magnetic  drag. 

Pigs*  25*4-  -  25*10  illustrate  variation  of  the  flow  parameters  with  distance 
and  magnetic  field.  Except  in  the  case  of  F'g.  25**0  the  load  factor  has  a 
constant  value  of  0*8.  On  thess  graphs  the  eu*  for  B  -  0  represents  the 
cases  both  of  zero  magnetic  field  and  equilibrium  at  the  gas  temperature,  since, 
as  stated  above,  the  results  agree  to  within  a  f ew  percent* 

Pigs.  25.11  ~  25*13  demonstrate  the  variation  with  distance  of  the 
electrical  parameters  for  equilibrium  at  the  gas  temperature.  In  Pig.  25*13  a 
value  of  1.0  was  used  for  B;  however*  since  in.  this  case  the  conductivity  and 
velocity  are  virtually  independent  of  magnetic  field,  the  variation  of  power 
output  with  distance  may  be  obtained  directly  for  any  &; ue  of  B  by  scaling 
according  to 

P  =  o  V2B2K  ( 1  -K)  d 

The  results  for  elevated  electron  temperature  are  shown  in  Figs.  25*14  - 
25*18  as  a  function  of  distance  and  magnetic  field.  Since  the  length  for 
supersonic  flow  varies  so  widely  in  this  case,  the  graphs  are  plotted  against  a 
reduced  distance  which  is  the  ratio  of  the  distance  from  ini  at  to  channel  (x) 
to  the  length  of  supersonic  flow  (x_)  (i.e,  di*tan -e  at  which  M  =  '*}. 

25-2*5  Discussion 

25*2*5*1  Supersonic-subsonic  transition 

The  system,  expansion  nozzle  -  MPD  channel  ~  diffuser,  for  the  DID  generator 

A  ; 

is  shown  in  Fig.  25.19(a)  and  the  originally  predicted  performance  (namely, 
acceleration  from  subsonic  to  supersonic  in  the  expansion  nozzle,  supersonic 
flow  in  the  channel  with  is  1  at  the  exit  and  subsequent  subsonic  diffusion  of 
the  flow)  of  curve  1  in  Fig*  25.19(b).  The  present  calculations  indicate  that 
when  B  =  0  or  when  the  gas  temperature  and  electron  temperature  are  equal  and 
the  effect  of  power  extraction  on  the  flew  is  small.  H  -  1  will  occur  at  about 
14  cm  from  the  beginning  of  the  channel  (see  Table  25*3)=  A  lowing  for  the 
probable  inaccuracy  in  the  estimation  of  the  fricti an  and  hsat  transfer 
coefficients,  the  present  1 5  cm  channel  is  satisfa  tory  for  these  conditions. 


The  case  where  T 


considerably  different  results. 


>  T  however,  gives 


Fig.  *5*4=  indicates  tnat  H  =  1  will  occur  well  before  the  end  of  the  channel | 
a  typical  result  of  this  type  is  shown  by  curve  2  of  Fig.  25.19(b)*  M  =  1  being 
prediuwed  a1-  the  point  A.  This  is  physically  impossible  in  a  constant  area 
channel  with  boundary  layers  of  constant  thickness.  It  is  probable  under  these 
u^nditions  that  the  boundary  layer  will  become  progressively  thicker  as  the  Mack 
number  decreases  to  unity  and  thinner  when  the  hulk  of  the  flow  becomes  subsonic. 
The  flow  in  tnis  region  may  be  very  complex,  with  many  shocks  being  formed,  and 
f  unstaoig.  Choking  in  the  channel  couid  cause  the  mass  fuow  to  be  reduced, 
so  that  supersonic  flow  in  the  nozzle  would  no  longer  be  achieved;  an  unstable 
system  giving  alternate  choking  in  the  nozzle  ana  channel  could  possibly  occur. 

If  this  situation  is  to  be  avoided.  Fig.  25.9  (which  shows  the  distance  for 
which  M  =  1  as  a  function  of  magnetic  field)  may  be  interpreted  as  indicating 
the  maximum  permissible  channel  length  as  a  function  of  magnetic  field,  neglectir g 

the  effects  of  boundary  layers. 

25.2.5.2  Flow  regime 

The  type  of  flow- laminar  to  turbulent  -  occurring  ever  a  flat  plate  or 
in  a  tube  is  determined  mainly  by  the  Reynolds  number,  which  is  a  measure  of  the 
ratio  of  the  inertial  force  to  the  viscous  force.  For  a  tube  this  is 

V  D  VI 

Re  =  p  ,  JIL 

*1  V1 

where  Vm  is  the  mean  velocity  across  the  tube,  D  is  the  hydraulic  mean  diameter, 
is  the  dynamic  viscosity  and  v ,  s  ^Vp  the  kinematic  viscosity.  The  flow  is 
usually  laminar  for  Re  <  2500  and  turbulent  for  larger  values.  However,  this 
condition  is  not  strictly  obeyed:  for  example,  when  extreme  care  has  been  taken 
in  the  design  and  construction  of  experimental  systems  laminar  flow  has  been 
observed  at  far  higher  Reynolds  numbers.  A  transition  region,  in  which 
unsteady  turbulent  flow  occurs,  usually  exists  up  to  Re  *  1o\  above  which,  for 
many  cases,  fully  developed  turbulent  flow  may  be  assumed.  Fully  developed 
turbulent  flow  and  fully  developed  laminar  flow  are  amenable  to  a  mathematical 
treatment  but  there  13  no  adequate  model  for  the  transition  region. 

Using  the  continuity  equation  (25.5)  with  the  mean  velocity  V 
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.....  (25.25) 


and,  taking  p,  to  be  that  for  pure  helium  and  B  to  be  the  hydraulic  mean 
diameter  of  the  channel, 


D  * 


4A 


. « . » •  (25*26) 
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The  Reynolds  number  is  shewn  in  Fig,  25,20  as  a  function  of  temperature.  Fully 

developed  turbulent  flow  may  not  be  set  up,  since  the  Reynolds  number  is 

I 

considerably  less  than  10^. 


To  estimate  the  velocity  profile  and  thickness  of  the-  laminar  boundary 

sublayer  in  the  channel,  fully  developed  turbulent  flow  will  be  assumed  and  the 

15 

relations  given  by  Eckert  "  used.  The  velocity  profile  may  be  divided  into 

three  regions  depending  on  the  distance  away  from  the  tube  wall.  Using  a 

+  + 
dimensionless  velocity  V  and  a  dimensionless  distance  from  the  tube  wail  y 

these  regions  are: 
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In  equations  (25*27) 

yt  =  JjT 

N  w 


where  7^  is  the  wail  shearing  stress  given  by 
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so  that 


y!  =  0.1  Re8 


V'  ^  6.2  Re 
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8  I_ 
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.....  (25.29) 
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where  /2  is  the  hydraulic  mean  radius  ai^d  y  the  distance  from  the  wall,  and 
is  the  axial  velocity,  given  by  Vo  =  -2 .  This  velocity  distribution, 

V  *  Oil 


/V  ,  is  shown  in  Fig,  25*21  as  a  function  cf  VD  for  a  Reynolds  number  of  5000* 

In  discussing  turbulent  flow  and  velocity  profiles  it  has  been  assumed 
that  the  cross-section  of  the  channel  is  circular,  taking  a  mean  diameter  from 
equation  (25. 2o).  Most  MPD  channels,  for  magnet  economy  and  other  reasons, 
have  a  rectangular  cross-section  and  the  boundary  layer  thickness  will  not  be 
uniform  over  the  channel  perimeter.  Curves  of  constant  velocity  for  a  pipe 
of  rectangular  cross-section  (as  measured  by  Nikuradse ,  in  Schliehting  J)  are 
shown  in  Fig.  25.22  as  an  example  of  this  effect. 
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the  channel. 


An  additional  complication  concerns 


en  trail ie  onditions  to 


If  a  fluid  enters  a  channel  from  a  large  container,  the  vel 
in  the  inlet  region  varies  with  distance  along  the  channel. 


city  distribution 
hear  the  inlet 


region  the  velocity  distribution  is  almost  uniform, 
effects  until  a  fully-developed  velocity  profile  is 
constant  further  downstream.  It  is  this  constant 


changing  through  viscous 
established  and  remains 
elocity  profile  which  has 


been  discussed  above  and  which  is  shown  in  Pig.  According  to  Nikursdse 

the  fully-developed  turbulent  velocity  profile  for  a  pips  of  circular  cross- 
section  exists  after  an  inlet  length  of  25  tc  40  diameters  from  entrance.  The 
value  for  a  rectangular  channel  is  not  known,  even  less  the  velocity  profile 
for  a  supersonic  flow  emerging  from  a  convergent-divergent  nozzle.  Thus  the 
above  assumption  of  a  fully  developed  velocity  profile  may  well  be  incorrect, 
but  is  limited  by  present  knowledge  in  fluid  mechanics. 


The  above  discussion  concerning  the  velocity  profile  in  a  channel  refers 
vO  conditions  in  the  absence  of  a  magnetic  f ' exd»  For  flow  in  a  magnetic  field, 
insufficient  data  is  available  at-  present  to  permit  definite  predictions  of 
turbulen*  or  laminar  flow  for  a  given  configuration.  A  magnetic  field  acting 
on  an  electrically-conducting  fluid  tends  to  suppress  turbulence  through  the 
internally  induced  electric  fields  in  small  regions  of  turbulence,  the  overall 
effect  being  to  flatten  the  velocity  profile.  Thus  in  the  channel  of  the  IRD 
device  the  flow  may  be  in  the  turbulent  transition  region  when  the  magnetic 
P— aid  is  zero  but  laminar  for  non—sero  magnetic  fields  above  a  certain  value* 

Look*  found  that  a  transverse  magnetic  field  has  a  stabilising  effect  on 
the  laminar  flow  of  a  conducting  fluid  between  parallel  plates.  The  experiments 
of  Lehnert  ,  however,  showed  a  destabilizing  effect  due  to  the  magnetic  field. 
This  illustrates  that  at  present  the  interaction  of  a  magnetic  field  with  a 
moving  electrically-conducting  fluid  is  not  completely  understood, 

IQ 

Wu  has  considered  Hartmann  flow  (laminar  and  incompressible  with  viscous 
and  magnetic  forces)  to  investigate  the  effects  of  reverse  currents  in  the 
boundary  layers  on  the  power  output  of  an  MPD  generator*  The  reverse  currents 
arise  because  the  velocity  in  the  boundary  layer  may  be  xess  than  the  ratio  of 
the  electric  to  magnetic  field  (-  /b)  so  that  the  plasma  is  accelerated  in  this 
region  rather  than  being  slowed  down.  Thus  part  of  the  power  generated  in  the 

main  stream  is  lost  in  the  boundary  layer.  Wu  shows  that,  for  the  same  pressure 
gradient: 
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where  P  *  F(y)  is  a  factor  determining  the  velocity  distribution  across  the 
channel.  For  Hartaann  flow 

<F> s  1  -  .....  (25.31) 

« 

for  large  values  of  the  Hartmann  number*  For  the  IRD  generator f  taking 

typical  values  in  the  channel,  60  ari  K  =  0.5.  the  power  ratio  i^  —  9^. 

>20  **  -  .  —  - 

Blackman  and  Jones^  estimate  that  the  power  ratio  may  be  as  Ivw  as  2%  for 

small  shook  tube  channels. 

25.2.5.3  Channel  wall  losses 

< 

The  friction  factor  at  the  wall  of  the  alumina  channel  has  ‘been  estimated 

assuming  helium  only  to  be  present.  Cesium  affects  the  viscosity  appreciably 

2i 

owing  to  the  large  mass,  as  has  been  shown  theoretically*”.  No  experimental 
data  are  at  present  available  and  no  quantitative  estimate  of  the  effect  of 
seeding  has  been  attempted.  The  tendency  will  be  for  the  friction --factor  to 
increase  with  increasing  cesium  seeding. 

Heat  loss  through  the  walls  has  been  estimated  for  alumina,  assuming  only 
convection  and  conduction  in  the  gas.  There  may,,  in  addition,  be  appreciable 
radiation  from  a  hot  cesium  plasma,  which  should  be  taken  into  account  ir.  the 
energy  balance:  no  experimental  data  az-e  available  upon  which  quantitative 
estimates  may  be  based. 

Another  assumption  implicit  in  neglecting  wail  properties  is  the  neglect 
of  the  power  _ost  by  current  leakage  through  the  alumina.  Consideration  of 

^  ~~  y  .  ~  1  1 

the  highest  measured  electrical  conductivity  for  alumina  {a  =  10  mho/ cm 
at  1b00°K)  shows  that  the  loss  of  power  is  small  ly-  watts)  compared  with  other 
losses  (see  below  for  end  leakage  and  eddy  current  losses), 

25.2.5.4  Electrode  processes 

23 

Pain  and  Smy’s  power  generation  experiments  with  a  combu  tion-driven 

argon  shock  tube  showed  that  large  power  outputs  could  be  obtained  through  an 

MFD  interaction  without  any  apparent  evidence  of  emission  from  the  electrodes. 

24 

These  results  are  contrary  to  those  of  Maycook  et  al  ,  who  have  observed 
considerable  deterioration  of  electrodes  in  an  open-cycle  combustion  device. 

There  appears  to  be  no  conclusive  evidence  on  the  electrode  emission  mechanism 
in  MPD  generators. 

25.2*5*5  End  loop  and  leakage  losses 

At  the  entrance  to  aniexit  from  the  magnetic  field  region  there  is  a 
reverse  flow  of  electric  current,  known  as  eddy  current  or  end  loop  loss. 


Tills  has  bsen  oonsidsred  in  59*e  detail  by  Sutton  ^  and  Dzang^  amongst  others* 
Sutton  found  that  the  losses  increase  as  the  loading  factor  increases  or  as 
the  aspect  ratio  of  the  channel  decreases.  The  losses  can  best  be  offset  by 
extending  the  magnetic  field  region.  Szung  considered  the  Hall  effect  in  his 
evaluation  and  found  that  the  losses  increase  rapidly  with  for  a  continuous 
electrode  system,  but  are  not  so  severe  for  a  segmented  configuration.  Numerical 
values  showing  the  loss  in  pow*r  output  may  be  found  in  these  papers. 


In  addition  to  these  losses,  the  power  output  of  a  generator  may  fee 

considerably  reduced  by  leakage  currents  from  the  generator  region  to  an  earthed 

on 

portion  of  the  system.  Way  et  al**  have  estimated  this  effect  for  constant 
conductivity  throughout  the  considered  length.  Their  theory  has  been  extended 
here  for  non-constant  conductivity*  Under  these  conditions  the  total  power 
output  cf  a  segmented  electrode  generator  when  end  leakage  cc "nrs  it 

P  ^  c  vVa  (1  -  K)  (i  -  C)  LA  .....  (25.32) 

where  c  and  7  are  the  average  conductivity  and  velocity  in  the  generating 
region  and  G  is  the  end  leakage  factor  given  by  (see  Appendix  25.6); 
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For  a  given  §  in  the  generator  region  this  expression  shows  that  end  leakage 
decreases  when  /2s  decreases  and  when  S/V  decreases*  For  the  present  HD 
system  /2s  is  large  (-*  vl)  and  z/e  in  the  downstream  part  of  the  generator 
vrill  be  close  to  unity,  swing  to  the  small  amount  cf  ze combination  occurring 
over  the  distance  s(~  i  cm.)j  thus  the  end  losses  may  be  expected  to  be  high. 

Fig.  25.23  ohows  C  as  a  function  of  V23  for  several  values  of  /?  and  when  ffs/o»1 , 
while  Fig,  25.24  shows  C  as  a  function  of  "s/c  for  several  values  of 

To  reduce  end  leakage  there  must  be  a  considerable  length  of  free  insulator 
at  each  side  of  the  active  generator  region  and  the  conductivity  should  be  small 
outside  the  active  region.  Xf  the  magnetic  field  is  correctly  designed  there 
should  be  little  leakage  upstream  when  L/2s  is  small.  Downstream  of  the 


generator  region  the  same  applies,  with  the  addition  that  recombination  will 
also  act  to  decrease  the  electron  density,  which  however  is  high  in  the  generator 
i‘*gion  (see  Fig.  25.25(c)).  If  the  downstream  leakage  is  appreciable  and  the 
length  of  insulator  is  limited  by  external  pressure  requirements  it  may  be 
necessary  to  introduce  some  electron  attaching  eie.nent  _r 
the  generator  region  to  reduce  the  conductivity. 


radical  dow. .stream  of 


25.2.5.6  Cesium  dissociation 


The  effect  fji  the  dissociation 


Cesium  tv  the  reaction 


Cs  #  Cs  +  e”  .....  (25*34) 

on  the  t pacific  heat  of  an  equilibrium  mixture  of  oesium  is  shown  in  Fie.  25.2b 
(from  Ciifton  }•  The  ranges  of  temperature  ana  pressure  in  the  channel  for 
the  conditions  of  interest  here  are  900°  -  ISOCTK  and  40-120  -mm  of  Hg  (total 
pressure),  that  is.  cesium  pressure  0.4  ”1*2  mm  of  Hg.  For  these  conditions 
the  effect  of  dissociation  is  small  (Fig*  25.26)  and  the  snecifie  heat  of  cesium 
is  close  to  the  value  m/2  ~  4*90  tal/vK  mole  for  any  monatomic  gas,  including 
helium. 

Fig*  25*27  (also  from  Clifton  ')  shows  the  percentage  of  Cs:J  in  cesium 
vapour  es  a  function  of  temperature  for  various  total  cesium  pressures.  The 
lowest  cesium  temperature  and  pressure  anticipated  in  the  chcnnel  will  be  $00°K 
and  0*4  mm  Hg  (at  the  channel  entrance).  For  these  conditions  less  than 
0,01  atomic  percent  Cs2  will  be  present  (Fig,  25.2?) j  thus  there  will  be  no 
noticeable  contribution  to  the  specific  heat  of  cesium  from  the  reaction 

2Cs£  Cs2  .....  (25*35) 

Figs*  25*26  and  25.27  have  been  derived  assuming  thermodynamic  equilibrium.  For 
a  non-aqui librium  system,  as  for  example  when  electron  heating  occurs, 
dissociation  of  cesium  by  the  above  reactions  (equations  25,34  -ana  25.35)  may 
have  a  pronounced  effect  on  the  specific  heat  cf  the  gas  mixture* 

25.3  IONIZATION  KOK -EQUILIBRIUM 

25*3.1  Introduction 

In  the  previous  section  the  results  obtained  when  the  MFD  equations  for 
constant  area  flow,  taking  into  account  wall  friction  and  heat  transfer,  were 
solved  for  the  eases  of  ionization  at  the  temperature  and  at  the  electron 
temperature.  The  computer  programme  used  was,  however,  not  sufficiently  fast 
or  versatile  to  permit  all  the  required  investigations  to  be  undertaken. 
Preliminary  results  obtained  with  a  faster  (machine  orders)  rrograsae,  in  which 
several  important  modifications  and  generalisations  of  the  previous  equations 
have  been  made,  will  now  be  discussed.  One  of  tne  most  important  modifications 
is  the  use  o 1  an  equation  governing  magretically-induced  ionisation  non- 
equilibrium. 


25.3.2  Equations 

The  energy  equation  for  s  compressible  fluid  in  -ha  absence  of  magnetic 


(25.36) 

(25.37) 


The  heat  transfer  coefficient  (h)  in  equation  (25*36)  is  obtained  from  the 
friction  factor  (f )  using  Reynolds  analogy: 

b  *  pVCptf/2,  .....  (25.38) 

iihich  is  accurate  to  a  few  percent  for  fully  developed  turbulent  flow*  Using 
the  friction  factor  (see  below) 


fifes 


f  =  0.079/S^l/4 
h  =  O.C395PVC 


(25.39) 

(25*40) 


An  alternative  way  of  expressing  the  heat  transfer  coefficient  is  by  means  of 
the  Stanton  number: 

h  =  pVC  St 
P 

Using  this  expression,  together  with  the  relation 

St  =  Nn/PrRe, 

tl»  empirical  relation  obtained  by  Durham  et  al^  for  heliiss  (neglecting  entrance 
effects)  is 

Hu  =  Q.034Re^*^  Pr^,i*‘ 
and  a  Prandtl  masber  of  0.75  gives 

h  =  0*040 5pVC j/He0,2  , 

which  is  in  fair  agreement  with  equation  (25.40)*  Equation  (25*40)  is  used 
in  tlie  present  calculations. 


The  momentum  equation  for  a  compressible  fluid  in  the  absence  of  magnetic 
interactions  is 

S££2jM£  . (25.41) 

where  the  friction  coefficient  (?)  is  defined  as  the  ratio  wall  shearing  stress 
to  the  ^rnamis  head  of  the  stream: 

f  =  Zr/pV2  .....  (2  .42} 

(Note:  the  friction  coefficient  is  sometimes  defined  as  four  times  greater  than 


25*16 


.  ^  J 

this)*  According  to  Eckert  and  Drake  for  fully  developed  turbulent  flew  the 
friction  coefficient  is  given  by  Blasius's  las 

1/4  .....  (25.43) 


where  a*  is  a  coefficient  equal  to  0,079  and  the  Reynolds  number  is  calculated 
with  the  mean,  velocity  in  the  channel  and  the  hydraulic  diameter  (D) . 

Combination  of  equations  (25.3),  (25.36)  and  (25*40 )  with  the  addition  of  the 
electrical  energy  extraction  term  yields  the  MPD  constant-area  energy  equation 


dx 
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• . *  * .  (23»44) 


Combination  of  equitisns  (25.3),  (25.41),  (25.43)  and  (25»A,l)»  together  with  the 
magnetic  braking  force  tens,  yields  the  MPD  constant  area  momentum  equation 
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For  a  coordinate  system  with  a  plasma  flowing  in  the  positive  x  direction  and  a 
magnetic  field  in  the  positive  2  direction 
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and 
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where  for  a  solid-electrode  generator  =  0,  for  a  sepsented-eleetrode  generator 
Kg  =  1 -K,  and  for  a  Hall  generator  K  -  0%  in  the  present  calculations  a  segmented- 
electrode  system  has  been  evaluated. 

One  of  the  ©ore  important  differences  between  these  investigations  and 
those  reported  in  Section  25.2  is  the  use  of  an  equation  giving  the  rate  of  change 
of  electron  temperature  due  to  electron  heating,  coll*sional  energy  exchange  and 
ionization  or  -deionization  processes 
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In  equation  (25*48). 
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using  a  three-body  electron-eleotron-ion  recombination  coefficient  a  and 
asaming  that  ionisation  is  determined  from  this  by  the  condition  that  the  two 
rates  must  be  equal  at  equilibrium, 
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The  total  rate  of  change  of  electron  concentration  is  obtained  from  the  rate  of 
change  due  to  ionisation  and  recombination  plus  the  rate  of  change  due  to  overall 
changes  in  the  gas  density!  thus 
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From  the  enthalpy  equation  (25*5) 
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and  from  the  continuity  equation  (for  a  constant-area  generator) 
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the  six  simultaneous  differential  equations  (25.44),  (25.45),  (25 .48),  (25.51) 

(25*52)  and  (25*53)  are  integrated  using  the  auxiliary  equations  (25*46),  (25*47) 
(25*49)  and  (25.50). 


The  Saha  fractional  ionisation  is  found  from 
2  5/2 
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where 
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*nd  d  ■=  eV^/k 

The  equilibrium  electron  and  atom  concentrations  are  given  by 
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fhe  non-equilibrium  at®  concentrations  are  obtained  from 
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fhe  scalar  electrical  conductivity  ^.s  found  from  the  usual  expression  (see 
Section  25.2*2)  and  the  Hall  and  ion  effects  are  found  using 
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fhe  remaining  quantities  required  for  the  solution  of  the  six  differential 
equations  are  given  in  Section  25*2. 


25.3*5  Numerical  solution  of  equations 

The  set  of  equations  given  in  Section  25*3*2  say  be  solved  by  a  numerical 
step  method  to  give  the  variation  of  flow  and  field  parameters  along  the  channel 
length.  Initially  the  values  of  the  sain  variables  (T  ,¥,T,T  ,n  and  p)  are 

O  0  6 

kuewn  at  the  channel  inlet.  From  equation  (25*59),  knowing  the  mass  flow,  the 
channel  cross-sectional  area  and  the  gas  constant  the  total  pressure  (p^)  may  be 
found.  Knowing  and  the  seeding  fraction  (defined  in  (25.60))  the  cesium  atom 
concentration  (n  )  can  be  obtained  from  equation  (25*57)  and,  using  (25*56),  the 
belli®  atom  concentration.  Knowing  the  electron,  ce  c£  tu£  and  helium  atom 
concentrations  the  total  electronic  and  ionic  collision  frequencies  can  be  found 
using  equations  (25.14)  and  (25.64);  all  collision  cross  sections  are  assiffled 
constant  in  these  preliminary  calculations  the  exception  of  that  for  electron 

ion  collisions  (equation  (25.1?)).  Using  equations  (2 5.62),  (25.63)  and  (25.11 ), 
Pmt  ai‘id  a  can  now  be  evaluated.  Knowing  the  seeding  fraction  and  the  electron 
and  cesium  a  tom  concentrations  the  non-equilibrium  fractional  ionisation  (X)  is 
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found  fret  equation  (25*61 )*  Then,  using  the  equilibrium  fractional  ionisation 
tmm  equation  (25.54),  the  equilibri®  election  and  cesium  atom  concentrations 
*r*  found  froa  equations  (25. 55)  and  (25.56)* 

This  rather  inelegant  procedure  ^outlined  above)  is  required  if  the 
asa^ption  of  low  fractional  ionization  is  not  valid  (as  is  often  true  for  the 
case  of  extrathersal  ionization)  and  when  non-equilibrium  extrathermal  ionization 
is  investigated,  that  is  when  equation  (25*46)  is  used  rather  than  the  simple 
expression  derived  froa  it  (equation  (25.22))* 

It  is  row  possible  to  integrate  the  six  simultaneous  differential 
equations  (25*44),  (25*45),  (25.48),  (25*51 ),  (25*52)  and  (25.53).  This  is 
carried  out  using  the  fourth-order  Runge-ICutta  process  which  r  equires  evaluation 

of  each  of  the  derivatives  at  the  point  under  consideration j  these  are  obtained 
as  follows. 


From  equation  (25*46)  B.ji  is  calculated  using  the  values  of  0  t  and  c 
found  previously,  then  using 
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the  Reynolds  maber 

16  "3^  .....(25*66) 
*ay  he  obtained.  Using  a  recoveiy  factor  of  0,85  in  equation  (25.37),  T  is 
obtained  and,  since  Tw  is  supplied  as  data,  dT^/dx  (equation  (25,44))  is  obtained. 
Uaing  dT^/dx  aid  evaluating  equation  (25-47)  dff/dx  is  found  from  equation  (25*45). 
Staining  j. J*  from  elation  (25*49)  and  the  values  of  (n  }  .  (n  )  and  v 
calculated  above.  dT^/dx  is  found  free  equation  (25*48).  Slewing  dT ^dx  and 
dS/dx,  equation  (25*52)  gives  dT/dx;  equation  (25*53)  gives  d#/dx  and  thus 
dr.y'dx  can  be  evaluated  using  equation  (25.51). 

The  integration  process  gives  values  of  T_,  V,  T,  T  ,  p  and  n  at  x  -  h 
Aere  h  is  the  step  length.  The  above  procedure  is  then  repeated  along  the 
channel  until  the  exit  distance  is  reached  or  until  the  Mach  number  becomes  unity, 
whereupon  the  calculation  stops. 


25*3*4  Discussion  and  results 


Several  important  differences  between  the  investigations  reported  in 
Sections  25*2  and  25.3  have  already  been  men' loned.  In  addition  to  these 
differences  the  gap  constant,  specific  heat  a*  constant  pressure  and  ratio  of 
•pacific  heats  used  in  the  calculations  of  section  25*3  are  those  for  a  helium- 
essiua  mixture  rather  than  for  pure  helium  as  in  Section  25.2*  Thus,  although 


25*20 


the  channel  entrance  data  are  stlf -consistent  and  closely  approximate  the  nozzle 
exit  conditions  in  most  respects,  there  if  not  complete  agreement  since  the 
present  negate  calculations  employ  data  for  pure  helium*  The  manner  by  which 
the  gas  constant  for  a  cesium-helium  mixture  is  calculated  is  given  in  Chapter  23 
of  this  report*  for  a  1a/e  cesium-helium  mixture  the  gas  constant  is  about 

_  I 

75  per  sent  of  that  for  pure  helium. 

The  essential  difference  between  the  present  calculations  and  the  earlier 
ones  (Section  25*2)  is  that  non-equilibrium  values  of  electron  temperature  and 
concentration  may  be  evaluated,  Thus  the  effects  on  the  flow  and  electrical 
parameters  of  entrance  to  and  exit  from  the  generator  region  may  he  examined 
in  some  detail.  For  example,  if  non -equilibrium  flow  occurs  in  the  nozzle,  the 
effects  on  generated  power  in  the  chancel  may  be  examined?  the  decay  of  electrical 
conductivity  on  exit  from  the  magnetic  field  region  and  its  effect  on  end  leakage 
can  also  he  investigated*  The  continuation  of  non-equilibrium  flow  from  the 
nozzle  into  the  channel  in  the  absence  of  magnetic  field  can  also  be  examined. 

In  the  present  programme  the  magnetic  field  is  a  step  function*  being  zero  outside 
the  generator  region*  Later  calculations  may  enable  the  magnetic  field  strength 
to  be  varied  as  a  function  of  channel  distance*  (In  the  present  programme  only 

the  wall  temperature  (T  )  can  be  expressed  as  a  function  of  channel  dis  Lance i 

w 

when  sufficient  experimental  data  on  wall  temperature  are  available  this  will 
be  used*) 


Under  typical  operating  conditions  in  the  channel  the  heat  transfer 
coefficient,  h  ^ equation  (25*40) ,  derived  from  Reynolds  analogy),  used  in  Section 
25*3  is  about  three  times  greater  than  that  of  Section  25*2;  this  will  not  make 
a  great  difference  to  the  flow  conditions  when  power  extraction  occurs. 


The  only  calculations  carried  out  to  a  ate  have  been  with  a  segmented-electrode 
system;  however,  the  programme  can  readily  be  used  for  either  a  solid-electrode 
or  Hall  generator  configuration,  but  no  account  is  taken  of  finite  segmentation* 

In  later  versions  of  the  HD  MPD  channel,  electrodes  are  placed  only  in  the  central 
region  of  the  channel,  and  long  blank  sections  are  left  at  either  end  to  reduce 
end  leakage.  The  flow  in  these  three  separate  regions  of  the  channel  can  be 
examined  by  varying  the  loading  faster  appropriately,  that  is,  Kw!(and  therefore 
K^=0)  in  the  channel  region  where  no  electrodes  exist* 


In  equation  (25*48),  §  is  a  factor  account 
non-Maxweilian  distribution,  given  by 
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where 
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In  the  present  calculations  1*1 


j 

is  assumed* 


but  thii 


take  into  account  such  effects  as  inelastic  collisions  due 


readily  be  changed  to 
to  molecular  impurities. 


In  equation  (25*50)  the  three  body 
coefficient  used  is  that  given  by  Hinnov 

a  =  1.1*10  ~T  ; 
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this  has  been  shown  to  be  the  most  appropriate  under  MPD  generator  conditions^, 
being  close  to  the  recent  value  obtained  by  Bates  at  al^&. 


The  data  used  for  the  calculations  reported  here  is  tabulated  in  Table  25*4 
it  is  similar  in  several  respects  to  that  given  in  Tables  25*2  and  25.3  for  the 
calculations  of  Section  25*2.  Pigs.  25.28  to  2 5.35  inclusive  show  the  results 
obtained  using  this  data  for  a  magnetic  field  strength  of  0.3T,  (Note  the 
semi-logarithmic  scales  of  Figs.  25*32,  23.33  and  25*34.)  These  curves  may 
be  compared  with  the  corresponding  ones  for  the  calculations  of  Section  25*2. 
While  the  general  form  of  most  the  the  curves  correspond,  there  are  several  points 
worthy  of  note,  as  outlined  below. 


Fig.  25-32  shews  the  electron  temperature  as  a  function  of  distance  down 
the  channel.  The  rapid  rise  of  electron  temperature  (within  the  first  centimetre 
from  channel  inlet)  is  clearly  visible.  This  finite  rise  time  for  the  electron 
temperature  is  governed,  through  equation  (25 .48), by  the  rate  at  which  energy  is 
supplied  to  the  electrons  by  the  electric  field,  the  rate  at  which  it  is  lost  by 
elastic  collisions  and  by  ionising  collisions.  When  these  processes  equilibria te 
the  electron  temperature  is  close  to  the  value  given  in  Fig.  25.14  which  was 
obtained  using  the  equilibrium  form  of  equation  (25.48)  (that  is.  equ-tion  (25,22)) 
for  the  same  magnetic  field  (the  minor  differences  are  due  to  the  use  of  slightly 
different  constants). 

However,  although  the  electron  temperature  rises  rapidly,  the  electron 
concentration,  which  is  obtained  from  equation  (25,51),  follows  this  variation 
much  more  slowly,  as  may  be  seen  from  Fig,  25*33  (compare  with  Fig,  25.16), 

Taking  an  average  flow  velocity  of  3000  s/sec  in  the  channel  the  rise  time  for 
electron  concentration  is  *.30  psee,  compared  with  3  uses  for  the  electron 
temperature.  This  rate  of  rise  of  electron  concentration  is  determined  from  the 
reverse  process  of  three  body  electron-eiectrou-ion  recombination  as  given  by 
Hinnov  and  Hirsc hberg*  While  this  procedure  may  be  in  some  error  far  from 

equilibrium,  an  unpublished  survey  of  several  quoted  eleutron-atom  ionization  rates 


25.22 


snows  tnai  they  are  considerably  further  from  agreement  than 

rates  _ 


recombination 


The  alow  rate  of  rise  of  electron  concentration  will  probably  not  be 
significant  in  a  large-scale  MPD  generator;  but  the  influence  on  the  HD  channel 
nay  be  very  important* 

Fig.  2p,J4  shows  the  scalar  electrical  conductivity  of  the  helius«*oesium 
pusssis  determined  from  equation  (25.11 J*  As  expected e the  conductivity  follows 
the  general  form  of  the  electron  concentration  curve,  being  initially  low,  and 
rising  slowly  to  a  value  of  *^130ano/s  (compare  Fig=  25*17)*  The  small  drop  in 
conductivity  in  the  first  centimetre  channel  distance  is  a  real  effect  and  occurs 
because  the  electron  temperature  here  has  increased  more  rapidly  than  the  electron 
concentration* 

The  specific  power  (power  output  per  unit  channel  volume)  is  shown  in 
Fig*  25*35  as  a  function  of  distance  down  the  channel.  As  with  the  electrical 
conductivity  there  is  a  small  drop  in  power  in  the  first  centimetre  of  channel, 
for  the  same  reason.  The  specific  power  then  increases  rapioly*  being 
determined  primarily  by  the  rate  of  increase  of  electrical  oond.uotiyity.  A 
maximum  value  (--  32  MWe/ )  is  reached,  beyond  which  the  specific  power  decreases, 
-^he  r  sa^son  fur  this  is  that  the  electrical  conductivity  has  become  sensibly 
constant  but  the  flow  velocity  is  decreasing  at  an  increasingly  rapid  rate.  The 
area  under  the  curve  multiplied  by  the  channel  cross  sectional  area  (4.8  cm^) 
yields  an  ideal  power  output  of  about  G.SkWe,  (neglecting  effects  of  finite 

segmentation)  for  a  magnetic  field  (B)  of  0*5  T  the  newer  cutout  increases 

2  '  * 
approximately  as  B  * 

Comparison  with  the  previously  determined  corresponding  figures  for  B=Q«5 
shows  that  a  Mach  number  of  unity  for  the  present  calculations  occurs  at  x*14cm 
from  inlet,  compared  with  about  10  cm  for  calculations  of  Section  25.2  (see 
Fig.  25*9  for  example).  This  arises  mainly  because  of  the  lower  average  magnetic 
interaction  in  the  latest  calculations,  mainly  due  to  lower  levels  of  electrical 
conductivity. 


While  only  a  few  simple  calculations  using  the  latest  programme  have  been 
reported  here  several  other  investigations  can  and  will  be  made*  Although 
the  programme  is  by  no  means  ideal,  omitting,  for  example,  the  effects  of  finite 
electrode  segmentation,  MPD  boundary  layer  formation,  electrode  sheath  formation 
and  flow  inhomogeneities,  it  is  considerably  more  general  and  versatile  than 
programmes  hitherto  used  for  IK?  investigations. 
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APPENDIX  25. A 

DERIVATION  OF  STAGNATION  TEMPERATURE 


Prom  equations  (25.3)  and  (25.4/ 


mR  dT 
dP  =  -j-v" 

while  from  equation  (25.5) 


mRTdV 

9 

A\r 


VdV 

dT  =  dT  - 
o  C 

P 


Combining  these  equations  gives  an  expression  for  dp  in  terms  of  T  and  V  only: 

o 

dp  =  _  m(x-l)dV  .  +  m(y-1  )dV 

AV  Ay  ^y2  2AX 


. (25.A.1) 


Using  equations  (2?,2)  and  (25*3)  to  eliminate  p: 

mdV  ,  2mfVdx  2/  \ 

+  dp  +  +  oB*(l  -K)Vdx  =  0 

aiid  using  equation  (25.A.1 )  to  eliminate  dp  gives  an  expression  in  which  the  only 
variables  are  x,  V  and  T: 


|  si£±; ' . 

L2Alf 


mRT 


dV  + 


mRdT  [~  - 

AV" 


Vox  =  0 


(25.A.2) 


When  this  is  combined  with  the  energy  equation  (25.1 )  to  eliminate  dr.  a  linear 


differential  equation  of  first  order  for  T  in  terms  of  V  is  obtained; 

o 


*1^,,  ytK/tic,) 

tr/  w  vj2  •.  \  /  „  «*2  \ 


r<K5v2-K6) 


(K5'd-Kg) 


(25.A.3) 


in  which 

^  =  mC/  A 
=  rB2K(l-K) 

5  ■  W* 

=  mS/A 

_  2f  o  J2,  « 

*5  *  XT  -  5B  (’-*) 

Kg  =  *k(T./V  <*-»  )/«pD 

Kj  ~  m(y+l)/2Ay. 
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Equation  (25.A.3)  may  be  written  in  the  usual  form: 


sr  -  vVv)  *  «i<v> 

where  f^(v)  and  g^V)  are  functions  of  V  given  by: 


(25.A.4) 


f,(v)  = 


w2^) 

^VCl^V2- Kg) 


gi(v),T.^~y 

The  solution  to  this  equation  is  obtained  by  multiplying  through  equation 
(25.A.4)  by  the  integrating  factor,  exp  l/^VjdvL  and  then  integrating  the 
resultant  equation.  This  operation  yields  an  expression  far  stagnation  temperature 
in  terms  of  velocity: 


T  = 
o 


where  Xg  =  K^/K^ 

4h(T  -Tw)(^Q 

rCpD  ”aI  "  *B2(l-K]  -  1 


«s  *0 

-S'  )  dV 


(25^.5) 


oB  K(l-K) 


-  oB2(l-K)  (—■  K  -  l)_ 


J_1 

f-1  J 


“d  K1 0=  S'1 


Provided  V2  >  Kg,  (1  -  Kg/V2)  1C  may  be  expanded  binomially  and  the  integral 

in  equation  (2 5. A. 5)  then  becomes: 

f  2K  +2  2K,_  2  K  n-2 

Kj  (V  10  -  .  7  10  +  a  V  10  ...  )  dV 


h\  (V 


,*  2*  2x  -4 

^0  _  10 


-‘iv 


2K  -4 

+  a,V  ...  )  dV 


(25.A.6) 


where 


4a  = 


-  Vr 
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K10(K10"l)if82 

&2  ~ 

21 

K10CK10-1>(K10-',’'K83 

a3  =  3! - 


When  the  magnetic  field  is  non-tero  there  are  no  singularities  so  that  the  integral 
may  be  evaluated  to  give 


T  = 
o 


K  V2*3^  0 
11 


(V2-^) 


S 


V*  + _ _ _ 

2K1q+3  2K1q+1  (2Ki|0-l)V2  (2K1q-3)V4 


'1 


+  ♦ 


const 


+2*10  j 


.....  (25.A.7) 


where 


b1  =  ”  Vi  +  s 

b2  =  ^2a2  “  Mi 
b3  ~  ”  *2*3  *  *3a2 

\  =  h\  ~  *3*3 


etc. 

and  K11=  K/k^. 


From  this  equation  is  obtained  as  a  function  of  velocity. 
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APPENDIX  25. B 


END  LEAKAGE  LOSSES 

The  power  output  of  an  MFD  generator  is 


P  «  1  .  J  LA 

=  E  J  LA 

z  z 

for  the  co-ordinate  system  of  Fig.25*  25(a) 


(25.B.1) 

(25.B.2) 


For  a  segmented  electrode  system  with  no  end  leakage 

Bg  =  “  KVB  .....  (25.B.3) 

Jg  =  c  VB(I-K)  . (25.B.4) 

thus  P  =  o  vVltfl-lOlA  .....  (25.B.5) 

27 

Way  et  al  suggest  that  end  leakage  losses  may  be  accounted  for  by  use  of  an 
end  leakage  factor  which  is  contained  in  the  current  term: 

Jg  =  (l-C)(l-K)orVB  . (25.B.6) 

so  that  P„  *  ovVkO-KKi-CJlA  .....  (25.B.7) 


According  to  their  calculations 


-  gf— 
1+02 


1  - 


2s 


+  1 


(25-B.8) 


L  is  the  active  length  of  generator  and  s  the  distance  between  the  outer 
electrodes  and  the  adjacent  earthed  portion  of  the  (symmetric)  system  (Fig.25*25(a) ) 


Way’s  end  leakage  factor  has  been  estimated  assuming  a  constant  conductivity 
from  x-Q  to  x»L  +  2s.  However,  with  electron  heating,  the  electrical 
conductivity  is  much  higher  in  the  magnetic  field  region  than  it  is  outside  it; 
thus  this  assumption  is  no  longer  time.  For  this  case  an  alternative  expression 
may  be  obtained  for  C,  as  shown  below.  The  magnetic  field  is  assumed  constant 
in  the  region  L  and  zero  outside  these  limits,  corresponding  to  the  ideal 
practical  case. 


For  the  co-ordinate  system  of  Fig.  25*25(a) 

J  s  —2?  (B*  +  0  E*)  . (25.B.9) 

X  1+T  x  z 


3  s  oE* 

y  y 


(25.B.10) 
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r  -  — ^  (r  -  ^e*) 
2  i+r 


.  (25.B.11) 


where 


E*  =  E  +  V  x  B  (see  Harris  and  Cobine^) 


.  (25.B.12) 


V  =  V  ;  V 
x  y 

B  =  B  ;  B 
y  x 


•••»»  (25.B.13) 


so  that 


E*  =  E 

x  x 


(25.B.14) 


E*  =  E 

y  y 


(25.B.15) 


E*  =  E  +  VE 
2  z 


(25.B.16) 


J  =  — ^  (K  +  /SE?  +  /?VB) 
1  +£*-  x  2 


(25.B.17) 


J  =  oE 

y  y 


(25.B.18) 


Jz  *  2  ( E  +  vb  +  pE  ) 

i+p2  z  x 


(25.3.15) 


For  a  segmented  electrode  generator  with  no  end  losses  J  =0,  and  J  =  0;  thus 

x  y 

\  «  -  P(Ez  ♦  VB)  .  (25.B.20) 

When  end  leakage  does  occur  a  similar  form  for  may  be  assumed: 

Ex  =  -  p&{Ez  +  VB)  . (25.B.21) 

where  ’a'  is  a  constant  to  be  evaluated,  then 


(25.B.21 ) 


J  =  -^2  (1-a  )(E_  +  VB) 

X  1 40*  2 


j  x  O  w  (E  +  VB) 
2  2 


(25.B.22) 


(25.3.23) 


Bx  and  Ez  refer  to  the  fields  within  the  generator  region  L.  For  a  segmented 
electrode  geneiator  Ex  is  negative  (see  equation  (25.B.20))  thus  the  potential 
along  the  axis  of  the  generator  is  as  shown  in  Fig,25,25$>)<  The  total  voltage 
ueveioped  axially  along  the  generator  is  -  LEx>  thus  the  field  in  the  end  region 


E  =  - 
s  2s 


and  the  current  density  is 

J  ~  c  E 
x  s  s 


(25*3.24.) 


(25.B.25) 


25.3': 


where  is  the  conductivity  in  the  end  region,  which  differs  from  that  in  the 
active  region  a. 

Using  equations  (25»B.2l),  (25.B.22),  (25«B.2h)  and  (25.B.25)  gives 


cal(H02) 

a  2s 


(25.B.26) 


Prom  equations  (25.B.3),  (25.B.6)  and  (25.B.23) 

C  g  slhzil 


.  (25.B.27) 


(25.B.28) 


2  r 

thus  C  n  ■ ••  1  - — - -  (  . (25.B.28) 

U  )LgS  ''  . 

?sa  T  1 

Inspection  of  this  equation  shows  that  it  reduces  co  Way*s2^  expression  (equation 
(25*B.8))  when  er  =  a .  The  equation  may  be  re-written  as 


-2  L  ^3 

*  a  a 

(1^2)L  au 


(25.B.29) 


s  —  +  2 
a 
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TABLE  25.1 

INPUT  VALUES  TO  CHANNEL  (SEMI -FROZEN  FLOW  IN  NOZZLE ) 


N 

e 

= 

y  * 

xIO18 

m 

"c. 

= 

4 

xIO21 

nh. 

= 

4 

xIO23 

m-3 

'  + 
e-Cs 

x 

1 .2 

Q 

xIO" 

-i 

sec 

e-He 

= 

4.8 

xIO9 

-1 

sec 

e-Cs 

= 

2 .4 

xIO9 

sec  ' 

TOT 

= 

8.3 

xIO9 

sec 

0 

— 

10 

mho.] 

-1 

a 

p  =  5  xlO5  newton  2  (-  0.0.5  ata) 

T  =  874°K 

V  *  4100  3.  sec  ^ 

T  =  2500°K 
o 

TABLE  25*2 

CONSTANTS  OF  CHANNEL 

2  -2  o  -1 

h  “10  watt,  o  .  K 

T  -T  =  1 000°K 
aw  w 

ar  ~  1 .658 

m  =  5.6x10  3  kg.  sec 

f  =  10  fc 

A  =  4.8x1 0"4  m 2 

C  =  5.1 92x1 O3  m2.sec  2.  °K  ^ 

P 

=  8.31  Cxi O3  joule, kgs  .  °K  ^ 

W  =  4.003 


K 


0.8  unless  otherwise  stated 


DISTANCE  FROM  INIS?  AT  WHICH  1  (for  conditions  of  Tables  25.1  and  25.2) 


B  s  0 


X 

s 


14.25  x 


% 


0.5 


1.0 


i;s 


0.5 


1.0 


10.90  ) 

) 

4.043  ) 
1 .62?  ) 
14.18  ) 

14.06  ) 


elevated  T 

e 


equilibrium  T  =  T 


TABLE  25.4 

BATA  FOR  CALCULATIONS  OF  SECTION  2 


helium-cesium  MIXTURES 
CHANNEL  CONSTANTS 

X  =  0.01 

^MIX  =  ^69 
'hix  *  1  *6657 

,  2 

3  =0.5  webers/m 

K  =  0.6 

a’  =  0.079 

CpMIX=  5926  joule/kgm.°K 
a  ~  0.00645  kgw/sec 
D  =  0.01890  a 
A  =  0.00048  m2 

T  =  1250°K 

w 

STARTING  DATA 
Pj  =  0.0376  kgn/m^ 

T1  =  875°K 

n  =  3.1 01^/^ 
e1 

Te1  =  8  ?5°K 
¥1  =  3570  m/' sec 

?  4  =  2500°K 


Kn 


ET 


CONSTANT  AREA  GENERATOR  CHANNEL 


CUCCTRO*  T«MP*RATU«C  (Tfc), 


CkCCTRfl n  TinPCRATUItC 


!  Read  in  constants,  initial  values 
!  of  parameters,  values  of  v 

!  Calculate  constants  independent 
1  of  a 

[  - - - 


Put  *^<r  in  ^  a 


(as  first  approximation) 


‘P 


Calculate  constants  containing 


Calculate 

as  f(V)  [Eqn^5.A.7 

Use  j+1T0,  J+1V 
p,  T,  M  etc. 

to  obtain 

F  Set  e  =  1000,  T, 
o 

*  1 

— y - - 

-  <  — - 

Tem  =  Tem-1  + 
Calculate  tj  (Te) 

*  Z  [Eqn.  25.24 

1 

["Test  if  Z  >  0  “| 

YSS 


NO 


Subtract  rffi  from  Tem 

fm+1  =  cn/10 
Test  if  fffi>10‘4 


YES' 


NO, 


Calculate  J^nen+1 ,  J+\+1 
Test  if  -  ^+1 


n+1 


a 

n 


NO  —  > 


n=n+1 


YES 

t 


f  Calculate  A  x  [Eqn.25,71 
Print  results 


.  Js.J+1 


STOP 
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?b.1  INTRODUCTION 

Direct -current  I4PD  generators  are  generally  of  two  types  according  to  the 
working  fluid:  the  open-cycle  combustion  generator  and 


’■’e  closed-cycle  generator, 
probabiy  with  a  nuciear  heat  source.  High  gas  velocities  are  advantageous  in 
each.  However,  the  losses  associated  with  transonic  flow  will  probably  be 
prohibitive;  gas  velocities  consistent  with  high  subsonic  Mach  numbers  should 
lead  to  greatest  overall  economy. 

The  open-cycle  generator  is  suitable  for  use  as  a  topping  unit  for  an 
otherwise  conventional  steam  cycle  but  it  seems  likely  that  a  certain  amount 
of  oxygen  enrichment  will  prove  to  be  necessary.  Consequently,  the  specific 
heat  (Cp)  of  the  gas  will  be  high  and  the  effective  ratio  jf  specific  heats  ()  ) 
will  not  be  much  greater  than  unity, 1  For  a  given  pressure  ratio  across  the 
generator  the  temperature  drop  will  thus  be  relatively  small  and  constant  velocity 
flow  will  not  differ  greatly  from  flow  at  constant  Mach  number. 

The  closed-cycle  generator  will  almost  certainly  use  an  inert  gas  having  a 
high  ratio  of  specific  heats.  In  this  case  constant  Mach  number  flow  will  almost 
certainly  be  more  economical  than  flow  at  constant  velocity.  It  is  proposed  to 
establish  a  general  treatment  of  the  case  of  constant  Mach  number  flow  of  a 
monatomic  gas,  so  that  the  various  types  of  generator  for  example,  as  described 
by  Lindley”)  become  special  cases. 

26.2  THEORY 

The  one- dimensional  flow  equations  are 


energy:  p  UCp  ~  — 


dx 


_ff  T 

AJ 

X  X 


y  y  hh(  T-To  , 


•  •  .  *  V  £  •  1  } 


,  „  dU  dp 

momentum:  pu  -r~  +  -e*- 
dx  dx 


-1  E  - 


nil 


f  o  C  9  N 

\  •  C  f 


state : 


y 


P  RF 


*  •  *  •  •  ^  * 


Op  A 

iO,  i 


•  •  •  *  * 


continuity t 


V  Uii 


(26.4) 


k  i  number : 


M 


f?6  ■=,) 


Hifferentiating  equation 
(26.tX- 


i  dU 
U  dx 


1_  dT 
2T  dx 


and  substitution  of  ( 


]  gives 
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Ohm’s  law  can  be  written  in  the  form, 
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The  maximum  value  which  E  can  take  is  UB,  therefore  it  will  be  assumed  that, 
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Substituting  equation  (26.14)  into  (26./) 
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A  number  of  possible  cases  may  now  ie  disci 
f  1 )  Segmented  electrodes 


K 


inis  case  K  will  be  specified.  I  will  be  zero  and  thus 
y  ■  x 

A  V 
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(2)  Continuous  electrodes 


"n  this  case  K  will  be  specified  and 
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e  zero. 
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Hall  generator 
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of  series  connection  with  K  a 


in  the  case  (1)  above. 
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Angled  current  generator 
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( 6}  Angled  potential  generator 


mis  case  is  interesting  :  ro~  tne  engineering  point  or'  view  as 


inclined  tubular  duct  walls  having 
the  x  direction  can  be  used,  K  , 
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the  rate  of  energy  dissipation  in  the  gas,  namely*. 
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In  tl.e  general  case  K  .  -  <(1  -  K  ,)  and  thus  K  is  computed  from  (26.13) 

xi  y  i  x 

for  all  other  x.  Furthermore  the  tangent  of  the  tube  angle  (a)  will  always  be 
given  by  £  .  Only  in  case  (6)  above  will  this  angle  be  constant  as  p  will 

undoubtedly  vary  with  x. 


Generally  K  may  be  written  *(1  -  K  )  so  that  by  equation  (2b. 13) 
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In  this  notation  equations  (2b. 15),  (26.16)  and  (26.1/)  reduce  to: 
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The  above  equations  are  now  sufficiently  general  to  treat  cases  (1)  to  (5 
simply  by  choosing  K  and  c ^ .  However,  case  (6)  presents  some  difficulty  in 

that  K  is  a  function  of  K  and  is  thus  a  function  of  x.  In  this  case  equation 
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(26.13)  is  sti1!  true 
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As  a  result,  equation  (26.18)  becomes: 
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■'nils t  in  case  1 6)  (equation  (26.22)  )  becomes : 
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In  the  above  equations  £  =  j~ 

Certain  difficulties  remain, 
ssed  in  dimensionless  form  in  equations  (?6.1?) 
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Every  one  of  the  above  quantities  is  dimensionless.  The  coefficients  a^# 
which  are  also  dimensionless,  are; 
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a,  to  a  inclusive  are  coefficients  which  depend  on  the  gas  and  seed. 
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Referring  now  to  equations  (26.19)  and  (26.20),  an  expression  for  polytropic 
efficiency  will  be  developed.  In  an  adiabatic  expansion; 
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Thus  if  1>  c  ,  then  n  >  K„.  However,  comparison  of  efficiencies  should 
always  be  made  on  the  basis  of  similar  power  densities 
to  obtain  higher  efficiencies  by  reducing  the  rate 


e  <  '•  Alternatively,  comparison  of  power  densities  may  b 
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Referring  to  equation  ^26.19  ,  power  density  is 
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Making  the  above  substitution,  the  following  inequality  is  obtained; 


f  ( 1  +  p')  >  1  + 


It  would  sees  that  f  > 1  is  a  possible  answer  but  the  cancellations  made 
above  are  legitimate  if  1  -  €  is  negative,  therefore  the  answer  is  in  the  negative. 

Another  possibility  remains:  it  is  that  although  the  power  density  wou’f. 
seem  to  be  reduced  on  the  above  argument,  equation  (26.21)  shows  that  f  >  1  v^.ild 
mean  greatly  increased  electron  temperature.  Following  a  Similar  procedure  to 
thar  adopted  above,  the  possibility  that 

2 

2  2  9 

O’Ky)  (1  +  f  5  >  (1-7?)  (1+£")  will  be  investigated. 

Substitution  for  rj  from  the  expression  above  leads  to 
0  >  (1  -e)2 

which  is  clearly  impossible. 


Thus  it  is  clear  that  if  efficiency  is  to  fee  sacrificed  in  order  to  increase 
power  densities  (by  virtue  of  increased  electron  temperature)  there  is  no  better 
way  to  acnieve  this  end  than  by  use  of  the  segmented  electrode  generator,  whether 
this  be  series  connected  or  not. 

26.3  HAIL -MODE  OPERATION 


In  the  case  of  the  Hail  generator,  or  ether 
it,  certain  interesting  effects  (which  may  settle 
temperature  cne  way  or  the  other)  may  be  observed, 
qualitative  explanation  of  such  effects; 
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the  issue  of  enhanced  electron 
The  following  is  a 


For  a  Hail  generation  equations  (26.21)  and  (  6.2*3)  become 
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and 
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and  if  the  variation  of  @  with  Te  is  assumed  to  be  slow  compared  with  “Jr  which  varies 
&t >  rox^mately  as  Te  then  the  fora  of  the  above  equations  between  Te  and  €  is  as 
follows ; 
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and 
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b  (T  -  1) 


.....(26.33 


Possible  solutions  are  represented  o y 
t  versus  Te  which  are  drawn  below  with  Hall 
as  parameter. 


intersections  of  the  above  curves 
current  (proportional  to  ( 1  -  < . ) 


of 


1  .0 


f 


0 


From  the  diagram  above  there  are,  in  general,  one  or  three  solutions  of  the 

equations  (  (aj,  (b) ,  (c)  )#  however  it  is  found  that  solutions  at  low  T  or  f 

e 

are  always  unstable  and  only  the  highest  are  stable.  However,  for  best 
efficiency  a  Hall  generator  should  always  have  e  <0.5  so  the  result  of  this 
instability  is  high  power  density  (high  T  )  but  very  poor  efficiency. 

26.4  DISCUSSION 


The  equations  developed  in  this  report  are  in  a  convenient  form  for  forward 
integration  on  a  digital  computer.  In  fact  they  have  already  been  programmed  and 

a  number  of  preliminary  calculations  have  been  carried  out.  It  is  intends  that 
the  scope  of  the  programme  will  be  extended  further  to  include  magnet  and 
compressor  power.  When  this  has  been  done  it  will  be  possible  to  calculate  (by 
trial  and  error  the  optimum  values,  from  the  thermal  efficiency  standpoint,  of  all 
the  independent  parameters  of  the  system. 


For  example,  if  optimum  3eed  fraction  is  calculated  by  maximising  <?U"  or  by 
some  similar  device,  the  result  will  inevitaoly  be  a  function  of  pressure, 
temperature,  faux  density  etc.  However  if  simultaneous  optimisation  of  all 
the  independent  parameters  is  carried  out  as  proposed,  the  optimum  seed  fraction 

then  be  the  best  for  the  optimum  pressure  and  temperature  range,  and  therefore 
for  the  unit  of  highest  overall  efficiency. 


26.9 


*he  j-re^inun&ry  calculations  earned  out  on  the  computer  have  drawn  attention 
in  vhe  ease  ViL  Hall  type  generators  to  a  phenomenon  akin  to  voltage  breakdown 
whicr  nas  not  yet  been  fully  investigated;  it  is  &  consequence  of  electron 
temperature  elevation  due  to  the  presence  of  the  magnetic  field  and  its  effect 

is  to  reduce  the  efficiency  of  the  expansion  considerably.  This  could  preclude 
the  use  of  Haul  generators  altogether. 
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CHAPTER  2  ? 


OPTIMIZATION  OF  LARGE-SCALE  NUCLEAR  l 


B.C.  Lindley  and  I.R.  M.  Nab 


27.1  INTRODUCTION 

Research  related  to  the  eventual  use  of  magneiopiassadynamic  (MFD)  generators, 
for  large-scale  electrical  power  production  ,  d  for  other  purposes,  Is  currently 
m  progress.  The  coal-;  oil-  or  gas-fired  open  cycle  system  utilises  the  alkali 
metal  seeded  combustion  products  as  the  working  fluid  while  the  closed  cycle 
system,  with  a  nuclear  reactor  as  the  heat  source,  uses  an  alkali  metal  seeded 
inert  gas. 

*  e  ^esiUii*  ( n  i  tn  i.- ne  vowes  ^  ionisation  potent  tar  cd  any  material } ,  although 
expansive,  is  adopted  for  closed  cycle  systems.;  because  of  the  difficulty  of 
complete  seed  recovery  in  the  open  cycle  system,  a  cheap  potassium  compound  is 
employed.  For  thermal  ionisation,  the  electrical  conductivities  of  cesium-seeded 
inert  gases  are  comparable  witn  those  of  potassium-seeded  combustion  products  at 
temperatures  about  500°C  higher,  so  that  the  closed  cycle  system  has  an  immediate 
advantage  over  the  fired  system.  In  inert  gas  plasmas  there  also  exists  the 
possibility  of  achieving  high  electrical  conductivities  through  magnetically- 
induced  extra  thermal  ionization1  at  quite  moderate  gas  t  emperatures  (1200°-1800°C). 
If  operation  at  such  low  temperatures  car  be  proven,  the  direct  was— cooled  nuclear 
reactor  closed  cycle  scheme  appears  attractive. 

In  a  nuclear  cycle  selection  of  the  operating  pressure  level  presents 
difficulties:  the  reactor  should  operate  at  as  high  a  pressure  ( typically  10-20atm) 
as  is  technically  and  economically  feasible  (to  improve  heat  transfer  for  a 
permissible  pressure  3  ;ss  and  pumping  power)  whereas  the  MPD  generator  benefits 
in  plasma  electrical  conductivity  by  operation  at  ow  pressures  (l  atm  or  less). 

27.2  GENERATOR  PERFORMANCE 

27*2.1  General 

In  optimizing  the  specific  power  in  an  MPD  generator  the  usual  procedure 


*  International  symposium  on  magne tchydrodynaiPic 
Paris,  6  -  1C  July  I9ch, 


generation 


is  to  maXiiu ije  at*  assuming  thermal  ionization  conditions;  subsonic  Mach  numbers 
(about  u.n  for  helium— cesium)  result. 

With  *»S»ticaUy-induced  ionisation  the  ratio  of  electron  to  gas  temperature 
(T0/T)  increases  witn  increasing  Mach  number  and  Hall  number  of  the  flow.  The 
most  important  ionisation  process  in  an  MPD  generator  is  through  electron 
collisions;  consequently,  if  Tq/T>  1  oar.  be  achieved,  the  restriction  on  flow 
velocity  (imposed  by  thermal  ionization  at  the  low  static  temperatures  associated 
with  high  Mach  numbers)  no  longer  applies,  and  there  is  the  possibility  of 
operating  supersonically  with  high  specific  power.  In  this  case  the  expansion 
pressure  ratio  in  the  accelerating  nozsle  will  be  igh,  so  that  the  reactor 
operates  at  a  much  higher  pressure  level  than  the  MPD  duct. 

27*2.2  Theory 

Assuming  that  the  flow  velocity  is  small  at  entrance  tc  the  nozzle  (M  «1  ) 
the  static  temperature  and  pressure  are  closely  approximated  by  the  stagnation 
conditions.  After  an  isentropic  expansion  in  the  nozzle,  the  static  temperature 
and  pressure,  and  velocity  are  as  follows  (and  see  Fig.  27.1 ): 


T  =  T  1  +  u  2 

2  °L  2  *2 


P2  =  YYV 


r/r-i 


(27.1) 


(27.2) 


V2  =  M2UB?2)r 
In  the  generator  duct, 

r-i 


2k  .fc 


*  T 


(27.3) 


(27.4) 


Selecting  a  typical  value  of  may  be  evaluated,  and  the  average 

generator  conditions  may  be  approximated  by  using 

_  _  P4  +  P2 

P3  "  2 

ai.d  hence,  from  equation  (27.4),  the  corresponding  T  . 

.  i  3 


Then  M3  =  M2  (T^/T  )l/2 


♦  (27.p) 


np3  =  Py  3  *.*.*  (27.6) 

Magneti^aily-induced  extra-thermal  ionization  may  be  investigated  using 
the  equation  derived  by  Hurwitz  et  al  for  a  segmented-electrode  generator  (all 
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the  quantities  given  below  apply  to  region  3  of  the  generator) : 


T@  r(i-K)2lO?  d 

ST“  =  1  +  — - . .  {  >  ■>  7  i 

An  iterative  method  is  required  for  the  solution  of  this  equation  since,  in 
general  and  in  the  cases  considered  here,  is  a  function  of  T  ,  excluding 
analytic  solution  of  equation  (27.7). 

Then 
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Having  obtained  T  ,  the  corresponding  electron  concentration  is  obtained 
x  e 

from  Saha’s  equation": 


n  =  T 
e  e 


in  rl 'oh  it  is  assumed  that  the  fractional  ionisation  (ne/n  )  -  ne/^  )  is 
determined  by  Saha's  equation  at  the  electron  temperature  and  also  that  the 
fractional  ionization  is  such  that  (ne/r.  )  <<1. 

Jr 

For  thermal  equilibrium  ionisation,  equation  (27.11 )  is  used  to  give  the 
electron  concentration,  the  electron  temperature  being  replaced  by  the  gas 

temperature , 

The  scalar  electrical  conductivity  is 
2 

n  e 

o  -  — —  (is 


m  v 
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(27.12) 


and  the  specific  power  for  a  segmented-ele  trode  generator  is 


P,  = 

) 


P*V  2b2k(i-k) 

J  c. 


(27.13) 


27.3  SPECIFIC  POWER 

27.3.1  -aerating  conditions  and  data 

The  effects  of  reactor  operating  pressure  on  the  specific  power  (electrical 
power  output  pei  unit  volume  of  generator  duct)  in  a  segmented-electrode  constant- 


57  s 

Z-  X  •  V 


velocity  subsonic  MPD  generator  have  been  examined  using  the  simple  quasi -one 
dimensional  theoretical  model  described  in  Section  27*2.  Helium-cesium  and 
argon-cesium  mixtures  are  investigated  for  thermal  equilibrium  and  magnetically- 
induced  non-thensal  ionisation.  Wall  friction  and  heat  transfer  losses  are 
neglected.  In  a  first  series  of  calculations,  a  generator  inlet  Mach  number 
of  0.6  is  assumed;  in  a  second  series  the  Mach  number  is  increased  to  supersonic 
values. 


Data  for  the  calculations  are  given  in  Table  27*1*  The  electron-helium  atom 
collision  cross  section  is  obtained  from  a  previous  survey  of  experimental 

measurements  .  The  electron-argon  atom  cross  section  is  obtained  from  an 

6 

approximate  curve  fit  to  Erodes  results  .  Although  several  measurements  of  the 
electron-cesium  atom  cross  section  have  recently  been  made,  there  still  exists 
considerable  discrepancy  and  an  average  is  adopted  (see  chapter  22),  The  ionic 
mobilities  are  from  the  measurements  of  Chanin  and  Biondi  (Ce+  -  He  and  Cs+  -  A) 
and  the  recent  measurements  of  Chanin  and  Steen^  (Cs+  -  Cs);  the  total  mobility 
is  obtained  from 


1 


i 


(27.14) 


2?«3»2  Results 

Using  equations  (27*1 )  to  (27,14)  and  the  data  given  in  Table  27.1,  a 
ftgasus  Autocode  programme  has  been  constructed  to  perform  the  calculations. 

Specific  powers  for  a  generator  inlet  Mach  number  of  0.6  are  shown  in 
Fig».27.2  and  27-3 ,  and  as  a  function  of  Mach  number  in  Figs. 27.4  and  27*5* 

The  ratio  (enhancement  in  specific  power)  is  shown  in  Fig.  27.6  for  the 
parameters  cf  Figs.  27*2  and  27.3* 

27*4  DISCUSSION 

27.4.1  Constant  Mach  number 

In  Fig.  27.2,  both  for  thermal  equilibrium  and  elevated  elctron  temperature 
conditions,  maximum  specific  powers  are  obtained  at  particular  value*  of  reactor 
pressure,  depending  on  the  applied  magnetic  field.  The  maxima  occur  because 
of  two  opposing  effects:  as  pQ  decreases,  and  other  parameters  remain  constant, 
the  electrical  conductivity  increases;  and  as  pQ  decreases  ion  slip  becomes 
increasingly  important  in  the  specific  power  equation  (27*13),  giving  a  reduction 
in  power. 

The  maxima  for  equilibrium  and  non-equilibrium  conditions  occur  at  different 
pressures  because  of  the  different  form  of  dependence  of  the  electrical  conductivity 
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and  ion  slip  terms  in  each  case.  Over  the  central  pressure  range  the  specific 
power  for  nor -equilibrium  conductivity  is  very  much  higher  than  for  thermal 
equilibrium,  but  at  very  high  and  very  low  reactor  operating  press'*"  ■'s  non- 
equilibrium  effects  are  relatively  unimportant.  The  curves  for  different  magnetic 
fields  cross  at  certain  pressures  because  of  the  dominance  of  ion  slip  so  that  an 
optimum  magnetic  field  exists  in  this  low  pressure  regime.  At  typical  reactor 
operating  pressures  (>  10  atm)  the  specific  power  increases  monotonically  with 
magnetic  field  over  the  range  considered.  The  loci  of  points  are  shown  for  which 
<p  m  10  and  100$  at  high  operating  pressures;  as  the  magnetic  field  is  increased 
the  reactor  may  be  operated  at  higher  pressures. 

Por  a rgon-cesium  a  lower  range  of  pressure  and  a  seeding  fraction  of  0.01 
are  considered.  In  the  argon-cesium  mixture  the  electron-argon  atom  cross 
section  varies  appreciably  with  electron  energy  («  0.02  represents  an  optimum 
value  for  helium-cesium  in  which  the  electron-helium  atom  cross  section  is  nearly 
independent  of  electron  energy);  x  *  0.01  represents  an  avexage  value  to  give 
maximum  conductivity  over  the  relevant  energy  range.  Comparing  Fig.  27.3  with 
Pig.  27.2  the  maximum  specific  power  is  lower  and,  at  high  reactor  operating 
pressures,  <f>  is  greater  for  argon-cesium  than  for  helium-cesium,  because  of 
the  smaller  electron-argon  atom  collision  cross  section.  However,  even  with 
non-equilibrium  ionization,  the  specific  power  at  10  ata  for  helium-cesium  is 
greater  than  for  argon-cesium,  primarily  on  account  of  the  lower  velocity  for 
a  given  flow  Mach  number  in  argon. 

It  has  been  assumed  that  the f ractional  ionization  (n Jn  )  is  small, 
allowing  electron-ion  collisions  to  be  neglected  and  adoption  of  the  simplified 
form  of  Saha's  equation  (27.11 ).  The  regions  for  which  n ^n&  >  0.01  are 
indicated  in  Figs.  27.2  and  27.3  and  it  is  probable  that  specific  powers 
attainable  in  practice  a  re  lower. 

27.4.2  Variable  Mach  number 

The  effect  of  varying  the  generator  inlet  Mach  number  is  indicated  for 
helium-cesium  and  argon-cesium  for  reactor  operating  pressures  of  1  ,3  and  10  ata; 
for  all  pressures  the  maximum  specific  power  for  equilibrium  conductivity  occurs 
at  a  flow  Mach  number  of  about  0.5. 

In  Fig.  27.4,  at  the  lowest  pressure,  non-equilibrium  ionization  has 
considerable  effect,  and  up  to  a  Mach  number  of  1  .4  the  specific  power  increases 
monotonically  with  Mach  number.  At  10  ata  and  low  Mach  numbers  the  non-equilibrium 
curve  is  close  to  the  equilibrium  but  at  Mach  1  .4  the  specific  power  for  non¬ 
equilibrium  is  about  seven  times  greater  than  for  equilibrium,  although  the 
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absolute  magnitude  is  less  than  the  maximum  both  for  equilibrium  and  non-equi librium* 
At  intermediate  pressure  (—3  ata)  initially  the  non-equilibrium  curve  is  close  to 
that  for  equilibrium  (because  the  static  temperature  decreases  more  rapidly  than 
T Jl  increases  in  equation  (27*7)  but,  at  high  Maob  numbers,  non-equilibrium  effects 
become  important  and  the  specific  power  increases-  rapidly?  For  3  ata  pressure 
there  is  no  maximum  in  the  non-equilibrium  specific  power  over  the  range  considered 
but  for  3  <P  <1 0  ata  the  non-equilibrium  power  will  follow  the  maximum  in  the 
equilibrium  power  and  then  diverge  to  a  greater  value  at  high  Mach  numbers. 


For  argon-cesium  (Fig.  27*5)  there  are  several  significant  differences.  In 

the  range  of  parameters  being  considered  ion  slip  is  appreciable,  and  has  the 

effect  that  greater  specific  power  is  achievable  with  pQ  =  3  ata  than  with  1  or  10 

ata  for  equilibrium  ionization  (also  apparent  in  Fig.  2/.3)*  With  non-equilibrium 

ionization  at  p  =  1  and  3  the  specific  power  initially  increase  rapidly  with 

Mach  number.  For  p  =  1  the  rate  of  increase  of  specific  power  falls  at  high 

o 

Mach  number,  eventually  showing  a  maximum.  The  maximum  arises  for  two  reasons; 
as  the  Mach  number  increases  the  static  gas  pressure  and  temperature  in  the  MPD 
generator  decreases  but  ratio  T increases;  at  low  Mach  numbers  (for  pQ  =  1  ata) 

T Jl  increases  sufficiently  rapidly  to  produce  an  overall  increase  in  electron 
temperature,  whereas  at  high  Mach  numbers  the  static  gas  temperature  becomes  very 
low  and  the  rate  of  increase  of  T is  restricted  by  ion  slip.  For  pQ  =  10  ata 
a  maximum  and  a  minimum  occur  in  the  non-equilibrium  power.  A+  Mach  1.4-  the 
specific  power  is  approximately  the  same  as  fcr  the  maximum  at  -Mach  0.5,  but  is 
continuing  to  increase  with  Mach  number.  The  actual  magnitude  of  specific  power 
is  in  this  case  greater  than  that  for  helium  at  10  ata,  even  though  the  flew 
velocity  is  considerably  smaller. 


27.3  CONCLUSIONS 

1  In  helium-cesium  mixtures  at  typical  reactor  operating  pressures  (10  ata) 
magnetically- induced  extra-thermal  ionization  will  not  be  of  consequence 
in  a  typical  subsonic  segmented-electrode  MPD  generator  unless  magnetic 
fields  of  greater  than  5T  are  used.  Lower  reactor  operating  pressures 
would  enable  greater  specific  powers  to  be  achieved  with  the  same,  or 
lower,  magnetic  field  strengths. 

2  In  argon-cesiimi  mixtures  considerable  non-equilibrium  effects  will  occur; 
typically,  for  pQ  =  10  ata,  B  =  5T,  the  non-equilibrium  specific  power 

is  approximately  twice  the  equilibrium  specific  power.  The  actual 
magnitude  of  th°  specific  cower  in  argon-cesium  is  less  than  that  in 
helium-cesium  at  a  given  Maoh  number,  because  of  the  -aw  attainable 
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flow  velocity  in  argcn, 

3  As  the  operating  Mach  number  of  tne  MPD  generator  is  increased  non* 
equilibrium  efl  ©cts  become  increasingly  important,  out  for  helium- 
cesium  mixtures  the  effects  are  not  significant  for  a  reactor  operating 
pressure  of  10  ata  even  »?ith  a  Macn  number  of  1  .if.  For  a  rgon— cesium 
mixtures  the  specific  power  at  Mach  numbers  greater  than  1  .if  increases 
beyond  the  maximum  at  M  «  0*5 $  and  the  actual  specific  power  is  greater 
than  for  helium— cesium,  In  both  cases,  if  the  reactor  operating 
pressure  is  reduced,  for  example  to  3  ata,  the  specific  power  increases 
very  rapdily  with  Mach  number,  the  effect  being  most  significant  in 
argon-cesium  mixtures. 
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TABLE  27*1 


DATA  (MKS  UNITS) 


Universal  Constanta 

a  =  9.108.10  31 

*  -19 

e  -  1.601.10  7 

h  a  6.625.1 0-34 

k  =  1.36.10"23 


All  calculations 

=  0.85 
k  =  o.eo 
8  =  1 


Vi  =  3.893 

Tq  =  2000°K 

W1-3  20 

%-Cs*  300.1 0‘20 

Fig.  27.2 «  27.3  and  27.6 

B  =  5, 3, 2, 1,0.5 

^  =  0.6 

=  0.01  -*■  20  ata  (helium-cesium) 
pQ  *  0.05  •*  20  ata  (argon-cesium) 


Helium-cesium 


r  S 
R  = 


X  * 
"i  = 
^e-He 


1 .658 
2 .075.1 03 
0.02 


i.73-1022/^ 


=  6.10 
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Argon-cesium 


r  =  1  .66 8 

R  =  2.079.1 02 

X  =i.io"2 

hi  =  5.92.1022/nA 


r  -22 

1.12.10 


^e-A 


—t 


1  .2.10 


-20 


Fig.  27.4  and  27.5 
B  =  3 

P0  =  1,3,10  ata 

Mg  a  0.2  -  1  .4 
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SPECIFIC  POWER  OF  A  SUBSONIC  (m=  0  6*),  SEGMENTED  ELECTRODE 


HELIUM  CESIUM  MPD  GENERATOR  AS  A  FUNCTION  OF  REACTOR  OPERATING 
I  PRESSURE 
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NOMENCLATURE  (CHAPTER  27) 


B  magnetic  field  strength 

e  electronic  charge 

Sg4a  statistical  weight  of  electron,  ion,  atom 


h  Planck’s  constant 

k  Boltzmann’s  constant 

K  load  factor 

m  electronic  mass 
© 

M  Mach  number 

n&  atom  concentration  of  seed  gas 
ne  electron  concentration 

atom  concentration  (of  parent  gas) 
p  static  pressure 
PQ  stagnation  pressure 
P  specific  power 

q  collision  cross  section 
R  gas  constant 
T  static  ga^  temperature 
Tq  stagnation  gas  temperature 
electron  temperature 
V  flow  velocity 

ionization  potential 
electronic  Hall  coefficient 
P^  ionic  Hall  coefficient 
X  ratio  of  specific  heats 
a  inelastic  collision  factor 
<  Te/1l605 

T).  local  isentropic  efficiency 

XTt)  -  HD 

*  pTt) 

ve  electronic  collision  frequency 
ionic  mobility 
electrical  conductivity 
seeding  fraction 


Subscripts 

o  stagnation  conditions 

1  nozzle  inlet  conditions 

2  MPD  generator  inlet  conditions 

3  average  MPD  generator  conditions 

4  MPD  generator  exit  conditions 
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CHAPTER  28 


POWER  CONVERSION  IN  SPACE 


by 


I.R.  McNab 


28. 1  INTRODUCTION 

Electrical  power  is  undoubtedly  the  most  important  single  subsystem  of  any 
space  vehicle,  being  vital  to  communications,  data  gathering  and  processing,  life 
support  and,  directly  or  indirectly,  to  propulsion.  In  many  cases  the  constraints 
imposed  on  the  vehicle  design  and  operation  by  the  power  system  are  severe.  ?or 
example,  the  power  system  can  account  for  a  large  fraction  of  the  payload  weight, 
it  can  influence  the  thermal  control  system  and  it  may  introduce  crew  hazards. 

In  addition  it  may  also  require  solar  attitude  control,  involve  complicated 
prelaunch  procedures  and  recti y  limit  the  missions  the  vehicle  can  perform. 

Criteria  to  be  considered  in  the  selection  and  operation  of  space  power 
systems  include  the  power  to  weight  ratio,  reliability,  compatibility  and  degree 
of  integration  with  the  spacecraft,  launch  and  operation  techniques,  resistance  to 
space  environment  and,  for  manned  spacecraft,  the  possible  crew  hazards. 

Present  and  past  power  requirements  for  satellites,  spaceprobes  and  manned 
missions  have  been  filled  by  the  use  of  solar  cells  and  batteries,  either 
separately  or  together,  and  radioisotope-thermoelectric  devices;  the  power 
levels  generally  being  much  less  than  1  kWe.  While  such  systems  have  performed 
satisfactorily  they  become  less  attractive  when  mere  advanced  missions  are 
considered,  mainly  because  of  power  to  weight  ratio  considerations.  Estimated 
power  requirements  for  space  vehicles  over  the  next  few  ye ars^,if,  (Pig.  28.1)  show 
that  electrical  propulsion  systems  will  eventually  require  large  amounts  of  power. 
Missions  utilizing  electric  propulsion  may  be  arbitrarily  divided  into  four 
classes,  typical  examples  of  which  are,  in  probable  order  of  occurrence; 

(A)  earth  satellite  attitude  control  and  station  keeping,  and  earth  satellite 
orbit  transfer; 

(B)  solar,  non-ecliptic,  planetary  and  interplanetary  probes; 

*  I.R,  McNab.  International  Research  and  Development  Compary,  Newcastle  upon  Tyne,  f 
(Paper  presented  at  the  I.E.E,  Colloquium  on  ’Electrical  Methods  of  Propulsion  in 
Space',  London,  February  13th,  195.'+.) 

i  It  must  be  anphasissed  that,  in  the  absence  of  a  British  space  programme,  most  of  the 
information  contained  herein  has  been  obtained  from  American  sources,  particularly 
references  1 ,  2  and  3* 
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(C)  large  spacecraft,  Banned  or  unmanned,  of  the  lunar  ferxy  type;  and 

(D)  manned  planetary  missions. 

The  electrical  power  supply  requirements  for  these  missions  are  typically 
5-30  kWe,  30  -  300  kWe,  1  -  3  MWe  and  20  -  40  MWe  respectively.  Mission 
durations  are  unlikely  to  be  less  than  a  week  and  will  generally  be  several  months, 
(it  is  mainly  this  condition  which  leads  to  superiority  of  electrical  propulsion 
compared  with  chemical  or  chemical/nuclear  propulsion  systems  for  interplanetary 
missions.)  Plight  tests  of  electric  propulsion  systems,  and  possibly  of  Class  A 
missions,  are  likely  to  occur  during  1964.  The  early  introduction  of  electric 
propulsion  for  Class  B  missions  (by  about  1970)  is  currently  anticipated,  although 
this  is  dependent  to  a  considerable  extent  on  the  availability  of  electric  power 
supplies. 

At  present  a  great  diversity  of  space  power  systems  are  under  investigation, 
covering  many  different  power  levels  and  mission  times.  Pig.  28.2  ahows 
approximately  the  power  levels  and  mission  times  for  different  power  systems^; 
for  electric  propulsion  systems  with  high  power  levels  and  long  mission  durations 
the  mid  and  upper  right  hand  region  of  this  figure  is  of  greatest  importance. 

To  give  the  background  to  the  development  of  the  advanced  conversion  systems 
currently  under  consideration  a  brief  review  of  some  of  the  systems  shown  in 
Fig.  28.2  will  now  be  made;  some  of  the  advantages  and  disadvantages  of  each 
system  will  be  mentioned. 

28.2  PRESENT  POWER  CONVERSION  SYSTEMS  AND  ASSOCIATED  RESEARCH 

28.2.1  Solar  cell  systems 

Solar  cells,  with  or  without  associated  storage  batteries,  have  been  the 
most  widely  used  power  systems  to  date.  These  have  grown  rapidly  from  the  small 
wattage  devices  used  in  early  satellites  to  arrays  such  as  that  for  the  OAO 
satellite  (Orbiting  Astronomical  Observatory)  which,  if  deployed  so  that  all  cells 
were  simultaneously  illuminated  would  generate  about  1,5  kWe,  the  average  power 
being  about  600  We. 

The  main  problem  of  solar  cells  is  the  performance  degradation  occurring 
due  to  the  nuclear  particle  radiation  associated  with  the  earth  radiation  belts 
and  solar  flares.  A  change  has  recently  been  made  from  the  diffused  boron  P/N 
silicon  cells  originally  used  in  U.S.  Satellites  to  the  more  radiation  resistant 
diffused  phosphorous  N/P  cells.  The  most  probable  method  of  increasing 
performance  will  arise  from  the  development  of  new  materials.  One  possible 
method  of  increasing  the  power  level  of  present  systems  is  to  use  reflectors  to 
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increase  the  amount  ei 


offset  the  advantages  gained. 


28.6.2  Batter-/  systems 
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Batteries  will  undoubted!/  continue  to  be  of  importance  for  suxilia; 


power 


requirements  although  the  introduction  of  fuel  cells  may  well  reduce  the  number  of 
applications  of  batteries.  Fuel  ce1is  are  scheduled  for  the  Gemini  and  Apollo 


spacecrj 


.ft.  Attractive  features  of  these  self-contained  svftt.sms 


ems  induce  tnexr 


high  efficiency,  availability  of  water  as  a  by-product  (for  hydrogen/ oxygen  cells) 
and  lack  of  need  for  attitude  control  or  shielding.  Regenerative  fuel  cells  and 
batteries  will  undoubtedly  be  important  in  lunar  and  planetary  exploration  for 
life  support  and  vehicle  power.  Used  in  conjunction  with  a  main  energy  depot 
involving  a  nuclear  or  solar  power  plant,  these  systems  will  provide  high  mobility 
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and  efficie ncy ' * " , 

28.2.3  Radioisotope  systems 

Radioisotope  thermoelectric  generators  have  been  used  and  are  planned  for 
several  space  missions.  A  2,?  W  plutonium  238  fuelled  thermoelectric  generator 
(SNAP  3)  was  launched  in  Transit  4A  in  June,  1$6l,  and  was  still  operating 
satisfactorily  two  years  later.  A  larger  generator  (SNAP  9A;  25  We)  was  launched 

in  September,  1963,  in  Transit  5;  this  has  a  design  lifetime  of  five  years. 

Larger  devices  (200  -  500  We)  are  planned  for  Voyager  and  Nimbus  missions  and 
design  studies  for  3  -  5  kW  systems  for  use  in  manned  space  stations  are  under  way. 
Most,  if  not  ail,  of  the  radioisotope  systems  require  the  relatively  scarce  and 
expensive  a  emitters,  because  of  spacecraft  and  handling  considerations.  In  view 
of  their  high  cost  isotope  systems  will  probably  only  be  considered  where  their 
special  attributes  -  invulnerability  to  external  radiation,  long  life  and  potential 
reliability  -  are  important.  Isotope  system 3  current ly  use  lead  telluri.de  or 
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germanium -bismuth  telluride  conversion  elements,  operating  at  450’  -  500  C. 

Future  converter  technology  will  be  devoted  towards  the  more  efficient  use  of 
isotope  heat  sources  through  the  development  of  thermoelectric  systems  having 
higher  working  temperatures  and  efficiencies  (such  as  germanium-silicon  alloys) 4 
ultimately  thermionic  conversion  will  be  employed. 

28.2.4  So la r- iynaml c  sy st en s 
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mentioned.  A  more  efficient  use  of  reflectors  is  for  solar-thermionic  or  solar- 
dynamic  conversion  systems.  In  solar-thermionic  systems  a  mirror  is  used  to 


focus  the  solar  ener-;  into  a  small  cavity  formed  b-  the  emitters  oi 
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thermionic  diodes. 


Devices  of  this  type  have  bee.',  ground  tested  and  have  yielded 


efficiencies  of  about  1%. 


In  solar -dynamic  systems  a  mirror  is  used  to  focus  the  solar  energy  into  a 
cavity  containing  the  working  fluid.  Ranking  cycle  systems  are  generally 
considered  for  such  devices,  thus  the  working  fluid  is  evaporated,  produces 
electricity  in  a  turbo-alternator  and  Is  then  condensed  in  a  condenser/ radiator. 
Considerable  effort  is  currently  being  put  into  such  systems  in  the  United  States, 
particularly  on  the  NASA  Sunflower"  and  USAF  "Astec"  systems,  the  power  outputs 
being  about  15  kWe,  The  Sunflower  system,  which  is  based  largely  on  SNAP  2 
t-urbo^ Icemator  technology,  uses  a  32*5  ft  diameter  petal-type  mirror  and  mercury 
as  the  working  fluid.  Some  of  the  major  components  of  the  system  have  already 
been  tested. 


To  date  the  accurate  focussing  required  for  solar-thermionic  and  solar- 
dynamic  systems  is  only  available  with  precise,  one  piece,  rigid  mirrors. 
Considerable  effort  is  presently  being  made  on  the  development  of  foamed, 
inflatable  and  unfurlable  petal  systems  which  can  be  launched  without  difficulty. 
Table  28.1  summarises  the  advanced  mirror  development  programmes  in  the  U.S.  ♦ 

While  mirrors  probably  eliminate  the  radiation  damage  problem  of  solar  cells, 
little  is  known  at  present  of  their  performance  characteristics  in  a  space 
environment,  although  it  may  be  anticipated  that  micrometeorite  damage  will  have 
a  significant  effect  on  long  term  mirror  efficiency. 

For  power  systems  up  to  about  100  kWe  solar  dynamic  systems  have  potentially 
higher  power  to  weight  ratios  and  lower  costs  than  other  systems.  The  absence 
of  radiation  and  shielding  problems  make  these  3olar  powered  systems  attractive 
for  manned  space  stations.  The  energy  conversion  equipment  for  the  solar- 
dynamic  systems  :  j  similar  to  that  under  investigation  for  nuclear-dyr amis  systems 
and  thus  offers  opportunity  for  in-flight  demonstration  and  testing  without  the 
costs  and  hazards  associated  with  reactors. 

28,2,5  Reactor  systems 

Reference  to  the  categorisation  of  eiectric  propulsion  missions  in  the 
Introduction  to  this  report  and  to  Fig.  28.2  shows  that,  of  the  systems  this  far 
discussed,  only  solar  collectors  are  likely  to  provide  the  power  levels  and  long 
lifetimes  required  for  electric  propulsion,  ana  then  only  for  Class  A  or  lower 
power  Class  B  missions.  for  long  term  missions  of  more  than  5C!  -  100  kfe,  nuclear 
power  sources  are  required. 

Several  power  conversion  systems  for  use  with  nuclear  reactors  are  currently 
under  consideration.  Probably  the  most  advanced  systems,  in  terms  of  the  expected 
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t light  test  times,  are  reactor-conversion  systems 
Auxiliary  Power)  series.  ^able  28.2  shows  the  current  status  of  several  SNAP 

develoument  programmes  .  The  early  SNAP  reactors  such  as  SNAP  3  and  9A  (already 

launched)  and  SNAP  10A  use  thermoelectric  conversion.  The  larger  SNAP  units  use 


dynamic  conversion  based  on  the  Rankin?  cycle,  the  power  output  being  a.c. 
cycle  for  the  SNAP  2  system  (Fig.  28. 3) b  use 


The 


-wo  fluid  loops 


-he  NaK  ( /8%  Na, 


22%  K)  reactor  coolant  loop  which  remains  liquid,  and  the  power  conversion  loop 
in  which  mercury  is  vaporised,  expanded  tn rough  a  turbine,  condensed  and  recycled. 
Onxy  one  moving  part,  the  combined  rotating  unit  (CRU)  is  used,  thus  both  pumps. 

The  ^ 


the  turbine  and  the  alternator  are  supported  on 
completely  separated,  hermetic  sealing  being  em] 
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Several  problems  follow 
from  the  use  of  the  CRU,  for  example,  the  lubricating  ability  cf  the  mercury,  the 
balance  of  component  efficiencies  at  a  common  shaft  speed,  thermal  distribution, 
shaft  alignment  and  balance,  and  protection  of  sensitive  parts  from  the  hot  mercury 
vapour.  The  start-up  procedure  as  complex  and  orbital  operation  will  be  further 
complicated  by  the  zero-gravity  condition  and  the  effects  this  will  have  on  the 
processes  in  the  boiler,  and  to  a  greater  extent,  the  radiator/condenser.  The 
SNAP  2  system  is  about  >3  ft  long  and  5  ft  diameter  at  the  base,  and  the  weight 
with  a  fully  shielded  reactor  is  about  1500  lb. 

The  more  powerful  SNAP  8  s^  stem  can  he  considered  as  a  logical  successor  to 
the  SNAP  10A/2  systems.  Several  components  of  the  smaller  systems  and  much  of 

the  accumulated  data  will  be  incorporated  in  the  design  and  construction  of 
SNAP  8.  The  aim  of  the  programme  is  the  development  of  a  long-life  generating 
system  capable  of  providing  at  least  35  ..  3  to  meet  anticipated  needs  for  large 
amounts  of  auxiliary  power  and,  possibly,  for  early  electrical  propulsion.  The 
basic  design  is  similar  to  SNAP  2,  thus  the  heat  produced  by  the  reactor  is 
removed  by  liquid  NaK  and  transferred  to  mercury  in  the  boiler.  The  mercury  is 
evaporated,  drives  the  turbo-alternator  and  is  then  condensed  and  recycled.  An 

additional  NaK  loop  is  provided  to  transfer  the  he^t  from  the  mercury  to  t-ht 

radiator. 

Many  problems  have  been  encountered  in  the  development  programme  for  SNAP  8 
and  the  current  design  emphasises  reliability  of  performance  at  the  expense  of 
weight.  One  severe  problem  is  that  of  turbine  operation.  Because  of  the 
requirement  of  low  specific  weights  for  space  power  systems  the  rotating  speed 
must  e  ve.j.  r.,gh.  ifu,&  condition,  together  with  the  low  v'.scosity  of  liquid 
metsx  lubricants,  creates  a  condition  of  bearing  operation  in  the  turbulent  region, 
which  is  not  we u  understood  at  present.  In  $uAR  8  the  problem  of  turbulent 
lubrication  by  liquid  metals  has  been  exchanged  for  that  of  seal  development  and 
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Si  i  ’•friction  bearings  lubricated  from  a  fourth  loop  circulating  an 

oti,'nic  fluid.  The  shaft  speed  has  been  dropped  from  the  original  40000  rpm 
to  12000  rpa. 

While  SNAP  8  will  provide  an  early  nuclear-electric  capability,  auch  higher 
powers  and  lower  specific  weights  are  required  to  realize  the  potential  of  eleotric 
propulsion  for  advanced  missions. 

SNAP  50  (for  which  some  data  are  given  in  Table  28.2)  is  one  approaoh 
currently  under  investigation  in  the  U.S.  to  the  problem  of  power  generation  in 
the  power  range  up  to  about  550  kWe,  The  system  is  basically  similar  to  the 
SNAP  2  and  8  projects  in  that  a  condensing  Ranldne  cycle  is  used  with  a 
turboelectric  power  conversion  system.  The  reactor  power  to  weight  ratio  is 
however  vastly  increased  by  using  a  uranium  carbide  fuelled  lithium  cooled  fast 
reactor  and  as  a  consequence  the  overall  power  to  weight  ratio  is  expected  to  be 
far  better  than  that  of  the  previous  SNAP  systems.  If  the  design  and  performance 
predictions  are  fulfilled  the  SNAP  50  system  will  enable  many  useful  missions  to 
be  carried  out  using  electric  propulsion  in  the  Class  B  missions  previously  given, 
although  the  cost  of  such  a  system  may  be  very  large. 

28*3  HIGH  POWER  (KEJAWATT)  SYSTEMS 

Per  power  levels  in  the  several  magawatt  (electrical)  range,  it  is  by  no 
means  clear  at  present  which  power  conversion  system  will  be  best.  Present 
hardware  commitments  for  nuclear  elects  c  spacecraft  powerplants  in  the  megawatt 
range  are  limited  to  liquid  metal  Bankine  cycles.  This  arises  from  early 

Compa  ion8  of  3uch  cycl«*  with  Biayton  closed  gas  cycle  power  systems2.  These 
early  comparative  studies  were  based  on  the  weights  of  the  components  inv  Ived 
and  mu so  to  a  large  extent  on  the  weight  and  type  of  radiator.  The  radiator 
represents  a  major  portion  of  the  weight  of  the  power  conversion  system  and  is 
also  the  largest  component,  being  therefore  difficult  to  launch  and  decreasing 
the  system  reliability  due  to  meteoroid  penetration.  (if  armour  is  provided  to 
decrease  the  probability  of  penetration  the  radiator  weight  increases.)  Thus  it 
is  des  rable  to  minimise  the  radiator  area.  Calculations  performed  for  Rankine 
end  Brmyton  turboeleetric  cycles  (Pig.  28,5  and  28.4  respectively  show  typical 
cycles)  show  the  superiority  of  Bankine  cycles  compared  on  the  basis  of  radiator 
area  (for  the  same  maximum  cycle  temperature)  or  on  the  basis  of  the  maximum 
working  temperature  for  the  same  area,  (Table  28.5  shows  results  due  to  Sutton7 
which  illustrate  this,  although  these  relate  specifically  to  powsr  conversion  by 
magnetoplasma dynamic  generators.) 

These  cycle  studies  have  spurred  the  development  of  practical  Bankine  cycle 
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of  the  many  difficult  problems  encountered  in 
devices  some  attention  is  being  given  to  the 
and  turbo— machine ry  have  been  proved  reliable 
whereas  alkali  metal  systems  have  yet  to  demo 
addition,  while  Table  28.3  indicates  that  gas 


the  development  of  alkali  metal 
Bray ton  cycle,  Gas  cooled  reactors 

over  several  years  of  operation 
nstrate  long  life  performance.  In 
cycles  require  higher  maximum 


temperatures  than  alkali  metal  cycles  for  the  same  radiator  area,  gas  turbines 
can  almost  certainly  operate  at  higher  temperatures  than  alkali  metal  turbines. 
Further  advantages  of  the  single  phase  gaseous  working  fluid  are  the  avoidance 
of  flow  problems  occurring  in  two  phase  systems  in  zero  or  low  gravity  operation 


and  the  elimination  of  freezing  du ring  snut-down  periods**. 


It  is  interesting  to  note  that  the  conventional  steam  cycle  has  been  proposed 
as  an  alternative  to  the  liquid  metal  Rankins  cycle.  While  s'  mi  is  limited  to 
low  cycle  temperatures  by  vapour  pressure  considerations,  its  .igh  heat-  capacity 
and  the  possibility  of  using  aluminium  construction  results  in  a  low  radiator 
weight  per  unit  area  and  a  reasonable  power  to  weight  ratio  of  the  overall  system. 


If  full  advantage  is  to  be  taker,  of  the  high  reactor  powers  available  for 
either  gas  cr  alkali  metal  cycles  it  is  necessary  to  use  one,  or  more,  of  the 
methoas  of  direct  conversion  of  heat  to  electricity  presently  under  investigation. 
The  two  methods  of  greatest  applicability  are  thermionic  diodes  and  magnetoplasma- 
dynamics  (MFD).  Both  these  methods  have  the  advantage  of  operation  at  higher 
temperatures  than  turbogenerators,  which  are  limited  by  stresses  in  the  rotating 
machinery.  While  both  systems  are  currently  in  the  research  and  development 
stage  it  is  probable  that  in  neither  case  will  the  efficiency  be  very  high.  For 
ground  based  power  stations,  both  diodes  and  MFD  are  currently  under  consideration 
only  as  ’topping’  devices,  working  in  conjunction  with  a  conventional  steam 
turboalternator  plant.  However  for  space  power  plants  the  ability  of  these 
systems  to  work  at  high  temperatures,  and  thus  to  reject  heat  at  high  temperatures 


may  reduce  the  radiator  size  to  such  an  extent  that 
competitive  as  the  main  power  conversion  system. 


thermionic s  and/or  MFD  become 


Several  different  cycle  arrangements  covering  different  diode  locations  and 

working  fluids  have  been  suggested  for  thermionic  systems.  probably  the  most 

general  categorization  of  thermionic  diodes  is  i«to  ’in-pile’  and  ‘out -of -pile1 
2 

systems  .  In  the  m-pils  configuration  the  heat  supplied  to  the  cathode 
(electron  emitter)  surface  comes  directly  from  the  reactor  fuel  and  heat  rejected 
by  the  anode  (electron  collector)  is  removed  by  the  coolant  passing  through  the 
reactor.  In  the  out-of-pile  systems  the  reactor  coolant-  heats  the  cathode  and 
the  rejected  heat  is  either  removed  by  a  second  flui1  or  radiated  directly  into 
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space. 

Thermionic  cell  performance  is  currently  temperature  limited  with  typic&x 
realistic  long-life  values  of  1750°C  at  the  cathode  surface  and  a  corresponding 
optimum  anode  temperature  of  about  700  G.  Thus  the  in-pile  systems  can  use 
reactor  coolant  temperatures  of  about  10QO°G  lees  than  the  out-of-pile  systems, 
for  similar  diode  performance  #  This  results  in  a  considerably  different  choice 
of  coolants  for  the  two  systems.  Inert  gases,  liquid  metals  and  boiling  or 
superheated  steam  can  be  considered  for  in-pile  systems,  For  out-of-pile  systems 
the  high  temperatures  necessitate  the  use  of  inert  gases  in  view  of  the  corrosion 
problems  involved  with  alkali  metals.  Several  possible  cycles  are  discussed  below. 

Fig.  28.5  shows  a  typical  single  stage  liquid  metal  vapour  system  employing 
in-pile  thermionic  diodes.  This  configuration  appears  capable  of  producing  a 
rather  better  specific  weight  than  a  turbo-electric  cycle  operating  at  the  same 
temperature.  The  main  problems  are  the  production  of  a  reactor  cooled  by  metal 
vapour,  thermionic  diode  reliability  in  the  radiation  environment  and  the  complex 
reactor  construction  required  to  contain  the  cells  and  their  associated  electrical 
connections  .  (Note  that,  although  a  single  stage  two  phase  cycle  is  shown  in 
Fig.  28.5,  a  practical  system  would  probably  use  a  two  stage  system  with  a  heat 
exchanger.) 

An  out-of-pile  liquid  metal  cycle  would  eliminate  the  disadvantages  of 
in-pile  diode  performance  and  reactor  design,  but  the  very  high  temperatures 
required  would  undoubtedly  cause  severe  corrosion  problems. 

In-pile  gas  thermionic  systems  of  the  type  shown  in  Fig.  23,6  suffer  from  the 
same  problems  of  diode  reliability  and  reactor  construction  as  the  liquid  metal 
vapour  systems  .  However  the  use  of  inert  gases  minimises  corrosion  effects  on 
the  cycle  components  and  may  enable  the  anode  surface  to  serve  directly  as  the 
coolant  wall.  (With  liquid  metal  vapour  or  steam  coolants  an  electrical  insulator 
which  is  also  a  good  heat  conductor  is  required  between  the  anode  and  the  coolant 
channel.) 

To  avoid  some  of  the  problems  inherent  in  the  in-pile  cycle  the  cycles  shown 

2 

in  Fig.  28.7  and  28.8  have  been  investigated  .  These  out-of-pile  gas  thermionic 
systems  have  many  advantages  (for  example,  fissionable  fuel  need  not  be  used,  the 
reactor  design  is  improved,  electrical  connection  problems  are  simplified)  but  a 
very  high  temperature  reactor  is  required  since  the  coolant  must  heat  the  diode 
cathodes  to  about  1750°C.  For  the  cycle  shown  in  Fig.  28.7  the  overall  efficiency 
is  improved  by  employing  a  turbine  to  drive  the  compressor.  The  results  of  a 
recent  theoretical  analysis  comparing  various  liquid  metal  and  ^s  turboelectric 
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and  thermionic  systems  are  shown  in  Table  2o.4,  The  temperatures  considered  in  the 
analysis  are 

liquid  metal  (maximum)  1130°C 
gas  (maximum)  1?30°C 
(turbine)  1280°C 

diode,  T  ,,  ,  -T  *  fit  n°n 
cathode  anoru 

The  criterion  on  which  the  systems  have  been  compared  is  the  specific  radiator 
weight,  which  is  the  single  item  of  greatest  weight  for  a  system  producing  more 
than  1  MWe,  The  results  show  that  for  a  thermionic  cell  efficiency  of  more 
than  about  13%  the  in-pile  thermionic  gas  system  will  be  superior.  It  is 
possible,  even  if  advances  are  made  in  the  efficiency  of  turboelectric  liquid- 
metal  systems,  that  if  corrosion  problems  cannot  be  solved,  gas  systems  will  be 
preferred  from  a  reliability  viewpoint,  even  though  less  efficient. 

The  second  direct  conversion  system  for  power  levels  above  1  MW«  is 
magne t op la sma dynami c  (MPD) ;  this  system  holds  promise  of  generating  larger  amounts 
of  power  with  better  efficiencies  than  thermionic  diodes.  MPD  generators  may 
be  divided  into  two  main  categories,  combustion  and  nuclear  reactor  powered; 
combustion  devices  are  usually  open  cycle  and  nuclear  devices  closed  cycle  systems. 
Direct  current  generating  combustion  devices  are  the  most  advanced  at  the  present 
time  so  far  as  power  output  is  concerned.  Substantial  amounts  of  power  have 
been  generated  by  several  systems,  for  example  the  AVCO  Mark  2  generator  has 
given  an  output  of  over  1  MWe  with  an  efficiency  of  about  7 %  Larger  devices 

are  expected  to  have  a  higher  efficiency  and  it  is  likely  that  land  based  MPD 
stations  generating  several  tens  of  megawatts  will  be  operational  within  the 
next  few  years.  However  for  space  missions  of  several  weeks  open  cycle 
combustion  systems  are  not  attractive  in  view  cf  the  large  fuel  supplies  required* 

Closed  cycle  nuclear  fuelled  MPD  systems,  while  conserving  the  working 
fluid,  suffer  from  the  temperature  limitations  of  present  and  foreseeable  reactor 
developments.  The  main  consideration  is  to  increase  the  electrical  conductivity 
of  the  working  fluid  to  a  high  level  (?•  30  mho/m) ,  In  both  combustion  and 
nuclear  reactor  heated  MPD  systems  the  main  working  fluid  is  seeded  with  an  alkali 
metal  to  enhance  conductivity.  Generally  only  a  few  atomic  percent  of  the  seed 
is  used  since  the  electrical  conductivity  has  a  maximum  with  respect  to  seed 
concentration  in  this  regime.  Thus  pure  alkali  metals  are  not  usually  considered 
for  rand  based  MPD  generators;  however,  radiator  si^e  and  weight  requirements  may 
alter  this  picture  in  the  3paee  field,  as  shown  in  Table  28.3. 
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Even  with  the  optimum  amount  of  alkali  m-^tal  seeding  it  is  unlikely  that 
reactor  temperatures  will  be  high  enough  to  make  ItPD  conversion  ieasibj.8  if 
the  electrical  conductivity  takes  its  thermal  equilibrium  value.  Consequently 
several  methods  of  producing  extra— thermal  ionisation  in  MFD  generat  rs  are 
currently  under  investigation.  These  methods  include; 

g 

ionization  by  high  voltage  electron  beams  ; 

8 

ionisation  by  high  current  arcs  ; 

q 

ionization  by  r,  f.  ; 

10 

photoionization  ; 
magnetically  induced  ionization 
ionization  by  p  rays  from  fission  products0;  and 
the  use  of  high  acceleration  to  produce  non-equilibrium  flows 
in  which  the  electron  density  does  not  decay  significantly 


11 


8 


12,13 

during  the  expansion  process  * 

At  present  none  of  these  methods  has  been  conclusively  demonstrated  in  an 
MFD  generator  producing  substantial  amounts  of  power,  although  several  devices 
are  under  construction. 

Fig.  28.9  shows  a  typical  seeded  inert  gas  nuclear  MPD  system  for  space 
purposes.  For  this  system  the  problem  of  corrosion  due  to  alkali  metals  is 
lessened  by  the  low  percentage  of  seed  present  and  the  fact  that  the  alkali  metal 
need  not  enter  the  reactor,  however  as  with  other  gas  cycles  the  radiator  area 
will  be  large  compared  to  the  metal  vapour  system. 

The  main  problem  of  MPD  generators  may  well  be  the  weight  of  the  magnetic 
field  coils.  Even  if  superconducting  magnets  are  available  the  weight  of  the 
associated  cryogenic  equipment  will  be  high  and  there  is  the  additional  problem 
of  mdiating  away  the  excess  heat^. 

To  overcome  the  weight  disadvantages  of  the  large  compressors  required  for 

gas  MPD  cycles,  two  phase  liquid-vapour  systems  are  currently  under  consideration 

in  which  the  reactor  is  cooled  by  a  liquid  metal  which  then  mixes  with  metal 

1 1 

vapour  and  expands  the  latter  into  the  MPD  channel  ,  It  seems  probable  that 
the  overall  efficiency  of  this  system  will  be  lower  than  that  of  the  usual  MPD 
gas  cycle,  although  tne  power  to  weight  ratio  may  be  similar  since  compressors 
are'  avoided. 

28.4  CONCLUSIONS 

For  power  levels  up  to  about  100  kWe,  solar  or  isotope  powered  systems  will 
probably  provide  the  necessary  power  for  space  vehicles;  power  levels  of  about 
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kWe  will  enable  early  interplanetary  missions  to  be  accomplished  using  electric 
propulsion  systems.  For  power  levels  ranging  1  ross  50  kWe  to  several  tens  of* 

MWe,  several  alternative  conversion  systems  based  on  nuclear  reactors  are  under 
consideration.  Of*  these,  alkali  aetai  R&nklne  cycle  turtoelectrio  systems  are 
receiving  the  most  attention  at  present,  although  many  problems  beset  the 
development  programmes.  Several  other  systems,  based  on  the  Bray  ton  gas  cycles, 
show  promise  of  greater  efficiencies  and  reliabilities,  thermionic  and  MPD  systems 
being  particularly  interesting  above  one  megawatt.  It  is  possible  that  the 
ultimate  choice  of  system  will  be  based  on  reliability  (particularly  for  manned 
missions)  and  the  use  OJ  hybrid  systems  involving  two,  or  more,  of  the  conversion 
methods  mentioned  above,  should  not  be  discounted, 
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lar  Brayton  cycle 
stem 
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m  temperature 
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formed  aluminium,  epoxy  fibreglass 
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RADIATOR  AREAS  AND  MAXIMUM  WORKING  TEMPERATURES 


FOR  VARIOUS  CYOIiES 
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same  area,  °K 
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15% 
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Brayton  (a) 

0.60 
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84 
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Tri -cycle 

wmm 
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KB 

46 
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(Grey  and  Williams"' 


Thermionic  cell  efficiency 
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Out-of-pile  gas  (radiation  cooled) 
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7 
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2 

Liquid  metal  turboelectric 

1.3 

1.3 

1.3 
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2 
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CHAPTER  29 


A  NON-EQUILIBRIUM  ELECTRON  HOPE  FOR  KILO WATT-HANG- 


MFD  SPACE  POWER* 


by 


l.R,  MeNab 


29 . 1  INTRODUCTION 

The  demand  for  electrical  power  sources  for  terrestrial  and  space  purposes 
is  increasing  rapidly  and  many  techni  ues  are  presently  under  investigate  on  in 
efforts  to  increase  the  efficiency  of  conventional  equipment  and  to  provide 
alternative  systems.  The  problems  associated  with  space  power  systems  are 
severe  so  far  as  the  power-to~weight  ratio,  reliability,  and  operating  conditions 
are  concerned;  it  is  these  systems  which  will  be  considered  here. 

Present  and  postulated  space  power  systems  include:  solar  cells,  solar 
collectors  (Brayton  or  Rankine  cycle  and  thermoelectric),  nuclear  reactor  - 
thermoelectric,  nuclear  reactor  -  turbogenerator,  nuclear  reactor  -  condensing 
fluid,  radioisotopes,  fuel  cells,  and  magnetoplasmadynamic  (MPD)  generators, 
i/hile  several  of  these  systems  (for  example  solar  cells)  are  already  well  proven, 
others  hold  promise  of  generating  much  larger  amounts  of  power  with  higher 
efficiencies:  one  such  system  is  MPD  generation, 

MPD  generators  may  be  divided  into  two  main  categories:  combustion  and 
nuclear  reactor  powered.  Combustion  devices  are  usually  open  cycle  and  nuclear 
devices  closed  cycle.  Direct  current  generating  combustion  devices  are  the  most 
advanced  at  the  present  time  so  far  as  power  output  is  concerned.  Substantial 
amounts  of  power  have  been  generated  by  several  systems,  for  example,  the  AVCG 
Mark  2  generator  has  given  an  output  of  over  1  MWe  with  an  efficiency  of  about 
7  per  cent.  Larger  devices  are  expected  to  have  a  higher  efficiency  and  it 
is  likely  that  land-based  MPD  stations  generating  several  tens  of  megawatts 
output  will  be  operational  within  the  next  few  years.  However,  for  a  space 
mission  of  duration  greater  than  about  one  week,  it  is  unlikely  that  open  cycle 
combustion  systems  will  be  attractive  in  view  of  the  large  fuel  supplies 
required.  At  present  the  prospects  of  developing  closed  cycle  combustion 
devices  are  remote, 
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Closed  eyw  nuclear**  failed  MPD  systems ,  while  conserving  the  working 
fluid,  suffer  from  the  temperature  limitations  of  present  and  foreseeable 
reactors,  thus  limiting  the  Saha  electrical  conductivity  and  consequently  the 
maximum  attainable  power  output.  Development  of  nuclear  reactors  to  yield 
high  working  temperatures,  while  resulting  in  a  higher  electrical  conductivity 
and  power  output,  will  introduce  more  severe  materials  problems  which  are  only 
capable  of  long  term  solution.  The  solution  to  this  paradox  is  to  increase 
the  electrical  conductivity  of  the  working  fluid  in  the  generator  without 
increasing  the  gas  temperature.  High  power  outputs  could  then  be  obtained 
without  severe  materials  problems. 

Several  methods  of  obtaining  this  extra-thermal  ionization  in  MFD 
generators  have  been  suggested  and  a  few  have  been  tried  experimentally. 

a 

Karlovitz  and  Halses  ,  who  started  work  on  MPD  in  1938,  investigated 
several  modes  of  extra- thermal  ionisation  in  their  natural  gas- fuelled 
combustion  generator.  Most  attention  was  paid  to  ionisation  by  high  voltage 
electron  beams,  but  glow  discharges,  arcs  and  r.f.  were  also  tried.  These 
attempts  were  unsuccessful  mainly  because  of  the  effects  of  dissociative 
recombination  which  rapidly  removed  the  electrons  created  by  the  ionisation 
process.  This  fast  recombination  mechanism  will  not  be  as  severe  in  a 
monatomic  gas. 

Kerrebrock  has  Tiggested  and  investigated,  theoretically  and  experimentally, 
magnetically- induced  ionization  in  MPD  generators.  In  these  processes  the 
electrons  in  the  plasma  gain  energy  from  the  electric  field  induced  by  the  magnetic 
field  through  which  they  pass.  The  energy  gained  in  this  manner  is  lost  through 
collisions  with  the  gas  particles,  but  if  the  parent  gas  is  monatomic,  there  are 
no  molecular  energy  leases  and  an  elevated  electron  temperature  can  result. 

Since  electron-atom  collisions  are  the  dominant  ionisation  mechanism  under  MPD 
generator  conditions,  the  electrical  conductivity  is  considerably  enhanced.  'Hie 
success  of  this  method  depends,  amongst  other  factors,  on  the  degree  of  purity  of 
the  gases  present,  as  may  be  seen  from  the  following  expression^: 

T.  .  k1**1 

_  ^  B  2 

36(1+a  > 

© 

where  8  is  a  factor  accounting  for  the  effects  of  inelastic  collisions  (other  than 
ionising  ones)  and  non -Maxwellian  electron  energy  distributions.  In  general,  8 
is  a  function  of  temperature  but  typical  values  for  MPD  generators  are  8  *  1  for 
monatomic  gases  and  8  ®  10^  for  molecular  gases.  The  theory  has  been  verified 
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experimentally  by  Kerrebrock?  with  an  electric  field  and  tentatively  by  Robben 

with  both  electric  and  magnetic  fields.  Devices  under  construction  to 

investigate  and  utilize  this  effect  include  those  at  IED  and  Martin-Marietta, 

5 

Recently  it  has  been  reported  that  the  device  at  Martin  has  proved  extra- 
thermal  ionization,  but  no  details  are  available.  Rosah  has  shown  that  the 
production  of  extra-thermal  ionization  by  this  means  is  quite  sensitive  to  the 
degree  of  uniformity  of  the  gas;  it  is  possible  therefore  that  the  effect  may 
prove  more  difficult  to  achieve  experimentally  than  originally  anticipated. 

7 

Sternglass  et  al  have  discussed  theoretically  extra -thermal  ionization 
for  closed  cycle  nuclear  powered  MPD  systems  by  electron  beams,  by  /S  rays  from 
fission  products  and  by  d.c.  arcs.  Of  these  methods  electron  beams  appear  the 
most  promising  in  view  of  the  high  efficiency  of  generation  of  beams,  continuous 
application,  ease  adjustment  of  energy  spectrum  and  lack  of  contamination 
caused  by  their  use.  The  main  disadvantage  appears  to  be  the  technical 
complexity  of  the  system.  Fission  product  ionization  is  less  efficient  and 
leads  to  radiation  hazards,  although  this  may  not  be  a  severe  limitation  for 
space  systems,  D.G.  arcs,  while  giving  high  conductivities,  cause  severe 
heating  of  the  electrodes  (which  must  be  cooled)  in  addition  to  Joule  heating 
of  the  gas  stream. 

c 

A  novel  extra-thermal  ionization  method  suggested  by  Maitland  is  to 
utilize  photoionization  caused  by  a  high  energy  light  beam.  Initially  a  laser 
was  suggested  as  the  light  source,  although  this  has  a  low  basic  efficiency  but 
the  recent  development  of  high  intensity  ultra-violet  lamps  may  make  this  a  more 
interesting  suggestion. 

9 

Gourdine  is  presently  investigating  the  production  of  non-equilibrium 
plasmas  by  r.f.  discharges.  The  main  difficulties  are  associated  with  the  low 
density  of  the  plasma,  the  difficulty  of  coupling  into  it,  and  diffusion  losses. 
The  technique  has  the  advantage  of  providing  high  flow  velocity  and  hi^ri 
conductivity  independently  -  a  useful  condition  for  high  efficiency  MPD 
conversion,  bearing  in  mind  the  temperature  limitations  associated  with 
materials. 

29.2  PROPOSED  NEW  MODE 

Although,  as  is  apparent  from  the  above,  irany  methods  of  producing  extra- 
thermal  ionization  in  closed  cycle  MPD  generator  systems  are  being  actively 
examined  at  present,  to  date  no  MPD  generator  producing  substantial  amounts  of 
power  has  yet  utilized  such  an  effect.  It  is  the  purpose  of  this  paper  to 
examine  another  mode  of  operation  which  may  permit  high  power  outputs  to  be 


29.5 


obtained  from  MPD  generators.  This  method,  suggested  independently  by  Lindley 
and  iischearoeder  and  Daiber  '  ,  consists  of  rapidly  expanding  a  high  enthalpy  gas 
through  a  large  area  ratio  nozzle,  so  that  the  rate  of  fall  of  gas  temperature  is 
eery  large.  When  the  temperature  falls  rapidly  in  a  gas  many  internal  adjust¬ 
ments  must  be  made  by  molecular  collisions.  For  example,  the  energy  of 
molecular  vibrations  must  be  reduced,  a  new  balance  must  be  found  between  atoms 
and  molecules,  chemical  reactions  must  take  place  at  a  different  rate,  and  ions 
and  electrons  must  recombine  to  form  neutral  particles.  All  these  adjustments 
require  a  large  number  of  collisions  between  particles  before  a  new  equilibrium 
is  established.  If  the  temperature  changes  faster  than  the  equilibrium  is 
established  then  non-equilibrium  flow  results. 

Many  studies,  both  theoretical^  ^  and  experimental^  4“",  on  non-equilibrium 
flows  have  been  published  recently,  most  are  concerned  primarily  with  molecular 
gases  or  gas  mixtures  where  vibrational  and  dissociative  non-equilibrium  are  the 
dominant  factors.  This  interest  arises  because  the  main  domains  in  which  such 
flows  are  encountered  relate  to  real  or  simulated  (shock  tube  and  wind  tunnel)  high 
altitude  hypersonic  and  missile  re-entry  conditions  in  air  and  si mi l«r  conditions 
in  combustion  gases.  In  general  the  theoretical  solution  of  such  problems  is 
complicated  since,  in  addition  to  solving  the  quasi-one-dimensional  flow  equations, 
the  population  densities  of  the  (many)  species  present  must  be  determined  from  the 
rate  equations  governing  the  appropriate  reactions.  In  many  cases  the  required 
rate  coefficients  are  not  fully  tabulated.  Solution  of  the  flow  equations  is 
generally  complicated  since  the  changes  of  energy  during  the  relaxation  phenomena 
usually  form  a  eignificant  fraction  of  the  total  enthalpy  of  the  gas. 

The  flow  of  a  monatomic  recombining  gas  is  more  amenable  to  analysis  since 
valid  simplifying  assumptions  can  be  made  concerning  the  number  of  species  present 
and,  consequently,  fewer  rate  equations  are  required  to  characterize  the  flow. 

For  example,  the  recombination  mechanism  may  be  assumed  to  be  controlled  by  only 
one  excited  state.  One  significant  difference  between  a  recombining  atomic  gas 
and  a  relaxing  molecular  gas  is  that  the  electron  temperature  may  be  considerably 
higher  than  the  temperature  of  the  heavy  particles  in  the  atomic  gas.  This  leads 
to  reduced  recombination,  amongst  other  effects. 

2*f 

Bray  and  Wilson  (following  the  methods  outlined  in  references  25  and  26) 
have  described  the  flow  of  an  atomic  recombining  gas  in  terms  of  the  equations 
of  mass,  momentum  and  energyT  the  thermodynamic  equations  of  state  and  enthalpy, 
the  rate  equation  for  production  and  loss  of  electrons  and  the  electron 
temperature  equation.  Even  such  a  simple  model  is  difficult  to  solve  and  again 
requires  detailed  knowledge  of  the  rate  coefficients. 
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To  examine  the  consequences  of  non-equilibrium  nozzle  flow  on  the  design 
a«d  power  output  of  an  MPD  generator  a  »auch  simplified  model  has  been  constructed 
which  retains  the  essential  features  of  a  non-equilibrium  flow* 


29-3  THEORETICAL  MODFX  AND  RESULTS 


In  concepts  of  closed  cycle  nuclear  reactor  MPD  generators  the  working 
fluid  is  usually  selected  to  be  an  inert  gas  seeded  with  small  amounts  of  an 
alkali  metal.  This  may  not  be  the  final  form  of  an  !PD  space  power  supply  (for 
example,  a  completely  condensing  cycle  may  be  superior  on  a  power- to-weight 
basis)  but  such  a  system  will  probably  be  the  first  to  yield  experimental  data 
fluid  it  is  this  type  of  generator  which  will  be  considered  here.  The  assumption 
of  small  fractional  seeding,  which  is  generally  valid,  enables  considerable 
simplification  to  be  made  in  the  analysis  since  the  energy  involved  in  the 
ionization  and  recombination  processes  in  the  seed  element  form  only  a  small 
fraction  of  the  total  enthalpy  of  the  working  fluid  and  thus  may  be  neglected. 
Consequently,  neglecting  boundary  layers  and  shock  phenomena  the  quasi-one- 
dimensional  flow  equations  for  the  parent  gas  only  can  be  assumed  to  approximate 
the  real  flow  system. 

To  simplify  the  model  further  the  isentropic  flow  of  a  perfect  gas  has  been 
considered  so  that  the  relation  between  area  ratio  and  Mach  number  in  the  nozzle 
is 
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which,  for  t  =  5/3,  is 
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....  (29.1) 


To  avoid  discontinuities  in  the  theoretical  flow  caused  by  a  discontinuous 
throat  area  the  nozzle  shape  has  been  assumed  hyperbolic,  of  the  form 


A  =  A* 


....  (29.2) 


where  x  is  the  distance  along  the  nozzle  axis  measured  from  the  throat  and  where 
the  half-angle  of  the  nozzle  at  infinite  distance  from  the  throat  is  given  by 

tan  d  =  . 

Dsing  the  non-dimensional  quantities  k  -  ^ K*  and  x  s  fyi*/fA*,  and  combining 
equations  (29.1)  and  (29.2)  yields  an  expression  for  Mach  number  in  teA-u3  of  the 
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dimensionless  distance  through  the  nozzle 
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*hich  has  the  form  shown  in  Fig,  29,1  The  temperature,  velocity  etc*  at  any 
point  in  the  nozzle  ay  then  be  obtained  from  the  usual  iaentropic  flow  relations. 


The  electrical  conductivity  at  any  point  in  the  nozzle  is  determined  from 
the  flow,  rate  and  electron  energy  equations.  In  view  of  the  difficulties 
associated  with  the  values  of  the  rate  coefficients  and  the  solution  of  these 
equations  it  is  assumed  that  the  electrical  conductivity  at  the  exit  of  the  nozzle 
is  related  to  the throat  value  by  an  exponential  function  of  the  form: 


V  *  "thp  *Jtp("  CV  •••• 

where  c  is  a  coefficient  depending  on  the  rate  coefficients  of  the  various 
processes  occurring  and  the  electron  temperature,  and  t^  is  the  time  of  flight 
through  the  nozzle. 


The  exit  conductivity  is  related  to  the  throat  conductivity  rather  than  the 
inlet  conductivity  since,  in  general,  the  gas  is  not  accelerated  sufficiently 
rapidly  in  the  subsonic  region  of  a  nozzle  to  cause  non-equilibrium  flow*  The 
throat  conductivity  may  be  obtained  from  the  inlet  parameters  through  the  use  of 
the  isentropic  relations,  which,  for  a  perfect  monatomic  gas,  are 
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....  (29.6) 


T0  and  pQ  being  the  stagnation  temperature  and  pressure  respectively* 
The  time  of  flight  between  two  points  (x^,  x^)  in  the  nozzle  is 
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For  the  isentropic  flow  of  a  perfect  gas  with  t  *  V3 
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Evaluation  of  the  integral  enables  the  time  of  flight  and  hence  conductivity 
to  be  calculated  at  any  point  for  given  inlet  conditions* 

The  power  output  per  unit  length  of  channel  for  an  MPD  generator  when 
Hall  effects  and  ion  slip  are  neglected  is 

P  =  erV^B^K(l-£)A 

If  the  MPD  generation  section  jcins  smoothly  on  to  the  exit  of  the  nozzle  this 
expression  gives  the  power  output  per  unit  length  at  the  inlet  of  the 
channel  when  <r,  V  and  A  are  the  exit  values  for  the  nozzle,  so  that 
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For  a  family  of  hyperbolic  nozzles  with  constant  expansion  angle  and  x^  *  -  $ex  » 
constant,  the  values  of  temperature ,  velocity,  etc.,  are  identical  at  each  value 
of  x  (for  the  same  inlet  conditions)  provided  the  mass  flow  i®  adjusted 
appropriately. 


Examining  such  a  family  of  nozzles  for  constant  magnetic  field  and  load 
factor  shows  that  if  the  electrical  conductivity  were  determined  only  by  the  Saha 
(equilibrium)  equation,  the  exit  conductivity  would  be  the  same  for  each  nozzle, 
and  thus  the  power  output  per  unit  length  of  channel  would  increase  linearly  with 
cross  sectional  area. 


If,  however,  extra-thermal  conductivity  occurs  owing  to  the  non-equilibrium 
effects  of  a  rapid  expansion  this  conclusion  will  no  longer  hold.  The  conductivity 
will  decrease  as  the  nozzle  dimensions  increase,  owing  to  the  increased  time-of- 
flight,  while  the  cross  sectional  area  will  increase;  consequently,  a  optimum 
nozzle  size  can  be  expected  for  any  given  set  of  inlet  conditions. 


To  examine  this  optimum  nozzle  size  for  the  family  of  nozzles  with  = 
^5L,.  =  constant,  equation  (29*8)  is  used;  V„  is  constant,  a  is  determined  from 

BA 

equations (29. *0  to  (29.7)  and  A  from 
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Combining  these  equations  yields 
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The  variation  of  I  with  *2  *  ^0X  (that  is,  from  the  nozzle  throat  to  the  exit)  is 
shown  in  Fig.  29*2. 


Differentiation  of  equation  (29.10)  yields  the  result  that  the  maximum  value 
of  P  occurs  when 


....  (29.14) 


The  value  of  k*  for  a  hyperbolic  nozzle  with  $  _  *  2  (that  is,  A  *  5A*), 
op  v  ex  ex 

0  *  10°  and  Cp  tar  pure  helium  is  shown  in  Fig.  29*3  as  a  function  of  Tq  and  c. 
The  values  of  To  used,  1000°  -  2500°K,  lie  in  the  range  of  interest  for  MPD 
generators.  There  is  little  information  available  at  present  which  permits 
accurate  evaluation  of  the  coefficient  c  in  the  above  expressions.  Tentative 
calculations  carried  out  at  IRD  indicate  that,  for  a  helium  plasma  seeded  with 
one  atomic  percent  of  cesium,  c  will  be  m  10^  /sec;  thus  the  range  of  c  shown  in 
Fig*  29*3  should  adequately  cover  the  error  in  the  independent  calculations. 

When  the  value  of  x _  is  altered  the  value  of  A*  ,  is  shifted  owing  to 

ex  opt 

the  change  in  time  of  flight  through  the  nozzle  which  causes  a  decrease  in 

conductivity  as  £  is  increased;  this  effect  is  partially  offset  by  the 

increase  in  exit  velocity*  The  effect  is  shown  in  Fig,  29*4  for  T  *  2000°K 

and  c  *  lO'Vsec,  where  A*  .  Is  seen  to  decrease  as  £  increases. 

cot  ex 

The  power  output  per  unit  length  of  the  MPD  channel  when  A*  is  optimized 

is 
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f or  helium  seeded  with  one  atomic  percent  cesium  is  -  32  k*W/a,  The  actual  power 
output  in  a  metre  long  channel  will  be  less  than  this  since  power  extraction  will 
cause  a  rapid  drop  in  the  total  enthalpy  of  the  gas.  This  will  reduce  the 
Telocity, leading  to  a  reduction  in  power  output  through  the  V2  term  and  also  a 
reduction  in  conductivity  because  of  increased  recombination.  Nevertheless 
this  calculation  does  show  that  the  power  output  of  an  MPD  generator  using  a 
non-equi librium  nozzle  expansion  will  lie  in  a  range  of  interest  for  space  power 
purposes. 

Z9,k  DISCUSSION  AND  CONCLUSIONS 

Although  nany  methods  of  producing  extra -thermal  ionization  in  closed  cycle 
MPD  generators  have  been  suggested  and  tried  experimentally  none  has  yet  been 
utilized  in  a  device  producing  substantial  amounts  of  power*  An  alternative 
method  of  gaining  a  high  electrical  conductivity  is  to  use  the  (already  proven) 
methods  of  producing  a  non-equilibrium  flow.  The  non-equilibrium  expansion 
through  a  hyperbolic  nozzle  of  a  plasma  with  small  fractional  seeding  has  been 
examined  here  using  a  simple  theoretical  model.  The  investigation  shows  that 
the  nozzle  exit  conductivity  will  be  high  compared  with  the  equilibrium  value* 

If  such  a  nozzle  immediately  precedes  an  MPD  generator  section  this  r.j  gh 
conductivity  enables  High  power  outputs  to  be  achieved  with  relatively  small 
devices,  for  example,  for  the  case  quoted  in  the  previous  section,  an  initial 
power  output  per  unit  channel  length  of  ^  32  kV/m  for  a  channel  of  7  cm2.  A 
power  source  producing  such  an  output  is  certainly  of  interest  for  space 
purposes;  more  refined  calculations  including  the  loss  mechanisms  neglected 
here  (heat  transfer,  wall  friction,  boundary  layers,  shock  waves,  Hall  effects, 
ion  slip)  will  be  worthwhile. 

The  optimum  throat  area  for  a  hyperbolic  nozzle  expansion  has  been 
evaluated  and  shown  to  depend  on  the  inverse  square  of  the  coefficient  c. 

This  coefficient  depends  on  the  electron  temperature  throughout  the  expansion 
and  the  rate  coefficients  for  all  the  atomic  processes  occurring.  In  order 
that  it  might  be  accurately  evaluated  the  complete  solution  of  the  non-equilibrium 
flow  must  be  found.  While  there  is  at  present  some  doubt  concerning  the  rate 
coefficients,  independent  calculations  carried  out  at  JRD  indicate  that,  for  a 
helium- 1  a/o  cesiua  plasma,  c  will  be  10^/sec  for  a  nozzle  inlet  stagnation 
temperature  of  2000°X.  Under  these  conditions  the  optimum  throat  area  is 

•A 

1  *5  c/,  a  size  which  falls  into  the  range  of  interest  for  space  purposes. 

It  may  be  concluded  that  the  effect  outlined  here  is  certainly  of  interest 
for  space  power  sources  and  is  worthy  of  a  more  accurate  evaluation  than  so  far 
performed.  It  should  be  noted  that  other  methods  of  producing  extra-thermal 
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ionics 3 on  in  the  MPD  generator  channel  (for  example,  electron  beam's, 

magnetically- induced  ionization,  etc.)  are  not  precluded  by  the  use  of  a 

non-equilibrium  flow. 
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OPTIMUM  NOZZLE  THROAT  AREA  AS  A  FUNCTION  OF  THE 
COEFFICIENT  C  AND  NOZZLE  STAGNATION  TEMPERATURE 
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OPTIMUM  NOZZLE  THROAT  AREA  AS  A  FUNCTION  OF  THE 
OIMENSONLESS  N02ZLE  EXIT  DISTANCE 


NOMENCLATURE  (Chapter  29) 
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area  of  nozzle  or  channel 
magnetic  field  strength 

coefficient  in  electrical  conductivity  equation  (29**0 

coefficient  defined  by  equation  (29*10) 

specific  heat  at  constant  pressure 

integral  defined  by  equation  (29*11) 

transverse  MPD  generator  load  factor 

axial  MPD  generator  load  factor 

nozzle  constant 

Mach  number 

pressure 

stagnation  pressure 

power  output  per  unit  channel  length 

characteristic  power  output 

time  of  flight  in  nozzle 

gas  femperature 

electron  tempera ture 

stagnation  temperature 

flow  velocity 

axial  distance 

electronic  Hall  coefficient 

ionic  Hall  coefficient 

ratio  of  specific  heats 

factor  accounting  for  inelastic  energy  losses 
hypersonic  nozzle  angle 
electrical  conductivity 


Subscripts 

ex  -  nozzle  exit 
in  -  nozzle  inlet 
opt  -  optimum  value 
thr  -  nozzle  throat 


Superscripts 

•  -  value  at  sonic  speed 

*  -  non-dimensionalised  q  entity 


CHAPTER  30 

NQN-EWUILIBRIUH  PIASMA5 


Dy 


I.R.  MoNab 


30.1  INTRODUCTION 

For  a  plasma  in  thermal  equilibrium  (that  is,  in  a  steady  state  of 
equilibrium  with  its  surroundings  and  having  no  interaction  with  them)  the  Saha 
equation  ,  wniuh  is  derived  on  a  purely  thermodynamic  basis,  serves  to  give  the 
particle  concentrations  of  the  various  constituents  present  (atoms  in  ground  and 
excited  states  electrons  and  ions).  If  however,  for  any  reason,  the  plasma  is 
disturbed  from  this  equilibrium  state  the  particle  concentrations  cannot  be 
obtained  from  the  Saha  equation  but  must  be  found  by  investigating  the  appropriate 
particle  creation  and  annihilation  mechanisms.  This  situation  applies  if  a 
plasma  is  momentarily  disturbed  from  the  thermal  equilibrium  state  and  then 
relaxes  back  to  it,  or  if  the  plasma  as  maintained  in  the  non- thermal  equilibrium 
state  by  a  continuous  disturbance.  In  the  first  case  the  rate  at  which  the 
plasma  returns  to  thermal  equilibrium,  and  the  particle  concentrations  during  this 
period,  must  be  obtained  from  the  appropriate  rate  equations  for  the  various 
mechanisms  involved;  in  the  second  case  the  non-thermal  equilibrium  particle 
concentrations  are  obtained  by  setting  the  appropriate  rates  equal  to  zero. 

At  present  it  appears  likely  that  if  closed  cycle  gas  nuclear  MFD  systems 
are  to  prove  economically  feasible,  non-thermal  equilibrium  plasmas  must  be  used 
in  some  form.  Two  typical  non-equilibrium  plasma  conditions  presently  under 
ligation  a*"  I R D  are!  non— equilibrium  f^,ow  during  a  rap-id  expansion,  and 
photoionization-induced  non -equilibrium.  In  both  these  cases  evaluation  of, 
for  example,  the  electron  concentration,  reauires  solution  of  f'! 


ne  appropriate 


rate  equations. 


To  study  recombination  and  ionization  mechanisms  it  is  necessary  to  consider 
tne  rates  at  which  all  these  processes  can  occur  in  all  the  energy  levels  of  an 
atom.  While  in  principle  this  requires  an  infinite  set  of  simultaneous  differential 
rate  equations  it  can  be  shown  in  practice  that  convergence  can  be  obtained  with 
a  finite  set.  Thus  Bates  et  ai  used  sets  of  about  twenty  simultaneous  rate 
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tna  calculation  of  recombination  and  ioni 


i-aiion  coefficients . 


equations 

A*«  8a  erna  */e  approach,  particularly  useful  as  many  of  the  required  rate 
coefficients  are  not  accurately  known,  is  to  simulate  the  many  energy  level 
system  of  an  atom  by  a  simpler  model.  For  example,  BonDaniel  and  Tamor  have 
investigated  a  four  level  model  of  the  cesium  atom  containing  the  ground  3tate, 
two  excited  States  and  the  continuum.  While  this  model  is  not  completely 
accurate,  in  many  respects  it  closely  resembles  the  many  level  system. 

Three  separate  lines  of  investigation  have  been  conducted  at  IRD  fit  must 
be  stressed  that  these  are  merely  of  a  preliminary  nature),  these  are; 

(i)  a  three  level  model  of  a  cesium  plasma  in  which  only  electronic  and  radiative 
processes  are  considered; 

(ii)  a  many  level  system  of  a  cesium  plasma  in  which  only  electronic  and  radiative 
processes  are  considered;  and 

(iii)  a  three  level  model  of  a  cesium-helium  plasma  in  which  atomic  processes  are 
included  in  addition  to  electronic  and  radiative  processes. 

The  aims  of  these  three  investigations  are  respectively: 


U) 


to  evaluate  the  deviation  from  Saha’s  equation  in  an  equilibrium 
^ non- thermal)  piasma  at  iow  electron  concentrations,  with  and  without 
resonance  trapping; 

(ii)  to  investigate  a  more  accurate  model  of  a  cesium  plasma  than  (i) 

and  to  consider  its  application  to  a  plasma  not  in  equilibrium  (thermal 
cr  otherwise);  and 

^iii)  to  investigate  the  influence  of  atomic  processes  ca  recombination  and 
ionization  in  plasma  mixtures  where  the  fraction  of  inert  gas  atoms 
is  high. 

Only  the  first  of  these  investigations  is  reported  here  since  the  remainder 
have  not  reached  the  stage  where  definite  conclusions  can  be  formed, 

30.2  THREE  LEVEL  MODEL 

For  a  three  level  model  for  pure  cesium  (ground  state,  excited  state, 
continuum)  in  which  electronic  and  radiative  processes  only  are  considered,  the 
ionization  and  recombination  mechanisms  investigated  are; 
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where  A  and  B  are  the  collisional  and  radiative  ionisation  and  ^combination 
rate  coefficients  respectively. 

The  equations  giving  the  rate  of  change  of  the  partible  densities  in  each 
of  the  three  levels  are: 

=  ~A(1 ,2)n(l  )n(c)  +  A(2 ,1  )n(2)n( c)  -  A(l ,c)n(l )n(c)  +  ACcpO^c)3 
-B(l,2)n(l)  +  B(2,l)n(2)  -  B(l,c)n(l)  +  B(e,l)n(c)2  .....  (30.7) 

~  A(1 ,2)n(l  )n(c)  -  A(2,1  )n(2  )n(c)  -  A(2  ,c)n(2)r.(c)  +  A(ct2)^n(c)^ 
+B(l,2)n(l)  -  B(2,l)n(2)  -  B(2,c)n(2)  +  B(c,2)n(c)2  .....  (30,8) 

:  A(l ,c)n(l )n(c)  -  A(c,l)1n(c)3  +  A(2,c)n(2)n(c)  -  A(c,2)1n(c)3 
+B(l,c)n(l)  -  B(c,l)n(c)2  +  B(2,c)n(2)  -  B(c,2)n(c)2  .....  (30.9) 

and  the  conservation  equation  is: 

N  =  n(l )  +  n(2)  4-  n(c)  _ _ (30.10) 

(N.B.  In  equations  (30.7),  (30.8)  and  (30.9).  A(c,l)  and  A(c,2)^  are  not 
functions  of  n(c),) 

In  equilibrium  the  rates  of  change  of  particle  concentration  with  time  are 
zero.  Making  dn(2)/dt  and  4u(c)/dt  equal  to  zero  in  equations  (30.8)  and  (30,9) 
n(2)  may  be  eliminated  (after  some  manipulation)  and  a  relation  between  n(l )  and 
n(c)  obtained: 


du(c) 

dt 


i(l)  =  - 


n(c)  [_n(  c )  2K^  +  n^)^  K^J 


n(c;Kj  +  n(c)K^  + 

where  I ^  =  -A(2 ,c)A(c,1  -  A(2 ,1  )A(c,1  J1  -  A(2,l  MCc^)1 

K2  =  -A(2,c)b(o,1  )  -  A(2,l)B(c,l)  -  A(2si)B(c.2)  -  B(2,c)A(c ,1 J1 
-B(2,l)A(c,l)1  -  B(2,l)A(c,2)1 
Kj  =  A(l,2)A(2,c)  +  A(2,c)A(l,c)  -  A(2,l)A(l,c) 

^  =  -B(2,c)B(c,l)  *  B(2,l)B(c,l)  -  B(2fl)E(cs2) 


(30.11 ) 
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